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t .EftgBAL. . 

Th« program for the DC8~D3 flight planning has been slightly modified 
for the three dimensional flight planning for DCIO aircrafts. 

Several test runs of the modified program over the North Atlantic 
and North America have been made for verifying the program 
in this project. 

Reference is made to pages and formulas in 'Optimal Track Selection 
and 3**Dimensional Flight Planning' uiritten by H. M. De Jong* No. 93. 

r^DlFlCATtON 

While Or. Oe Jong used geopotential height and temperature in his 
program as a meteorological data* the modified program uses wind 
(direction and speed) and temperature^eceived from National Weather 
Service. 

A scanning program has been written to collect re<iuired weather 
information from the raw data received in a packed decimal format. 

Two sets of weather data, the 12~hour forecast and 24>hour forecast 
based on COCO 6NT« will be used for dynamic processes in our testruns. 
In order to save computing time only the weather data of the North 
Atlantic and North America is previously stored in a PCF file and then 
scanned onif by one. 

The element name (for weather data) of the PCF file is 'NWS' in our 
test ri»nv.. 

Wind data is supplied in the form of grid point values of the vector 
wind 

Vector wind is decomposed into x and y components to be applied to 
interpolation schemes. 
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Th»n ii prepares spot values by interpolation uiith respect to the 
values of the uiind components in the surrounding vertices of grid. 
See Fig. i> Fig. 2 and Fig. 3. 


£ SLUL BAJLnJL&- 

Mhile Dr. De Jong used weather data depicted from the grid points in 
a Cartesian grid superposed on a polar stereographic chart projection 
with standard parallel at 60 degree N> the modified program uses 
weather data depicted from thd same grid points that National 
Weather Service uses. 

See 'Aviation Digital Forecasts Program' issued by National Weather 
Service in Aug» 1978. 

A blockette and grid points: 


* 7. 5 degree 
» « 

e * 

«•———» 0 degree 
5 0 


Each gf-id point contains its Ouin weather data, 
a. In order to compute a great circle distance on earth between two 
points the following equation is used: 

DiSTs60 arccos( sin LS sin LD ■*’ cos LS cos LD cos < RD RS> ) 
where !..S=lati tude of starts 

LO-loti tude (>f destination* 

HS=longitude of start* 

KD-*lonQitude of destination* 

DiST~great cicle distance between start and destination. 

(See equation 5.20. P. 64). 
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b. To curaputo a true course it uses the follouing equation: 

TC*arccos< (sin LD - sin LS cos(0l8T/60> ) / ( sin (DlST/60) cos L8> 
< sinCRD-RS) > 0 then (360*TC) > is used 
inhere rC^ true course in degree. 

c. Computation of a drift angle uses the following equation: 

OA«arcsin(MV sin(UO - TO / TA8) 
where DA»drift angle* 

UD=wind direction (true) 

WV^wind velocity* 

TAS^true air speed. 

d. Headings TC * DA (in degree). 

d. PJLVjtoJLm^ric^e data £o£. DCIO . 

a. Specific range table as a function of weight and latitude. 

This is an important economic index which is the air distance 
covered by a turbo-jet aircraft per unit of fuel consumption. 

This data is almost independent of temperature and depends on the 
aircraft weight only at a prefixed flight level. 

See table 1> P. IIA and Fig. 4 in this documentation. 

b. Climb time table as a function of temperature and weight. 

See Fig. 5 

c. Climb distance table as a function of temperature and weight. 

See Fig 6 

d. Climb fuel table as a function of tempveature and weight. 

See Fig 7 

e. Max. weight table as a function of' flight and temperature. 

The maximally allowable weight depends on flight altitude and 
temperature deviation from standard. 

S«?e Fig 4. 
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f. Priram^ter rl and r2 ( Sat tablt 2* P. 116). 

Parametvr rl and r2 art derivtd from tht tptciflc rangt tablt. 

Hare rl and r2 danott tht eotfficitnte of a ttraight lint fit 
in a ( uitight* 1/tntry value) graph at constant prtscurt altitude. 
According to the specific rangt table for DCIO performance data* 
rl and r2 art as folloei; 
flight level rl r2 


310 2. 966 0. 0000722 

350 0. 85 0. 0000792 

390 >0. S3 0. 0000890 

g! Extra burns for step-up and step-down. 

While extra burn parameters for DC6-D3 were 100 kg/ 4000 ft for 
step-up and —80 kg/ 4000 ft for step— down* the new program uses 
700 kg/4000 ft for ste-p~up and -28—kg/4000 ft for step-down 
for the DClO flight planning. 

h. The equations (6.31) in p. 121 and (6.33) in p. 123 will be used 
in our test runs# since precisely matched parameters for DCIO are 
not .available now. 

(See lines 874-879 and 1305-1313 of the new program). 


-I Ini ri al InHUS. 

FLURs* sequence number. 

TAX1» t ixi fuel < kg). 

OHW-- z»ro fuel weight (kg). 

RESERVE- reserve fuel (kg). 

TUW- If -1. then computation starting in end point of flight. 
MAX row* max. lakeuff weight(kg). 

MAXL.W- m<^ix. landing wsight(kq). 


DATE* d.:>tr/moiith/ye<ar. 
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10* If -I- th»n tatt bound flight. 

TTT» time of departure. 

ROUTE» route indicator. 

ST* starting graph point number. 

8T1* destination number. 

Mote: These input parameters are read from card images of the 
PCF file in our testruns (element name * DClOF). 

But these mag be read directly from cards. 


r See output listings. 

NO. *graph point number. 

HEAD* heading in flight (degree). 

FL* flight level in 100 ft. 

Trip* off standard temperature in centigrade degrees. 
TAS* true air speed (n. m). 

UIND*uiind along the track ( '*' * tail> -* * head). 
DIST* distance of flight segment. 

ACCD* accumulated distance f loufn(n. <n). 

TIIib=tinte along the segment flouin. 

Acer* accumulated time (hour* min). 

BURN* fuel consumed in segment <kg>. 

WEIGHT* uieight (kg). 

TOC- top of climb. 

TOD* cop of descent. 
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&tt.Y.u,ma. fta ImsJl 6Li4J1s.- 

Sincii Dr. DC Jong alroodg provod that • solution in S-spaco would 
roward as comparod with a 2~spaca solution^ wa will not traat tha 
factorized 2*apaca optimalization casa and tha i>rascribad track case. 
Thus wa will concantrata on tast runs about pura S^'dimansional flight 
planning in frea spaca. 

While averaga CPU tima for ona flight plan with the modified program 
is approximatalg 120 seconds* on UNIVAC 1108* tha original program 
takas approximatalg llO seconds or more for ona flight plan on 
UNIVAC 1108. 

Howavari for more than ona flight plan it would require lass than two 
minutes par one flight plan with tha modified program« since they can 
share precalculated results for the weather information. 

Note: In our test runs* two different sets of weather data and 

performance data will be used. 

Waatner iiaca. <i) qeopotential height and temperature* 

(3) wind (spaed* oirection) and temperature. 

Performance data: 

• 1 ) for DC9-D3: 


standard Ci 

uise* 

ORIGINAL PAGE IS 
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.flar h 

O. U034 


taxi fuel 

1500 kg 


Zero fuel 

75000 kg 


reserves 

10000 kg 



max. takeoff wt. 142900 kg 

I'lax. landing wt. 93000 kg 


iZi for DClO: 


r.t.'tndard cruita* 


t eiiuiarity 3 % 


Mach 


0.82 


taxi 

fuel 

862 

k9 


.:ai '0 

fuel 

177356 

kg 



3629 

k0 


max. 

takaoff 

lift. 

259458 

kg 

max. 

landing 

wt. 

190964 

kg 


a. Tffkt run* of tha original program (for DC6*-03>. 

‘1>- Jjssl rjtn. 

(a) . Input paramacar*. 

start*-da«tination: 
uiaathar data: 
parformance data: 
flight level rl 

JIOOO 4. 94 

1^5000 0 54 

,.■>9000 ~0 AS 

i>Ktra nurn. 

(b) iJutput 
* Sl* Fig a. 


cri ion 

trip fuel 

cost 

time 

diftanca 

m in. 

co*>t 

36398 kg 

9914 

6h 5lm 

3r.0S n. m 

min. 

1 1 Ilia 

39931 kc 

10359 

6h 43m 

3205 n. ffl 


Now York *> Amttardam (east bound) 
gaopotan9inal hoigiit and tamparatur# 
DC8-03 

r2 

0 0000693 
0. 0001043 
0 000117 

«tap up lOO kg/4000 ft 
• tap douin -SO kg/4000 ft 
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* routes taken: same route 


1 1 5-1 13-108-103-97-% J -87-77-66-52-37-31-11-6-2-0 
<c). Rtfmarkt. 

The original program written in fiurrough ALGOL hat been 
converted for wee in UNXVAC 11<^. 

Test run retultt which come from the converted program are jutt 
identical to those obtained frum the original program# assuming 
use of same input parameters. 


< 2 > 

<a). Input parameters. 


start-destination: 


New York - Amsterdam (east bound 


weather data: 
performance data: 


zero wind and standard temperature 
DC8-03 


rl and r2: 


same as in test run #1 


extra burn: 


same as in test run #l 


(b>. Output. 

* See Fig. 9. 


C'i cer ion 


min. fuel 


trip fuel rost 


time 


distance 


38327 kg 10357 7h OSm 3l93 n. m 


a route taken: 


11 5- 1 1 1 08- 1 03-97-9 1 -87-76-64-50-36-2 1 - 1 1 -6-2-0 

tc>. Remarks. 


To run tPe program under zero wind and standard temperature 


IS a kind of minimal check- 


The optimal track should then consist of the 'optimal operational 
distanca track' in the graph as a good approximation of the great 
circle ciack between the end points. 
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Thtt r«suxt» of tott run #2 tnotot « «li(|htlM diffcrottb rout# «• 
compared wtth that obtalntd in th« tott run #1. 

Tho rtouito ohould b« unaffdctntf whan thr program i« run for a 
flight in oppooita dlroctlon hotiuaon tho tamo pointt (Moot 
bound vortut oatt bound). 8oo table 11* P. 128. 

b. Tott runt of tho modifiod program. 

<I> lA^LLua Fa- 
ta). Input paramotort. 

ttart-'dottinabien; Now York - Amttordam (oatt bound) 
woathor data: soro wind and ttandard tomporaturo 

porf orman<:o data: D€8~03 


flight level 

rl 

r2 

31O0U 

4. 94 

0. 0000693 

39O00 

0. 54 

0. 0001043 

39000 

-0. 69 

0. 000117 


extra burn ftop up 100 Vg/4000 ft 

•trp down >-80 fcg/4000 ft 

(b). Otit|>ut.. 

♦ ao*» Kitj. 10. 

criterion trip fuel cost timo distance 

min. fuel 38220 kg 10330 Tb 09m .a1u 8 n. m 

♦ route taken: 

1 Ib-i lj-i0ti-l03-V7-9l-87-76-a4- 50-36-21-11-6-2-0 
tc>. Remarks 

Lot ut Compart test run 03 with tost run 02. 

These test runs with same input parameters should show close 
results for the optimal operational distance trad as shown 
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in Fig 37. P. 110. 

Boch test run #2 and #3 have same route bctb/een New York and 
/Nmi^terdam. 


program 

criterion 

fuel 

cost 

time 

distance 

original 

min. fuel 

38327 

10357 

7h 09m 

3193 n. m 

modified 

min. fuel 

38229 

10339 

7h 09m 

3186 n. m 


Differences in results may be attributed to the use of different 
formulas for computing segment distance and true course, and 
round of errors. 

Thus we know that the modified program works correctly under 
this minimal check. 


<a>. Test nun #4. 


(a). Input parameters, 
start-destination: 
weather data; 
performance data 
rl and r2; 
extra burn; 


New York - Amsterdam (east bound) 

wind and temperature 

DC8-D3 


same 

as 

in 

test 

run 

#3 

same 

as 

in 

test 

run 

#3 


(b>. Output. 

'^St; Fiy. 11. 

Ciiterion trip fuel cost time distance 

mm. time 34273 kg 9019 5h 9Sm 3192 n. m 
# route taken; 


115-113-1 08- 1 03-97-9 1 -87-77-6S-5 1 -37-2 1 - 1 1 -6-2-0 


(c>. Ri.marr. s 

Let us compare test run #4 with test run #i. 
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theftc tsftt run\> ui^rff •x#cut9d unddr conditions having standard 
cruise* Mach«0. S034* New York Amsterdam (east hound)* OCS^DS 
performance data* and min. time criterion. 


test run 

weather 

fuel 

cost 

time 

distance 

testrun#! 

geopotential 

39931 

10359 

6h 43m 

3205 n. m 

testrun#4 

wind 

34327 

9015 

Sh 58m 

3192 n. m 


Differences in results are attributed to the use of different 
weather data set. 

The difference in flight time implies the test run #4 used 
weather data set with almost east bound wind along the track. 
It is reasonable to have faster flight in the test run #4 than 
in the test run el since as gou see in Fig. 11* an aircraft to 
fly mag experience onlg tail wind*~«long the planned track. 

(3). Tef.t run #3. 

(a). Input parameter^;. 

start-destination: New York - Amsterdam (east bound) 

ttte.jtlier data: wind and temperature 

performance data. OClO 

rl and r2: same as in test run #3 

extra burn. same as in test run 43 

(b>. Oi»tput. 

:^er Fig. 12. 

criterion trip fuel cost time distance 
min. time 4876? kg 11129 5h 52m 3192 n. m 




« ruute taken: 


1 15-113-106-103-97-91-87-77 

< 4 > 

<a>. Input parameters, 
start-destination: 
weather data: 
performance data: 
rl and r2: 
extra burn: 

<b). Output. 

* St*e fig 13. 

A strange symptom occurs: A 

a negative fuel burn between 
Thus we temporarily add one 
thi'a negative magnitude. 

Let us make another test run 
same conditions ( test run 


65-5 1-37-21-11 -6-2-0 


Amsterdam - New York (west bound) 
wind and temperature 
DC 10 

same as in test run #3 
same as in test run #3 

negative segment flight time and 
graph point number 0 and 1. 
tatement in line #884 to avoid 

with this added statement under 


'5>. Feri-t run #7. 

(a) . Input parameters: same as in test run #6. 

(b) . Output. 

* Sue Fig 14. 


cr 1 ter ion 

trip fuel 

cost 

time 

distance 

mill, time 

68447 kg 

15383 

7h 52m 

3283 n. 1 
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* raute taken; 

0- t -4-8- 1 7-30-44-60-75-87-92-98- 1 04- 1 09- 113- 1 I 5 

<c). Remark t. 

Let us compare test run #5 uith test run #7. 

These test runs are made under following conditions: 

DCIO performance data* wind and temperature and same start- 
destination 


bound 

trip fuel 

cost 

time 

distance 

east 

48789 kg 

11129 

5h S2ffl 

3192 n. m 

west 

68447 kg 

15383 

7h 52m 

3283 n. m 


The west bound flight takes 2 hours more in flight time as 
compared with that obtained in east bound flight. 

This is caused bg the strong east~-bound wind trend in weather 
data and longer flight distance in west bound route taken. 

<7). Test E.UJ1 #a ini 
(a). Input parameters. 

start-destination: New York - Amsterdam (east bound) 

weather data; wind and temperature 

performance data: DClO 

■Plight level rl r2 

31000 ft 2 , 366 0. 0000722 

35000 ft 0. 35 0. 0000792 

39000 ft -O. 83 0. 000089 

I'Xtra Durn; step up 700 kg/4000 ft 

step down -28 kg/4000 ft 


(b). Output. 

:-sec Fig. 15 and 16. 


-rXt. L 


I 

I 

I 

1 

L 




i4 

crit«rlon trip #u«l cost tiAo 
min. fuel 49992 kg 11320 Sh 9^ 

min. time 91129 kg 11480 9h 92m 

o routes taken: sarao route 
119-1 13-108-103-97-91-87-77-45-91-37-21-11-6-2-0 
<c>. Remarks. 

These test runs show that the min. fuel flight plan gains 
1177 kg of trip fuel and loses 1 minute of flight time as 
as compared with the min. time flight plan. 

The gains may look inconsiderable. 

But this is the result of one test run only. 

Here another comment is made for optimal step altitude profiles. 
The fuel and cost profiles have both the normal appearance: 
Approximation of a continuous crae-se climb with one or two 4000 
ft steps. However* the min. time flight is more irregular. 

This is caused by the response of the aircraft's speed to the 
vertical temperature gradients: The level of highest 

temperatures is sought. 

As you see temperature values in Fig. 3> some higher altitude 
levels have unexpectedly lower temperature values* for example* 
around graph point numbers 87*77*37*21*2 etc. 


010 ( 18110 . MI # 

ofpookqumi^ 


distance 
3192 n.ffl 
3192 n. m 


‘ S> Is.'j. Eiia 4U1 SJUL 4JLL 


(a). Input parameters, 
star t-destinat ion: 
weather data; 
performance data: 
flight level 


New York - Amsterdam (east bound) 
wind and temperature 
DC 10 

rl r2 


31000 ft 


2 . 966 


0. 0000722 


iS 


35000 
39000 ft 
extra burn: 


<b>. Output. 


0. 85 
0. 83 


0. 0000792 
0. 000089 

ttap up 700 kg/4000 ft 
step doum -*28 fcg/4000 ft 


OFMORQIIMJfV 


* Sae Fig. 17 and 18. 


criterion 

trip 

fuel 

cost 

time 

distance 

min. 

fuel 

49948 

k9 

11319 

5h 53m 

3192 n. ffl 

min. 

time 

51118 

kg 

11479 

5h 52m 

3192 n. m 


* route taken: 

1 15-U3-108-103-97-91-87-77-65-51-37-21-1 1-6-2-0 
(c). Remarks. 

We used slightly different rl« -rS— parameter sets to see hom 
they affect flight plannings in this test run. 

When uie see test runs* #8* #9> #10 and #11 tee understand that 
there is no significant difference in results under slightly 
different rl . r2 parameter sets. 


9) 


Lm 9iS- 


<a). Input parameters. 
start*‘de$tination: 
weather data: 
perf or-mance data; 
rl and r2: 
e<tra burn: 

(b). Output. 


Amsterdam ^ New York (west bound) 
wind and temperature 
DC 10 

same as in test run #10 
same as in test run #10 


eSee fig 19 






criterion 
min. fuel 


Xk 


trip fuel cost time 
69928 kg 19667 7h 98m 


«9iipwii»i .. Mill III iJMlimmBnUpiRPKP 

ommKL Hmm 
ormmqmm 

dietanee 
3293 n. ffi 


(e). Remark*. 

In this test run we understand that there are significant 
differences in results mainly because of the east bound wind 
tendency in the weather data set used. 


criterion 

bound 

fuel 

cost 

time 

distance 

min. 

fuel 

east 

49948 

11319 

9h 53m 

3192 

min. 

fuel 

west 

69928 

15667 

7h 58m 

3283 


( 10). Test Tjfn #13. 

<a). Input parameters: same as in test run #12 except the following: 
60 % of vector wind for mean climb vector* 

4 minutes for take-off and acceleration. 

<b>. Output. 

♦ See Pig- 20. 


criterion 

trip fuel 

cost 

time 

distance 

min. fuel 

70404 kg 

15769 

8h 01m 

3283 n. m 


erouce taken; 

0- 1 -4-8- 1 7-30-4 4-60-75-87-92-98- 1 04- 1 09- 113-115 
tc). Remarks 

We made some modifications on climb distance computation to avoid 
the negative magnitude as described in the test run #6. 

We took 60 % of vector wind instead of 79 % for the mean climb 
vector wind (See line #882-863 of the program listing). 

The climb time is usually reduced by a few minutes to account for 


take-off and acceleration. 


We took four minutes for it in 
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this test run < tu>o minutet taken in test run #i2>. 


Then let us see how modi# icetions effect the flight plannings. 


test run 

criterion 

fuel 

cost 

time 

distance 

#12 

min. 

fuel 

69928 

15667 

7h 5Sm 

32B3 n. ffl 

#13 

min. 

fuel 

70404 

13770 

8h Olffl 

3283 n. m 


It shows that the test run #13 requires slightly more cost and 
time than the test run #12« hut not significantly. 

Graph point #1 has Oist*iO« Time»l< and burn«200 in this test run 
instead of zeros. 

EXTRA COMMENTS E3B, Ibi£ MODIFIED efiQfiBM 
54-334 

Procedure NDATA scans weather data received from National 
Weather Service and obtains the required weather information 
for North Atlantic and North America area. 

(see the weather data format in Operational Manual issued on 
9/19/78 and Aviation Digital Forecasts Program issued in 
Aug. 1978> by National Weather Service). 

The weather data received in packed decimal format is supplied 
in the form of grid point values of the vector wind and 
tempera r.ure. 

« 72-85 

Procedure CHOP cuts a piece of weather information from the 
whole weather data for a specifiad grid point. 

* 89-1 

Procedure STOR prepares a two dimensional array (120x34) for 
the weather information and stores the chopped weather data 




xi WPOOR 




in th« array according to tho atandard praaaura levala* i« a« 

300 mbar» 290 mbar» 200 mbar etc. 

« 119-133 

Procadura AREA gats tha baadar information of tha maathar data* 
i*a* area numbar* blocfcatta numbar ate. 

« 136-i42 

Procadura FINDNEXT ebaefcs tha and of tha waathar data block. 

* 149-213 

Procadura GETOATA storas tha waathar data of aach grid point 
in the prepared array according to GMT time* standard pressure 
level* area number* and blocketta number. 

* 216-225 

Procedure SCAN activates the scanning procedures. 

* 229-253 

Procedure WINDCOMP decomposes a vector wind into x and y 
components. 

» 257-262 

Procedure INTP contains a formula for linear interpolation. 
267-276 

This routine prints out the weather data which are stored 
in the array space. 

(One may delete this routine if not needed). 

281-313 

This routine converts the weather data in the standard 
pressure levels to the weather data in flight levels by using 
the interpolation formula. 

3J9-3;52 

(his routine prints out the converted weather data above. 
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(only X component of vector wind and temperature). 


One may delete this routine if not needed. 


« 336-340 


Initiai values; 


CO* conversion factor between degree and radian. 


MACH* 0. 82 for OCIO. 


P1.P2.P3* flight levels. 


Ql*<12.Q3* pressure altitudes for standard pressure levels. 


* 3S0-353 


The graph points are labelled by numbers running from 0 to 118 


(note that the gpraph point numbers shown in Fig. 37. P. 110 


runs from 1 to 116). 


The graph is such that the point sets consist of subsets of 


points arranged along 'Meridians whose indices run from 


0 to i5 


Array A( index; contains the graph point number assigned to the 


most northerly point on each meridian. 


Array A is used for the computation of zone index (See the 


proiedure ZONED. 


3S8-3/2 


Array V has values assigned for geographic coordinates for 


graph points and check points on continents: 4-digit latitude 


4-digit longitude in degrees and hundreths. 


(See Fig. 37. P. 110). 


376-3 wl 


This routine fills Boolean array BK and BKQ with 'TPUE' 


382-3VS 


This routine denotes blocking in continental ai’'ways. 




ao 
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(connectivity in the graph). 

Booloon array antriat of BK and BKO indicata which graph point! 
ara connactad or hlochad <TRUE«bloekad« FAU3S«eonnactad). 

Array BK it uaad for tha lona inda> < 3 and array BKO la utad 
for tha xona indax > 8. 

BK(a* b< c >»FaU 8E (naans point b on mariaian a is connatad with 
point c on (naridian a4>i. 

Hara point b danotas tha b-th point from tha bottom on tha 
maridian. Exampla: BK(8«2*2>*FAI.8E maant graph point numbar 
3 on maridian 2 is connactad with graph point 9 on maridian 
3 as you saa Fig. 37. 

* 409-414 

Procadure 20NEI datarminas tha indax of tha zona associatad 
with a graph point number. • ■*— 

Parameter O^graph point number. 

For example* a graph point number 49 gives Z0NEK49 )k6. 

* 416-424 

Procedure OT dissects latitude and longitude from tha 
compressed coordinates. LALAT and LAlONO denote latitude and 
longitude respectively. Parameter K is an indax of array V. 

* 426-4*39 

Procedure LiS detarinines actual latitude and longitude values 
for a graph point. ( See procedure OT above). 

The Sign of longitude value is changed in the east of Creenwhich. 
Parameter U-graph point number. 

« 461-4S0 

Procedure CTQ prepares a time instant array 00* zonal ly 
specified* to be used for tha time interpolation later on. 
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A«0. 4 d«nottfl «n tttioiAttd flight tiA# of 24 minuttt p«r 
coneinpntdl zont. 

Tlif instant valuts run froA indaz 0 to indoz 14. 

*484-498 

Proe«duro HH eomputoo grid point voluot according to SELECT 
valuo 

Xd ^LECT^i « HH goto K-componant value of vector wind 

If 8 ElECT* 2 * Md goto g-coAponant valuo of vector wind. 

If 8ELSCTO » HH goto toAperature value. 

These values are computed under consideration of time factor 
for dynamic process. 

Parameter V and W denote arrag index numbers for row and column 
respectively. Parameter P«index for flight levels (31000 ft«l* 
35000 ft-3 and 39000 ft»3). • <•— 

« 503-523 

Procedure (3E0P computes x and g component values of vector wind 
and temperature value for an arbitrary point bg using procedure 
HH above. Bi-linear interrpolation scheme is used in this 
computation Parameter C»flight level index. 

* 538-538 

Procedure (^EOOlSf computes a great circle distance between two 
arbitrary points Pl(LATi«LONi> and P2(LAT8« L0H2> . 

* 543-5b0 

Procedure OEOtlGRlD converts coordinates of two arbitrary 
points on sphere to their coordinates on the geometric grid. 

PP and 06 denote two grid point numbers. 

* 555-638 

Procedure PART26E0H computes length of segmti.t by summation of 
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contribution from tuboogioontt CLSNCTHd)) «nd o truo courso 
onglo of 0 flight fogmont. 

Cl ond E2 donoto unit voctort for tbo trwo courco. 

Normal 1 ltd coord Inatoo (running from 0 to I) art storad in 
arrag AB. 

Tnofto valuas oarvo also as MSigbtng factors labor on. 

Array LENOTH(I) contains langtfcs of swbsagmants on tha sghara. 
*A40~aS6 


Th 


s um mation formulas. 


« 6^-684 

Procsdura DRIFT computas a drift angla by using i«y componants 
of vactor wind. 

•»# 487-741 • *“ 

Procsdura rtETPROC procassos matsorological data bsforo 
optimization process labor on. 

Paramstsr 8«flight laval indax. 

* 495-705 

ihis rouCina computas avaraga waathar paramatars for aacb 
sub«aijmsnt. 

* 704-708 

This routine computas avaraga weather parameters for each 
segment in order to get tingle drift angle of the segment. 

■> 7JJ-714 

This routine computes flight time in each subsegment. 

- 717-740 

This routine computes temperature deviation* flight time* 
air distance* flight neading* x*y components of vector wind* 
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OP POOR quality 

«nd Xi V cumpontnttt of tri#« oir tpood in taeli fttgiMnt «nd 
fliitht. itvti. 

TUtof eomputtd v«lu*c «r« «ll «v«r«g« v«lut« along a oagaitnt 
ondor conoidaration uting praviootlg oalculattd eontridutiona. 
TEH^tofliporaturo. 

T0EV(8>«t»mparatura daviatlon. 

TlME(8}*f light tima in flight laval indax 8. 

AtRDX8T(8)*air diatanca. 

LE«di«tanca of flight aag«ant. 

WIX(8)"x coAponant of vaetor wind. 

TAX(8)»trut air tpaad componant. 

♦ 749-785 

Proeadura TABLE and CLXraCORT ara tab la lookup procaduraa 
for picking up parformanca data for<«d>CiO. 
a 787-830 

Procadura RSADO raadt in parformanca data and procadura 
MtXTEO writas out thosa inputad data. 

Ona mag dalata procadura URITEQ and linat 813-830 if 

not n<»a<icii. 

* 833-e:>2 

Thi» routine r«au< in input paramatars for a flight 
planning. In caat of prascrioad track flight (ROU'E >•* 

1111 ca&a)« tna routine reads in 0. 

In another ca»a« it reads in SfCctarting point nu^riber) 
and STKand point number) and decides 0 value. 

(!■ # of graph points on the prescribed route or ( a of 
xones between source and destination) - 1. 


* 8A1-8A4 
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Procedure SO end JQ are conversion formulas for indices 
used in subsegment algorithms bg using array A( index) and 
array F( index) respectively. 

« 8d»6»889 

Procedure OESCLIriB computes performance by using the table 
lookup procedure in climb (if MQ^O) or descent(if MQ not« 0). 

The calculation is somewhat complicated by the fact that the 
location of TOC is not fixed. The position eihere the flight 
levels off may vary due to the influence of the winds in the 
climb phase. 

For mean climb vector wind we may take 60 or 73 % of the vector 
wind and the climb time is usually reduced by a few minutes to 
account for take-off and acceleration and thus we may take 2 
minutes for it. 

Array DECLTIKE contains climb time or descent time. 

Array DEClDIST contains climb distance or descent distance. 
Array OECLFUEl contains climb fuel or descent fuel. 

Array AIROIST contains climb air distance or descent air 
distance. 

It IS noted that in the descent xone. the performance table 
is nob used and thus specific expressions are chosen according 
to company's planning policy. 

( See formula 6.31 in P. 120-121). 

892-933 

Procedure LIME prints out a flight plan. 


v3ti-l242 


Procedure SPACEOPT is a keg procedure for optimalization. 
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tn thit prdc«durt* crucial !• that soma paraoatars may nut 
«urpa>»f uppar bounds* dor asamplti aircraft waight. 

This raqutroscht inclusion of tavoral protactivo statoAonts 
in order to ensure proper functioning of the operational scheme. 
FBQ determines whether the flight plan computation will be 
performed baekward<FBQa^AL8E> or forward (FBQ«TRUC>. 

FB determines whether cost(FB*l)« fuel(FB«0> or flight time 
(FBe-i) will be optimised. 

BI determines whether the navigation regime is free in the 
horisontal<BI>BFALSE) or bounded bg one position(6I*TRUE>. 

BO determines whether the cruising alitude is free(B0>^AL8E) 
or bounded (BO bTRUE) 

RR istake-off weight or landing weight. 


♦ 993-961 


Procedure PREP computes distance and flight time* and prints 
out a line of flight plan. 

Al«segment distance 
Bl«'^egmrnt flight time. 

AAsiRuE denotes a climb phase and AA^FALSr a descent phase. 
963-1010 


Procedure EDITING prints out proper heading lines for a flight 
plan. 


10l;2-t0rt4 

Procedure (iQ computes segment contribution* fuel* etc. 

If 0»-l* computation is for in-flig*ht direction anc if 0*^1* 
then computation is for opposite flight direction. 

- tOia-1038 

Procedure CL computes the aircraft weight (GRUG) ano total fwel* 
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flight tim« or costft at and point of a flown •ogment. 

(the contribution of fuel in btepe included). 

ParaAeter Ml denotes segment fuel(kg) and M8 segment time(min). 

In the climb phase no step contribution is taken into account as 
required (BCLIMB»TRUE). 

GO denotes flight level index at the graph point when the computation 
starts. 6 denotes the flight level index to be investigated. 

Step consumption: UP«700 kg/4000ft 

DOUN—28 kg/4000ft. 

(surplus of segment fuel) 

ITERATION becomes TRUE when (1) back-tracking (2) first zone is 
reached. 

(3) TOW is unknown and Wb is known (see P. 120). 

In array ROW the aircraft- weight is- etored. 

CRWQ denotes aircraft weight at the end point of segment. 

QUANT denotes cost* fuel or time at the end point of flown segment. 
«1045-1052 

computation in climb or descent by procedure CL. 

If ED is TRUE then it prints out the results. 
rl055-l057 

saf ey uard ing statements. 

computes segment contribution of changes in aircraft weight. 

^ute that fuel consumption is expressed as loss of aircraft weight. 
The table for specific range <P 116) is used. 

Formula 6.27 (P.119) requires an estimate of weight halfway a flight 
Segment. For that purpose coefficients as shown in table 2* P. 116 

■319 required <rl* r2) 
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jumping back to labol JMl is for safeguarding < that th# actually 
derived aircraft weight does not surpass the maxifflally allowed 
weight. 

♦1071-107S ' 

computes the contribution for the last segment. 

If ED is TRUE then it prints out the results. 
eiOSl 

DRWQ and QUANT attain absurd values for safeguarding operational 
performance 
*1087 

EO becomes TRUE when all degrees of freedom in the horizontal and 
vertical are lost* for example when a call of procedure SPACEOPT takes 
place in the final mode for flight pten computation along the (optimal) 
track found in previous calls. 

•when navigation regimes is bounded by one point and cruising altitude 
is bounded). 

*1088-1113 

Initialization and preparation of all parameters needed for the 
algorithmic process of optimization. 

rBQsTKUE means flight plan computation will be performed forwards* 
while FBQ^FALSE backwards 

Array F is filled in line 1244-1248 (see procedure FF). 

The process can be activated for arbitray begin- and end-points. 

In order to reduce storage space a renumbering is made for all subset 
points. So array F is defined and functions completely analogous to 
array A for wuole graph point set. 

Array STOREl and STOREIO contain the assigned numbers of labelled 
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en<l-point& on meridians in the subgraph. ^ WOR QUMJTY 

Zone cycle runs from line 1119 to 1220. 

The tone cycle repeats as many times as the number of tones. 

A cycle for graph points along meridian runs from line 1122 to 1216. 
This cycle repeats as many times as the number of graph points in 
the current meridian. 

A cycle for graph points along next meridian runs from line 1127 to 
1214. 

This cycle repeats as many times as the number of graph points in 
the next meridian. 

A cycle for flight levels runs from line 1158 to 1212. 

This cycle repeats as many times as the number of flight levels used. 

» 1126 

conversion for indices- used in subsequent algorithms. 

» 1131-1151 

If the navigation is free (BI*FALS£> in the horizontals it checks 
whether a segment (II* JJ) is blocked or not. 

BK or BKQ s TRUE means airway is blocked. 

If blocked* it skips computation for the segment. 

* 1152-1156 

conversion for indices (point G4). 

determines geographic elements for a segment between points G3 and 64. 
See procedures GEOMGRlD and PART2QE0M. 

* I 162 - 1 1 VI 

In case of backtracking in a flight plan computation an iteration 
^rocec>s is put into action in the climb zone 

This IS done in order to determine the (unknown) take-off weight by 
iteration in such a way that the climb parameters do match with the 
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(laraffleter <• found during backtracking whan arriving in tha climb zone. 
The iteration etarts with a take-off weight 20000 kg below the 
(known) maximum take~off weight. 

Array ROW contain* the <optimal) value of aircraft weight found 
during optimalization> using a 2 one*~cycle. 

The iteration finishes when either the weight difference at the point 
of matching is less than 10 kg or the number of iteration steps 
exceeds 10 (poor convergence). 

In each step take'-off weight is adjusted# see line 1181-1191. 

* 1195-1203 
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similar computations but not backtrack ing case or not take-off 
weight adjustment case. 

# 1204-1210 


These statements are crucial optimalization criteria. 

« 1225-1241 

Preparation of elements which are required for a subsequent call 
for procedure SPACEOPT with lessened number of degrees of freedom 
including a final call for computation and presentation of flight 
plan data. In this final stage no degrees of freedom are left. 
Procedure SPACEOPT then merely operates along a predefined track 
I the solution of previous calls). 

1244-12-fS 


Procedure »*F determines a zone index array F for a subset of grapo 
points analogous to array A. (array F is a suo— zone array). 


1250-1351 

A call for procedure PROCES results in the production of a flight 
plan alony whatever route is desired. 


CETAluS 


♦ia79-12f9 


OTFOORQIMIIIV 

8p«Cific«tion o# limits in herisontm o# flrnpli points sleng 
m«ritfi«ns. 

Array 8T0RSI and array 8T0REIQ start loitftr limits and upptr limits 
rsspactivaly for tht numtar of grapt points on aacli maridian. 

Xf routs >■ liil* than production of flight plan along a prasribad 
routs. 

In this Casa tha procadura 8PACEOPT uorks in a daganarata mods. 

* i29i>l296 

Specification of limits in tha vertical. 

Array 8T0RE0 and array STOREOQ store loufcr limits and upper limits 
respactivalv for the number of flight levels on each meridian. 

« 1298-1303 

Procedure EP indicates- that the fbeght planning computation blows 

up. 

* 1305-1316 

See P. 123-124 for procedures 88< S88 and 8888. 

These procedures compute landing weights. 

NN contains landing weight or take-off weight. 

1318-1324 

Procedure TW is for safeguarding against overloading. 

- 1327-1351 

activates process as follow: 

If FACTORIZAION is TRUE< then the optimalization takes place first 
:n the horizontal and is followed by an optimaliiaion in the 

vi»rti«*al . 

Ttie fiillowing steps occur: 

a) optimal! zati on in the horizontal based on time. 
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b) optimdiliZdtion in tii* vertical* ucing th« tracit solution found 
in (a) jnd bosed on fuol* ttmo or costs. 
ic) criffiputatiun of flight plan along solution found in (b>. 

Xf thtf optimal ization takes place in free space* the sequence of 
calls results in: (three procedure calls) 

(i> optimali zatzon in space based on time* fuel or costs. 

Ui> optimalization in the vertical through horizonal track found in 
(i)« based on time* fuel or costs. 

(lii) compilation flight plan along track solution found in (ii). 

(•lote: (ii) could be bgpassed as the solution is found already 

in <i). But (ii) can be generated with a slightly different landing 
u<eight. The compilation of a flight plan along a prescribed route 
passes through all three procedure calls. 

This means that in fact tha computaaten is repeated threefold* however 
with properly tuned landing or take-off weights. 

In order to protect aqainst subsequent calls of procedure PROCES* 

FOUTE is assigned by 1000. 

1353-13*4 

Read statement in case that a flight plan compilation is desired 
along a pre-icrined track. 

This track is specified by graph points indicated by their numbers. 

•= 136:j-i369 

Computation of an estimate of flight time to be used for the 
imat i on of other parameters. 

“"he same* wir-h distance 

Procedure CiQ makes time mstanc array OD for each zone (mer i diar. : . 
-rray E is for storma the number of graph points used on a meridian. 


•' 1371 -l3S9 
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par#orm« « proctsa for ti)« pf^tfvction of « flight plan. 

Actual poromotoro whon calling procoduro « 

u*a landing lueight (BB*T1UIC>. 

cost optiflialiiation (FFB8«i>. 

FFfi0«O moans fuel opt* FFBB«*1 moans time opt. 

regularity percentage <PR»0. 03). 

fligh levels used botmoon 8i«i for 31000 and 82a3 for 39000. 

S3sl> the value is immaterial if FACTORIZATION is FALSE. 

» 

For required flight plans t 

1. ttihen a flight plan for the optimal cost track is needed* 
call ^proces ( 1. 3. 1* 0. 03* i* TRUE* FALSE) ^ 

2. when a flight plan for the optimal fuel track is needed* 
call 'proces ( 1. 3. 1* 0. 03* 0. TRUE* FALSE) ^ 

3. eihen a flight plan for- the least, heeie track is needed* 
call proces ( 1* 3* 1. 0. 03. -1* TRUE* FALSE) ^ 

For the flight plan 1* 2 or 3 above* FACTOR ZZATI0N«FAL8E. 

For other flight planning simulations* FACTORXZATION«TRUE. 



He*THeR data in 3 standard pressure levels 
3 AX 16 DATA OUT OF 3 AX 20 GRIOPOINTS • 
TENPERATUREt DIRECTION* SPEED • 
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ORIGINAL PACi: If5 
OF POOR QUALITY 



) 

) 






la. 


DETAILED FLIGHT PLAN ♦** 


MIN TRACK IN SPACE 

fltime 


CH 1190581 0.00 


NO. 

HEAD 

FL 

TOC 

46 

310 

113 

46 

310 

108 

66 

310 

103 

63 

310 

97 

62 

310 

91 

61 

310 

87 

51 

350 

77 

54 

310 

65 

62 

310 

51 

66 

310 

37 

72 

350 

21 

90 

350 

11 

89 

V'SlO 

6 

101 

'» 310 

/ 

2 

118 

350 

TOD 

0 

121 



TAS 

WIND 

DIST 

ACCD 

TIME 

ACCT 

BURN 

weight 








231393 


29 



20 

20 

5756 

225637 

479 

39 


162 

5 

25 

708 

224929 

470 

54 

232 

394 

26 

51 

3741 

221188 

477 

65 

150 

544 

17 

107 

2335 

218853 

477 

71 

142 

686 

15 

123 

2167 

216686 

476 

81 

99 

785 

11 

134 

1472 

215214 

475 

lOO 

180 

965 

19 

152 

2689 

212525 

480 

112 

216 

1181 

22 

214 

2914 

209611 

484 

116 

407 

1588 

41 

255 

5525 

204086 

486 

106 

364 

1952 

37 

332 

4898 

199180 

478 

104 

336 

2208 

35 

406 

4445 

194743 

470 

93 

336 

2624 

35 

441 

4297 

190446 

402 

55 

185 

2809 

21 

502 

2521 

187925 

480 

44 

93 

2902 

11 

513 

1330 

186595 

477 

41 

147 

3049 

17 

530 

2105 

184490 

477 

15 

51 

3100 

6 

536 

683 

183807 


12 


3192 

16 

552 

1203 

102604 


trip FUEU 40789 
C0ST11129.51 ' ' 


END OF RUN ♦♦♦ 


DETAILED FLIGHT PLAN 



MiN TRACK IN SPACE 

fltime 


CH 1190561 0*00 


NO* 

HEAD 

FL 

TMP 

TAS 

WIND 

DIST* 

ACCD 

time 

ACCT 

BURN 

weight 












251697 

TOC 

262 

310 



3 


136 

23 

23 

6611 

245086 

1 

262 

310 

-4 

477 

3 

0 

136 

(^10^ 

23 

3 

245088 

4 

262 

310 

-1 

480 

-22 

242 

. A76 

32 

55 

5065 

240023 

6 

279 

310 

2 

484 

-36 

115 

493 

16 

111 

2475 

237548 

17 

292 

310 

3 

465 

-73 

57 

550 

6 

120 

1313 

236235 

30 

266 

310 

6 

487 

-39 

397 

947 

54 

214 

8260 

227975 

44 

268 

310 

7 

489 

-55 

394 

1341 

55 

308 

8026 

219949 

60 

270 

310 

7 

466 

-64 

394 

1735 

56 

404 

7912 

212037 

75 

265 

310 

2 

464 

-75 

394 

2129 

59 

503 

7972 

204065 

67 

261 

310 

-2 

460 

-78 

179 

2306 

27 

530 

3568 

200477 

92 

257 

310 

-4 

477 

-74 

161 

2469 

24 

554 

3119 

197338 

96 

236 

310 

-5 

476 

-80 

66 

2535 

10 

604 

1283 

196075 

104 

239 

310 

-4 

477 

-68 

217 

2752 

32 

636 

4054 

192021 

109 

246 

310 

-3 

478 

-50 

189 

2941 

27 

702 

3330 

166691 

113 

219 

310 

-2 

479 

-26 

160 

3121 

24 

726 

2974 

185717 

TOO 

237 

310 

-1 

480 

-22 

77 

3198 

10 

736 

1262 

184455 

115 





-17 


3283 

16 

752 

1203 

163252 


TRIP FUEL 66445 
C0ST153b3.19 


m Mkt 




F''-'^'. \4r^ \ . 


♦♦♦ detailed flight plan ♦♦* 


MiN TRACK IN SPACE 

fltime 


CH 1190581 0.00 


NO. HEAD FL TMP TAS 

WIND 

DIST 

accd 

TIME 

ACCT 

BURN 

weight 

^2*J2-60r-KH — 0-rOOOO-r-K)0 







251699 

TOC 262 310. 

2-r328^+01* '0.0000 »>00 

3 


136 

23 

23 

6611 

245088 

1 262 310 -4 477 

2,3280»+01 1.9204»+04 

3 

0 

136 

0 

23 

0 

245088 

4 282 310 -1 480 

2.3280»+01 9.4894»+03 

-22 

242 

378 

32 

55 

5065 

240023 

8 279 310 2 484 

2.3280»+01 5.0680»+03 

-36 

115 

493 

16 

111 

2475 

237548 

l7 292 310 3 465 

2.3280»+01 3.2515»+04 

-73 

57 

550 

8 

120 

1313 

236235 

30 268 310 6 487 

2.3280»”*-0l 3.2736»+04 

-39 

397 

947 

54 

214 

8260 

227975 

44 260 310 7 489 

2,3280*'»-01 3.3479» + 04 

-55 

394 

1341 

55 

308 

8026 

219949 

60 270 310 7 488 

2.3280r+01 3,5169»+04 

-64 

394 

1735 

56 

404 

7912 

212037 

75 265 310 2 484 

2.3280»+01 1.632?»+04 

-75 

394 

2129 

59 

503 

7972 

204065 

87 261 310 -2 480 

2,3280»+01 1.4450»-^04 

-78 

179 

2308 

27 

530 

3588 

200477 

92 257 310 -4 477 

2.3280»-*-01 6,0011»403 

-74 

161 

2469 

24 

554 

3119 

197358 

98 238 310 -5 476 

2.3280»-^01 1.9l66»+04 

-80 

66 

2535 

10 

604 

1283 

196075 

104 239 310 -4 477 

-68 

217 

2752 

32 

636 

4054 

192021 


2.3280»+01 1.5906»+04 


J 


109 246 310 

2.3280»'^01 

-3 478 

1.4323»404 

-50 

189 

2941 

27 

702 

3330 

188691 

113 219 310 

1.5610»-^01 

-2 479 

6.1052»'»-03 

-26 

180 

3121 

24 

726 

2974 

105717 

TOD 237 310 

1.501O»+O1 

-1 480 

6,1052»+03 

-22 

77 

3198 

10 

737 

1262 

184455 

115 


-17 


3283 

16 

752 

1203 

183252 


TRIP FUEL 68447 
C0ST153B3.64 


*♦* END OF RUN ♦♦♦ 


Program completed in 
end of run: 137 cru 

BASE cru; 137.035 
CORE factor; 1.000 


130.1426 CRU. 

account; 510il55»l 
SEO* NO. 122425 
cate; 22 MAY 01 
limit; 7500 


cards: 4 in 0 out 

CPU time: 130.79 SEC. 

time: 19:59:28 


TOTAL 


5852 CRU 


DETAILED FLIGHT PLAN ♦♦♦ 


MlN TRACK IN SPACE 

Fuel 


CH 1190581 0.00 


NO. 

head 

FL 

TMP 

TAS 

WIND 

DIST 

ACCD 

TIvE 

ACCT 

BURN 

WEIGHT 












232421 

TOC 

46 

310 



29 

- 

121 

20 

20 

5796 

226625 

113 

46 

310 

-2 

479 

39 

41 

162 

5 

25 

694 

225931 

108 

66 

350 

2 

475 

55 

232 

394 

26 

51 

4622 

221309 

103 

63 

350 

1 

474 

66 

150 

544 

17 

108 

2407 

218902 

97 

62 

350 

1 

474 

74 

142 

686 

16 

123 

2210 

216692 

91 

62 

350 

1 

474 

85 

99 

785 

11 

134 

1496 

215196 

87 

51 


2 


102 

180 

965 

19 

153 

2598 

212598 

77 

55 

350 

4 

477 

114 

216 

1181 

22 

215 

3002 

209596 

65 

62 

350 

6 

479 

116 

407 

1568 

41 

256 

5493 

204103 

51 

66 

350 

6 

479 

iio 

364 

1952 

37 

333 

4814 

199289 

37 

72 


5 

^4^ 

106 

336 

2266 

35 

407 

4368 

194921 

21 

90 

.350 

4 


99 

336 

2624 

35 

443 

4317 

190604 

11 

89 

350 

5 

478 

65 

185 

2809 

21 

503 

2491 

188113 

6 

101 

390 

3 

476 

47 

93 

2902 

11 

514 

1989 

186124 

2 

118 

• 

390 

4 

476 

41 

147 

3049 

17 

531 

2023 

184101 

TOO 

121 

390 

4 

477 

22 

26 

3075 

3 

534 

361 

183740 

0 





17 


3192 

19 

553 

1271 

182469 


trip fuel 49952 
COST11320.37 


DETAILED FLIGHj PLAN **♦ 




MIN TRACK IN SPACE 

fltime 

CH 1190581 0.00 


NO. 

head 

FL 

TMP 

TAS 

WIND 

DIST 

ACCD 

TIME 

ACCT 

BURN 

WEIGHT 












23362U 

TOC 

46 

310 



29 


122 

20 

20 

5843 

227781 

113 

46 

310 

-2 

479 

39 

4j0 

162 

5 

25 

678 

227103 

108 

66 

310 

-4 

478 

56 

232 

394 

26 

51 

3787 

223316 

103 

63 

310 

-4 

477 

66 

150 

544 

17 

108 

2360 

220956 

97 

62 

310 

-5 

476 

73 

142 

686 

15 

123 

2185 

218771 

91 

61 

310 

-5 

477 

84 

99 

^785 

11 

134 

1488 

217283 

87 

51 


2 

^47^ 

102 

180 

965 

19 

152 

3331 

213952 

77 

54 

310 

-1 

480 

112 

216 

1181 

22 

214 

2986 

210966 

65 

62 

310 

2 

483 

116 

407 

1588 

41 

255 

5554 

205412 

51 

66 

310 

4 

486 

109 

364 

1952 

37 

332 

4930 

2004S2 

37 

72 

f 350' 

5 

^4^ 

106 

336 

2288 

35 

406 

5096 

195386 

21 

90 


4 

'.'478^ 

99 

336 

2624 

35 

442 

4331 

191055 

11 

89 

310 

1 

482 

62 

185 

2809 

21 

502 

2580 

188475 

6 

101 

310 

-1 

481 

49 

93 

2902 

11 

513 

1333 

187142 

2 

118 

350 

4 

477 

42 

147 

3049 

17 

530 

2717 

184425 

TOD 

121 

310 

-4 

477 

19 

51 

3100 

6 

536 

727 

183698 

0 





14 


3192 

16 

552 

1203 

182495 


trip fuel 51129 
COST11480.87 


*** END OF RUN *♦* 


CH 1190581 0.00 


NO. head fl tmp tas wind dist accd time acct burn weight 


232450 


TOC 

46 

310 



29 


121 

20 

20 

5797 

226653 

113 

46 

310 

-2 

479 

39 



5 

25 

694 

225959 

106 

66 

350 

2 

475 

55 

232 

394 

26 

51 

4622 

221337 

103 

63 

350 

1 

474 

66 

150 

544 

17 

108 

2407 

218930 

97 

62 

350 

1 

474 

74 

142 

--686 

16 

123 

2210 

216720 

91 

62 

350 

1 

474 

85 

99 

785 

■9 

134 

1496 

215224 

87 

51 

350 

2 

475 

102 

180 

965 

19 

153 

2598 

212626 

77 

55 

350 

4 

477 

114 

216 

1181 

22 

215 

3002 

209624 fl 

65 

62 

350 

6 

479 

116 

407 

1588 

■9 

256 

5493 

204131 1 

51 

66 

350 

6 

479 

110 

364 

1952 

37 

333 

4814 

199317 

37 

72 

350 

5 

478 

106 

336 

2288 

35 

407 

4361 

194956 

21 

90 

350 

4 

478 

99 

336 

2624 

35 

443 

4317 

190639 

11 

89 

350 

5 

478 

65 

185 

2809 

21 

503 

2491 

188148 1 

6 

101 

390 

3 

476 

47 

93 

2902 

11 

514 

1991 

186157 

2 

116 

390 

4 

476 

41 

147 

3049 

17 

531 

2023 

184134 

TOD 

121 

390 

4 

477 

22 

26 

3075 

3 

534 

361 

183773 

0 





17 


3192 

19 

553 

1271 

182502 

TRIP 

r" r ^ iL T 1 

FUEL 

'X T n 

49948 




















, 18 . 


*♦♦ detailed flight plan *•♦ 


N track in space 
FLTIME 


NO. head FL TMP TAS WIND 


TOC 46 310 


‘♦6 310 


3 63 310 


97 62 310 


91 61 310 


87 51 3 50 


54 310 


65 62 310 


66 310 


72 350 


90 350 


89 310 


6 101 310 


2 118 350 


TOD 121 310 


-2 479 


4 478 


-4 477 


-5 476 


-5 477 


2 475 


460 


2 483 


4 486 


5 478 


-4 477 


TRIP FUEL 51118 


• 22 


End of run 


CH 1190581 


DIST ACCD TIVE ACCT BURN WEIGHT 


233643 



I 


216 1181 


407 1588 


36 


336 2288 


20 5844 227799 


25 678 227121 


51 3787 223334 


108 2360 220974 


123 2185 218789 


134 1488 217301 


152 3331 213970 


214 2986 210984 


255 5545 205439 


332 4930 200509 


406 5096 195413 


51 3100 


3192 


536 727 183728 


552 1203 182525 


478 

99 

336 

2624 

35 

442 

4331 

191082 1 

482 

62 

185 

2809 

21 

502 

2580 

188502 

481 

49 

93 

2902 

11 

513 

1333 

187169 

^ 477 

42 

147 

3049 

17 

530 

2714 

184455 


li 


f 


> 



































*♦* DETAILED FLIGHT PLAN ♦** 


MIN track in space 


I 


NO. head tmp tas wind 


TOC 262 -^10 3 


1 262 310 -4 477 


8 279 310 


17 292 310 


30 268 310 


44 268 350 


60 270 350 


75 264 350 


87 260 350 


2190561 0.00 


dist accd time acct burn weight 


252998 


136 24 24 6678 2U6320 


24 0 246320 


242 378 32 56 5105 241215 


6 111 2490 238725 


1320 237405 


214 6306 229099 


310 9187 219912 


394 1735 57 407 8142 211770 


394 2129 60 507 7979 203791 


179 2308 28 534 3549 200242 


-36 

H 

93 ] 

16 


55 


9 


3 


l'> 

92 

257 

350 

1 

474 

-74 

161 

2469 

24 

559 

3063 

197179 1 

1 « 

96 

238 

350 

1 

474 

-78 

66 

2535 

10 

609 

1250 

195929 L 

1 

104 

239 

350 

1 

474 

-66 

217 

2752 

32 

641 

3926 

192003 1 

h 

109 

(VI 

350 

2 

475 

-46 

189 

2941 

27 

707 

3197 

188806 j 

I iW 

113 

219 

390 

2 

474 

-28 

180 

3121 

24 

731 

3590 

185216 1 


TOD 

236 

^0 

0 

473 

-27 

55 

3176 

7 

739 

875 

184341 ! 

1 > 
1 

115 





-20 


3283 

19 

758 

1271 

183070 |l 


TRIP FUEL 69928 
C OS T 1 5 b'66.9S 

n?67nu 


detailed FLI'^Ht plan *** 




























- 

1 



1 

**! 

I f ' 

t 1 

1 

i 20. 

**♦ detailed flight 

_PL^N 

1 

- 1 

L 





1 “IN track in space 

_^UEL 



i 

f 

J 



CH 21Q0581 0.00 ! 

K ^ 

1 • 

■ — - - 

NO# head Fl TvP 

TA«; WIND DIST 

ACCD 

time acct Burn weight 1 



TOC 262: 31 n 

__ l__2_62 ...31 n 

_ 4 282 _31n 

_.8 ^27.9 .310. 

17 292 310 


30 268 _31_n. 

44 268 35n. 

60 270 35n 








^ {^7 

:l...‘L8p 
J2-- 9C4 
3 485 


_8_JL87. 

7 480 

7 480 





2490 

238990 

1.2JL. 

JL3.2.3. 

_?.5JZ.66.7. 


397 9 47 
394 1341 
394 Tt35 


















215 8322 229345 

311 9203 220142 

409 8142 212000 












































/VO^XAV^ 






^ ELT DC10»l*8l052e»065609 

COhMENT KNMI/KLM 020169 SP6 NAV 20 JN6» 
NAvIGATIOU-FLIGHTPlANNInG-MODULE for the 
FLIGHTPLANS OVER ThE NORTH ATLANTIC. 


CDC Io). 


PRODUCTION OF 


THE MAIN features 

are: 

1. INCLUSION OF 6 METEOROLOGICAL PARAMETER FIELDS* 
temperature AND WIfiD DATA FOR 30U»250 AND 200MB 

valid for two standard times 12 HOURS apart. 

2. weather info, comes from NATIOfjAL WEATHER SERVICE* 

3. INCLUSION OF A DYNAMIC PROCESS BASED ON INTERPOLATION OF 


TWO PARAMETER FIELDS. 

4. OPTIMALISATION oF FUEL* COSTS OR FLIGHT TIME IN SPACE AS 
WELL AS IN the HORIZONTAL AS IN THE VERTICAL* 

5. IN - and outbound TRAFFIC (10= 1 OUTBOUND 10= -1 INBOUND) 

6. FLIGhTPLAN production for optimum tracks* SPECIFIC ROUTES 
ETC* E*G* ALTERNATIVE AtC MINIMUM ROUTES* 

7. blocking E*6. sector BLOCKING OVER THE NORTH ATLANTIC 
blocking In airways and ATC RESTRICTED AREAS* 

6. STANDARD CRUISE AND USE OF PERFORMANCE TABLES* 


AKW>WKCDaCa(;a 64 

• •♦•I 

♦ ♦•*2 

*«*W4 

*•♦*5 

*•♦*7 
♦ •♦•p 
**«*9 
**♦*10 
****11 
****12 
****13 
****14 

****15 

****16 

****17 

♦***18 

****19 

****20 


CONSTANT KQ= 14 » W= 3S 

integer P1»P2»P3»01»Q2»Q3»S»TTT» 

GltG2»G3*G4.Bl*B2fB3.MQ»G10»G2Q*Q» 

RES1»RES2»RES3»I»J»N1»N2»N3»M»G» 

M0.Ml»M2rM3»M4»M5.M6»K. lOr 
A1»A2»A3»A4»A5»T1»T2»T3»K0UNT» 

Kl.K2»Ll»L2»KlQ.K20»L:.0»L20»KKl.KK2»LL4-fLL2» 

FLUR.TAXI » TOW. GRW» reserve »MAXTOW»MAXLW» 

DATE. ST .STl .ROUTE. distance. FLTIME. BURN. 

grwq.towo.lw.cost.le.tel* 

REAL C0.CG.P.C1.X1.X2.Y1.Y2.X.Y.E1.E2.E3. 

GGl . GG2 » D . ANGLE . L AL At . LALONG . TT . LAI . LO l . L A2 . L02 . 
XX1»XX2,YY1>YY2.AID1.AID2.AID3.AID4.XSTER.YSTER. 
MINTIME.TC.TASM.WX.WY'TAS.ENDU.TEMS 

INTEGER array A(-1..kO+1). V(0..41>» 

TEMPTIME.TEMPDIST.TEmPFUELC 1. .99. 1..9). MW(1..3. 1..7). 
RANG£(1..3. 1..33)S 
ARRAY DD(0..KO)» 

ARRAY WCX.WCY.TrCd. .120. 1..34)S 

BOOLEAN ED» , , 

BOOLEAN array 13K(0. .2' I. .6. 1« .6) . BKQ ( 9. .KQ. 1 . .6. 1 . .6) » 
FORMAT FOKE2. SlO. d7.0» A5). 

PEJT (E2. • '. Al). 

FOTtXS. *16. A2). 

F0D(*I6. Al). 

F02(X5. *17. a2)S 
LOCAL LAUF.L ERR.EOPS 


COMMENT PROCEDURE dljATA* SCAtiS WEATHER DATA AND OBTAINS 
The REOUIREL weather info* FOR north AMERICA AND 

north ATLANTIC AREA % 


****21 

****22 

***♦23 

***4i24 

♦***2b 

***•26 

**••27 

•***28 

•***29 

***•30 

****31 

****32 

****33 

***•34 

****3b 

•***36 

****37 

***•38 

***•39 

***•40 

****41 

****42 

***•43 

****44 

****4b 

***•46 

•*•*47 

****46 

***•49 

***•50 

•**•51 

***•52 

•*•*53 



PROCEOUFiE KDATAS 

BEGIN IrjTEGER 1 *X»Y 1 .Y2»Y3#R»T»U»VS 
STRING ARRAY wD (7. . 1 , . 120 » 1 • .34 ) S 
STRING S(fll)S 

STRING TW(AG(2)» SP(3>» TP(2))» 

string BA(BD(1> »BL(3) »AR(2) »TM(D »B(3) > * 

GW ( UG ( 2 ) * W 1 ( 7 ) » WWM ( W 2 ( 7 ) » W 3 ( 7 ) » W 4 ( 7 ) ) » HH ( 2 ) » G ( 4 ) ) . 

SM <WT(2) »Ur(2) »WS(3) )S 

local label FINI» 

FORMAT F(A»SeO)4 
FORMAT F1(*S8» Al) » 

F2( ’BLOCkETTE «=»»S3»* AREA = NORTH AMERICA* » Al)» 

F4( ’BLOCkETTE «=*»S3»* AREA = NORTH ATLANTIC* » Al)» 

F5(3(X3»S7)» Al>» 

FW(*I6»A1) » 

PJ(E2»* *»A1)S 


PROCEDURE CHOPS 
BEGIN INTEGER K»LS 
IF I LEO 45 THEN 

BEGIN 6W(1»36)=S(I»36)S I=I+36S 

IF I GTR 80 THEN 

BEGIN ReAD(PCF(*NV»S*)»F»S»FINI)$ I = IS ENDS 

END 

ELSE BEGIfi K=81-I5 Ga ( 1 »K)=S( I »K ) S REAo(PCF ( *NWS* ) *F»S*FINI )S 
L=36-KS Ga(K-*1» L)=S(1*L)S I=L*1S 

ENOS 

IF S(D EOL *\* then 

begin READ(PCF( *fJAS* ) »F»S*FInI )S I = 1S 
ENDS 

ENDS comment chops 


Procedure stors 
BEGIN Integer ji 

local label L3S 

STRING ARRAY SA ( 36 . . 1 . . 8 ) S 

FOR J=1 STEP 1 until 8 DO 

BEGIN CHOPS SA( 1 » 38 * •U)=GHS ENDS 
FOR J=1 STEP 1 UNTIL 4 DO 
BEGIN Gifc “SA ( 1 » 36 » • w) ) S 

IF alphabet IC then 

BEGIN wf<lTE( *W-ERROR* )S 
EfJDS 

WDU»7..Y1.x)=^.2S 
kk0tl»7. .Y2»X>=*^3S 
WD<1»7. .Y3*x>=1i«4S 
Y1=Y1-1S Y2=Y2-1S Y3=Y3-1S 

ENDS 

X=X-1S Y1=Y1^4S Y2=Y2-*4S Y3=Y3+4$ 
FOR J=5 STEP 1 until 8 DO 
begin 6o=SA(1»36..J)S 


write(wwia»s Goto l3S 


•**•54 

*••♦55 

*♦••57 

*•*•58 

*•*♦59 

•***60 

*•••61 

*•••62 

*••♦63 

****64 

*•••65 

****66 

•♦•♦67 

****66 

*•••69 

••••70 

♦*♦♦71 

♦•••72 

•♦♦♦73 

♦♦♦♦74 

♦*♦•75 

•>♦♦•76 

♦♦♦♦77 

♦**•78 

♦♦♦♦79 

♦♦♦♦00 

♦*♦•81 

♦*♦♦82 

♦ *♦♦33 

♦♦♦*84 

♦♦*♦05 

♦♦♦•06 

♦ ♦♦•87 
♦♦♦♦88 
*♦♦♦89 
**♦•90 
*♦*•91 

*•**92 

♦*•♦93 

***•94 

*••• 9 ^ 

•*••96 

♦**•97 

••••98 

•**•99 

*« *'♦100 

*♦••101 

***• 10 ? 

*•♦*103 

♦•♦■.>104 

♦•♦*105 

♦♦♦•lot 

♦•♦•107 


IF ALPHABETIC (WWW) THEN 

GOTO L3* 

♦*♦♦106 

BEGIN WRITECW-ERRORMS WRITE(WWW)$ 


ENDS 


♦♦♦♦110 

WD (l»7«»Yl»x) “W2^ 


♦♦♦♦111 

WDd»7..Y2»x)=W3* 


♦♦♦♦112 

W0d»7..Y3»x)=W4S 


♦♦*♦113 

Y1=Y1-1* Y2=Y2-1S Y3=Y3-1* 


♦♦♦♦114 

END* 


♦♦♦♦115 

END* comment STOR* 


♦♦♦♦116 

♦♦♦♦117 

♦♦♦♦lie 

PROCEDURE AREA* 


♦♦♦♦119 

BEGIN Integer k»l* 


♦♦♦♦120 

IF I leg 71 THEN 


♦♦♦♦121 

BEGIN BAd»10)=S(l»10)* 1 = 1 + 10* 


****122 

IF I GTR 80 THEN 


♦♦♦♦123 

BEGIN READ(PCF( »NWS* ) »F»S»FINI)S 1=1* 

ENDS 

*♦♦♦124 

END 


♦♦♦♦125 

ELSE begin K=81-I* BAd»K)=S(I »K)* 


♦♦♦♦126 

READ(PCF(»NWS* ) »F»S»FIND* L=10-K* 


♦♦♦♦127 

BA(K+l*L)=Sd»L)* I=L+1* 


♦♦♦♦128 

end* 


♦♦♦♦129 

IF S(D EGL *\* THEN 


♦♦♦♦130 

BEGIN READ(PCF( *NwS’ ) »F»S*FInD* 1 = 1* 


♦♦♦♦131 

END* 


♦♦♦♦132 

END* comment AREAS 


♦♦♦♦133 

♦♦♦♦134 



♦♦♦♦135 

PROCEDURE FINDNEXT* 


♦♦♦♦136 

BEGIN 


♦♦♦♦137 

IF Sd) EQL • » THEN 


♦♦♦♦138 

FOR I = I WHILE Sd) NEO *♦’ DO 


♦♦♦♦139 

BEGIN READ(PCF(’NWS*) » F»S»FINI)S 1=1* 


♦♦♦*140 

END* 


♦♦♦♦141 

END* 


♦♦♦♦142 

♦♦♦♦143 

• 


♦♦♦♦144 

PROCEDURE GETDATAS 


♦ ♦♦♦145 

BEGIN INTEGER Gi..T»A»B»C»D»E»F» J»S»T* 


♦♦♦♦146 

LOCAL Label l5» nam» nat* 


♦♦♦♦147 

A=20S B=40* C=60* D=80* E=100* F=120S 


♦♦♦♦148 

AREAi 


♦♦♦♦149 

IF AK EOL *21 » THEN 


♦♦♦♦150 

BEGIN IF TM EQL *1* THEN 


♦♦♦♦151 

BEGIN GMT=0* X=l4* Y1=A* Y2=B* Y3=C* 


♦♦♦♦J52 

WPITECGMT = 0000*)* 


♦♦♦♦153 

end 


♦♦♦♦154 

else if TM EOL *3* THEN 


♦♦♦♦155 

BEGIN GMT=12S X=14* Yl= D* Y2=E* Y3=FS 


♦♦♦♦156 

WRITE (*GiV.T = 1200*)* 


♦♦*♦157 

END 


♦♦♦♦158 

ELSE BEGIN WRITE (* W-ERROR *) * GOTO L5* END* 

♦♦♦♦159 

GOTO NAM* 


♦♦♦♦160 

END 


♦♦♦♦161 


ELSE IF AR EQL »22* THEN 
begin if TM EQL »1* THEN 

BEGIN GMT=Oi X=3**» Y1=AS Y2=Bi Y3=C$ 
WRITE<*GM = 0000*)$ 

END 

ELSE IF TM EQL *3’ THEN 
BEGIN GMT=12$ X=34$ Y1=D$ Y2=E$ Y3=F$ 
WFITE<’6MT = 1200*)$ 

END 

ELSE BEGIN WRITE ( ’WERROR •) * GOTO L5$ ENDS 
GOTO NATS 
END* 

WRITE( *W-ERROR* )* 

GOTO L5S 

NAM.. FOR S=1 STEP 1 UNTIL 7 DO 

begin FOR T=1 STEP 1 UNTIL 5 DO 
BEGIN 

WRITE(F2»BD* 

STORS X=X+1* areas 
IF NUMERIC (BL) THEN 

IF integer tBL) GEQ 259 THEN 
FOR J=1 STEP 1 until 3 DO CHOP 
ELSE J=1 

ELSE BEGIN WRITE (* W-ERROR *) S GOTO L5S 
ENDS 

ENDS 

X=X-2S 

IF GMT EQL 0 THEN ^ 

BEGIN Y1=AS Y2=BS Y3=CS END 
ELSE BEGIN y1=D* Y2=ES Y3=FS END* 

IF S LSS 7 THEN AREAS 

END* 

GOTO L5S 

NAT.. FOR S=1 STEP 1 UNTIL 10 DO 

BEGIN FOR T=1 STEP 1 UNTIL 5 DO 
BEGIN 

WRITE tF4»BD* 

STORS X=X+1* areas 
IF NUMERIC (BL) THEN 

IF integer (BL) GEQ 253 THEN 

FOR J=1 STEP 1 UNTIL 3 DO CHOP 
ELSE J=1 

ELSE BEGIN WRITE ( ’W-ERROR’ )* GOTO L5* 
END* 

ends 

X=X-2* 

IF GMT EQL 0 then 
BEGIN Y1=m? Y2=BS Y3=C* END 
ELSE BEGIN Yl=D$ Y2=E$ Y3=FS END* 

IF S LSS 10 THEN AREAS 

END* 

l5.. end* comment GETDATAS 


♦*♦♦162 
♦•♦♦163 #1 
♦♦♦♦164 y 
• ♦♦♦165 II 
••♦•166 d 
♦♦♦♦167 ^ 

♦♦♦♦168 
♦♦♦♦169 y 
♦♦♦♦170 « 

♦♦♦♦171 j 
♦♦♦♦172 
♦♦♦♦173 
♦♦♦♦174 ! 
♦♦♦♦175 
♦♦♦♦176 
♦♦♦♦177 
♦♦♦♦178 . 
♦♦♦♦179 
♦♦♦♦100 
♦♦♦♦181 
♦♦♦♦182 
♦♦♦♦183 
♦♦♦♦184 
♦♦♦♦185 
♦♦♦♦186 
♦♦♦♦187 
♦♦♦♦188 
♦♦♦♦189 
♦♦♦♦190 

♦♦♦♦191 

♦♦♦♦192 

♦♦♦♦193 

♦♦♦♦194 

♦♦♦♦195 

♦♦♦♦196 
♦♦♦♦197 
♦♦♦♦196 H 

♦♦♦♦199 « 

♦♦♦♦200 E 

♦♦♦♦201 I 
♦♦♦♦202 I 

♦♦♦♦203 I 
♦♦♦♦204 I 
♦♦♦♦205 P 
♦♦♦*206 
♦♦♦♦207 
♦♦♦♦206 
♦♦♦♦209 
♦♦♦♦210 
♦♦♦♦211 
♦♦♦♦212 
♦♦♦♦213 f 
♦♦♦♦214 I 

♦♦♦*215 I 


procedure scans 

BEGIN 

local label BEGN»L0S 
BEGN. . READ IpCFt ’NWS* ) »F»S»FINI 
I = 1S 

LO.. if S(I»2) eql then 

BEGIN I=I+38S GLtDATAS FINDNEXTS GOTO LOS 

END 

ELSE GOTO BEGNS 
ENOS comment scans 

COMMENT computation 0^^ X»Y COMPONENTS OF WIND VECTORS 


PROCEDURE WlNDCOMP(X»Y)S 
VALUE XrYS 
rriTEGER X»YS 
buGiN real as 

IF X GTR 36 TK;N WRITER ’WIND ANGLE ERROR 
IF X LEO 9 THEN 

BEGIN A=(90-X*l0)*CGS 
WX=“Y*COS<A)S 
WY=-Y*SlNtA)S 

END 

ELSE IF X LEO 10 THEN 

BEGIN A=lX*10-90)*CGS 
WX=-Y*COS< A)S 
WY=Y*SIN(A)S 

END ^ 

else if X LEO 27 THEN 

BEGIN A=(270-X*10)*CGS 
WX=Y*COS(A)S 
WY=Y*SIN(A)S 

end 

ELSE BEGIN A= (x*10-270 ) *CGS 
WX=Y*COSU)$ 

WY=-Y*SIN(A)S 

ENDS 

ENDS comment WiNDCOMpS 


)S 


comment linear INTERPOLATIONS 

INTEGER PROCEDURE INTPtF»Fl »F2»FlQ»F2Q)S 
integer F»F1»F2»F1G»F2QS 
i INTP= ItF-Fl ) *F2Q •* (f 2-F)*F1Q) /{F2“Fl)S 

I comment intps 


♦♦♦*216 
*♦♦*217 
*♦**216 
*♦♦*219 
♦**♦220 
*♦**221 
*♦**222 
*♦**223 
****224 
*♦**225 
♦♦*♦226 
*♦♦*227 
****228 
*♦**229 
*♦**230 
*♦**231 
*♦*♦232 
♦♦♦♦233 
♦♦♦♦234 
♦♦♦♦235 
♦♦♦♦236 
*♦**237 
*♦**238 
*♦**239 
****240 
*♦*♦241 
***♦242 
*♦♦*243 
*♦**244 
****245 
♦♦♦♦246 
♦ ♦♦♦247 
♦♦**248 
♦***249 
*♦♦*250 
*♦**251 
****252 
****253 
*♦**254 
****255 
*♦*♦256 
***♦257 
***♦258 
*♦*♦259 
*♦♦*260 
*♦♦♦261 
*♦♦*262 
*♦*♦263 


SCANS 
FiNl* • 

FOR R=0 STEP 1 UNTIL 5 DO 
BEGIN T=R*20S V.RlTE<PU?S 

WRITEI’* weather data 


*♦*♦264 

*♦♦*265 

*♦*♦266 

*♦**267 

/ ♦♦♦*268 

IN 3 standard pressure LEVELS *’)* ****269 


I 


ENDS 

COMMENT 


WR1TE(»* 34X16 DATA OUT OF 34X20 
WRITE I •* temperature* DIRECTION* 
WRiTEt* *)4 WRITEt* ’)* 

FOR U=1 STEP 1 UNTIL 34 DO 

WRITE( Fl* 20* FOR V=T+20 STEP -I 
UNTIL T+5 DO WD ( 1 »7. . V»U> ) * 


CONVERSION OF WEATHER InFO. IN STANDARD 
BY interpolation SCHEME TO WEAjHER INFO. 


gridpoints ♦•)$ 

SPEED ♦’JS 


STEP 1 UNTIL 20 » 61 STEP 
J=1 STEP 1 until 34 DO 


1 UNTIL 80 DO 


BEGIN M1=M1-50S T1=”T1S END* 


FOR 1=1 
FOR 
BEGIN 

SW=WDd»7..I»J)* 

mi=integer(dr)s 
M 2=INTEGER(WS)S 
T1=INTEGER(WT)5 
IF Ml GEQ 50 THEN 
WlNDC0MPtMl»M2)* 

M1=WX* M2=WY* 

SW=WD(1»7..I+20*J)* 

M3=INTE6ER(DR)5 

M4=INTEGER(WS)5 

T2=lNTE6ERlWT)S 

IF M3 GEQ 50 THEN BEGIN M3=M3-50S T2="T2* END* 

W1NDC0MP(M3»M4)* 

M3=WX* M4=WY* — - 

SW=WD(1*7. .1+40* J)* 

M5= I NT EGER (DR)* f.',6=INTEGER ( WS ) * 

T3=INTEGER(WT)* 

IF MS GEQ 50 THEtj BEGIN M5=M5“50S 
WINDC0MP(M5*M6)* 

M5=WX* M6=WY* 

WCX( I* J)=INTP(P1 *0l*Q2*Ml»M3)* 

WCY(I* j)=INTP(Pl»0l'Q2*M2*M4)* 

TPC( I» j)=INTP(P1*q1*Q2»T1*T2)* 

WCX( 1+20* J)=INTP(P2»Q2»Q3»M3»M5)* 

WCY( 1+20* J)=lNTP(p2* 02*03* M4»M6)* 

TPC ( 1+20* J)=INTP(P2 *02*03 *T2*T3)* 
aCX( 1+40* J)=INTP (p3* 02*03 *M3*M5)* 

WCY( 1+40* J)=INTP(p3* 02*03 *M4»M6)* 

TPC ( 1+40* J)=lNTP(p3 *02*03 *T2*T3)* 

END* 


T3=-T3* END* 


***♦270 
****271 
****272 
****273 
**•♦274 
****275 
****276 
****277 

PRESSURE LEVELS ****278 
IN flight LEVELSS****279 

****260 
****281 
****282 
****283 
****284 
****285 
****286 
****287 
***♦288 
****209 
****290 
****291 
****292 
****293 
****294 
♦***295 
****296 
****297 
♦***296 
***♦299 
****300 
****301 
*♦♦*302 
***♦303 
****304 
****305 
****306 
**♦*307 
***♦308 
****309 
****310 
****311 
♦***312 
****313 


FOR A1=0 STEP 1 UNTIL 5 DO 

BEGIN A2= Al*20* WRITE(PEJT)* 

WRITE( ’*INTERP0LATED GRIDPOINT VALUES IN FLIGHT LEVELS*’)* 
WRITE ( ’*X-COMPONENTS OF WIND DATA*’)* 

WRITE(’ ’)* WRITE(’ ’)* 

FOR A4=l STEP 1 UNTIL 34 DO 
WRITE (FW» 20. FOR A5=A2 + 20 STEP -1 
until A2+1 DO WCX(A5* A4))t 


***♦315 

****316 

****317 

***♦318 

♦***319 

****320 

♦***321 

♦**♦322 

♦*♦♦323 


END* 


♦♦*♦324 

END* CO^^MENT EtjD OF WEATHER DATA COLLECTION* ♦♦♦♦SaS 

♦♦♦♦32b 

Pl=31000* P2=35000* P3=390001 ♦♦♦♦327 

01=30065* 02=33999* 03=38662* ♦*♦♦328 

TEL=0* ♦♦♦♦329 

C6=0. 0174532925* ♦♦♦♦330 

TASM=0.82*36.9826* ♦♦♦♦331 

♦♦♦♦332 

WRITE(»* SCANNING WEATHER DATA IN NORTH AMERICA AnD-»)* ♦♦♦♦333 

WRITE(* NORTH ATLANTIC AREA ♦♦)* ♦♦♦♦334 

♦♦♦♦335 

♦♦♦♦336 

WDATA* ♦♦♦♦337 

♦♦♦♦338 

♦♦♦♦339 

♦♦♦♦340 


COMMENT RELATIONSHIP BETWEEN GRAPH POINT NUMBERS AND ZONE INDICES*****34l 
A(-l)= -1* A(0)= 0* aU)= 2* A(2)= 6 * A(3>= 11* A(4)= 26* ♦♦♦*342 
A(5>= 41* A( 6 )= 56* A<7)= 71* A( 8 >= 86 * A(9)= 90* A(10)= 96* ♦♦♦♦343 
A(ll)= 102* A(12)= lo7* A(13)= 111* A(l4)= 114* A(15)= 115* ♦♦♦♦344 

♦♦♦♦345 
♦♦♦♦346 
♦♦♦♦347 
♦♦♦♦348 
♦♦♦♦349 
♦♦♦♦350 
♦*♦♦351 
♦*♦♦352 
♦♦♦»353 
♦*♦*354 
♦♦♦♦355 
♦♦♦♦356 
♦♦♦*357 
♦♦♦♦358 
♦♦♦♦359 
♦♦♦♦360 
♦♦♦♦361 
♦♦♦♦362 
♦♦♦♦363 
♦♦♦♦364 


COMMENT determination OF CONNECTIVITY AMONG GRAPH POINTS. ♦♦♦♦SbS 

FALSE=COfJNECTED» TRUE=BL0CKED * ♦♦♦♦3b6 

FOR 1= 1 STEP 1 until 6 DO ♦♦♦♦367 

FOR J= 1 STEP 1 UNTIL 6 DO ♦*♦♦368 

BEGIN ♦♦♦♦369 

FOR MQ= 0cl»2 DO EK(VQ»I»J) r TRUE* ♦*♦♦370 

FOR MiQ= 9 STEP 1 UNTIL KQ DO BKQ(MQ.I.J) = TRUE ♦♦♦*371 

END* ♦♦♦♦372 

BK(l»l»i)= BK(2»1»1)= BK(1.2»2)= 0K(2»2»2)= ♦♦♦♦373 

BKGt9»l»l)= BKG(10.1.1)= BKQ(11»1»D= PKQ(12»1»D= ♦♦♦♦374 

BKQtl3»l»D= BKG<14»1»1)= BKQ(9»2»2)= BKQU0»2»2)= ♦*♦♦375 

bKGlll»2»2)= BKG<12»2»2)= BKQ(13»2»2)= DKG(10»3.3)= ♦♦*♦376 

BK3lll*3»3)= BKGU3»3»3)= BKQ(11»4»4)= BK<2»3»4)= 8 k(2»4»5)= ♦♦♦♦377 


COMMENT COORDINATES OP CHECK POINTS* 

V(0)= 52290477* 

V(l)= 51850115* V(2) = 53640150* 

V(3)= 51000200* V(4) = 51990539* V(5) =--53580300* 
V( 6 )= 55000200* V(7) = 49980632* V( 8 ) = 51840849* 
V(9)= 52700892* V(10)= 53500630* VUD= 55490459* 
V(12)= 54271005* V(13)= 48905454* V(14)= 51375560* 
V(15)= 53735697* V(l 6 >= 58476263* V<17)= 47005800* 
V(18)= 48545856* V(l9)= 55506000* V<20)= 57006000* 
V(21)= 53286035* V(22>= 54836683* V(23)= 46156006* 
V(24)= 48006000* V(25)= 49846439* V<26)= 50506500* 
V(27)= 51006750* V(2a)= 52006900* V(29)= 45006300* 
V(30)= 46176459* V(3l)= 47506600* V(32)= 485768265 
V(33)= 50007050* V(34>= 43836608* V(35)= 44846867* 
V(36)= 45327179* V(37>= 45467385* V(38)= 41287003* 
V(39)= 42367099* V(40)= 43007250* V(4l)= 40647378* 



Erjo 

ELSE IF U LEO 4l THEN 
bEGIN 

LALAT= 20+Ul LALONG= 20 
End 

ELSE IF U LEG 56 THEN 
BEGIN 

LALAT= 5+Ui LAL0N6= 30 

END 

ELSE IF U LEO 7l THEN 
BEGIN 

LALAT= U-llS LALONG= 40 

END 

ELSE IF U LEO THEN 
BEGIN 

LALAT= U-26S LALONG= 50 

END 

ELSE DT(U-74) 

ENDS C0K.MENT LIS$ 

procedure CTqS 

comment this subroutine prepares the time instants at WHICH 

IN each zone THE N.ETEOROLOGICAL PARAMATER WILL BE 
DERIVED from both PARAMETER FIELDS» COMPOSITE CHARTS 
IN TIME are SIMULATEDS 
BEGIN 

INTEGER JO^MOS 
real A»SUM»RESi 
A= 0.4S 

JO = IF 10 EOL 1 THEN KO ELSE 05 SUM= 05 

FOR 1= KQ-JO step 10 UNTIL JO DO 

BEGIN 

SUM= SUM+(IF I GTR 3 and I LSS 8 THEN 2*A ELSE A)S 
DD(D= SUM 

MO = ZONEHST)$ RES= IF 10 EGL 1 THEN CD(MO> ELSE 
UD(M0-1)5 

FOR 1= 0 STEP 1 UNTIL KQ QO PD(I)= IF TT 6TR 12 
THEN TT-12 ELSE TT+A/E+DD ( I ) -RES5 
ENDS COMMENT CTGi 

COMMENT computation qF GRID POINT VALUES FOR DYNAMIC PR0CESS5 

real PROCEDURE HH (V * W * F » SELECT ) 5 
VALUE F »V»V rSELECTS 
INTEGER F»V»w»SELECTS 
BEGIN 

M5= (F-l)*201 
Mfa= (F+2)*20S 
IF SELECT EOL 1 THEN 

BEGIN M1=WCX ( V+M5* W) 5 M2=WCX < V+M6* W ) 5 

HH=n 12“DD<M) ) *M1 ♦ D0(M)*M2) /125 END 
ELSE IF SELECT EOL 2 THEti 

BEGIN M3=WCY (V+M5»'W)S ;.‘4=WCY ( V'^’MG » W ) 5 


*♦♦*432 

♦**♦433 

***«43U 

♦•**435 

♦**•436 

♦***437 

♦*•*438 

****439 

****440 

****441 

****442 

****443 

****444 

****445 

****446 

****447 

****443 

****449 

♦**♦450 

♦**♦451 

♦♦♦♦452 

♦*♦♦453 

♦•♦♦454 

****465 

♦**♦456 

♦***457 

♦♦*♦458 

♦**♦459 

*•**460 

♦**♦461 

♦•**462 

***♦463 

*♦**464 

♦♦*♦465 

****466 

♦*♦*467 

**♦*466 

♦♦♦♦469 

♦*♦*470 

♦•♦*471 

****472 

♦**•473 

*•*♦474 

*••*475 

•♦♦*478 

♦•**477 

***•478 

♦♦♦*470 

♦**♦480 

♦♦*♦461 

♦♦♦♦482 

♦♦♦*483 

♦*♦♦484 

♦♦•♦485 



^ — - 



^ HH=( (12“DD(M) )*M3 ♦ 0DCM)*M4) /12S END 

♦*♦♦486 


ELSE BEGIN Tl=TPCtV+M5» W)* T2=TPC ( V+M6» W)* 

♦♦♦♦487 

1 

Hh=( (12-D0(M) )*T1 ♦ DD(M)*T2) /12$ END* 

♦♦♦♦488 i 

; 

END* COMMENT HH* 

♦♦♦♦489 

*«*«49Q 

1 

comment computation oF grid point value by BI-LINEAR INTEPP0LATI0N****4ql J 

1 

FOR AN arbitrary POINT (TEMP.» DIRECTION AND SPEED)* 

♦***492 

' 1 


♦♦♦♦493 


REAL PROCEDURE GEOP < V» W»C »SELECT )S 

♦***494 

1 

VALUE C»SELECT»V»a* 

♦♦♦♦495 ' 

! 

INTEGER C»SELECT* REAL V»W* 

****496 


BEGIN 

♦***497 

1 

Integer vi»wi* real a»b»bi»b2»b3»B4S 

♦***498 


' Vl= ENTIER(V)S Wl= ENTIER(W)* 

****499 1 


B1=IF SELECT lQL 1 THEN HH ( VI » wl »C » 1 ) 

♦♦♦♦500 


ELSE IF SELECT EQL 2 THEN HhWI »Wl»C»2) 

♦♦♦♦501 


ELSE Hh(Vl»Wl»C»3)S 

♦♦♦♦502 


B2=IF SELECT EqL 1 THEN Hh ( V l+1 » Wl »C r 1 ) 

♦♦♦♦503 


ELSE IF SELECT EOL 2 THEN HH< Vl+1 » Wl »C»2) 

♦♦♦♦504 


1 ELSE HH(V1+1»W1»C»3)* 

♦♦♦♦505 


B3=IF select EoL 1 THEN Hh ( V l+1 » Wl+1 »C » 1 ) 

♦♦♦♦506 


ELSE IF SELECT EQL 2 ThEN HH ( Vl + 1 » Wl + 1 » C * 2 ) 

♦♦♦♦507 


. ELSE HH(V1+1»W1+1»C»3)5 

♦♦♦♦506 


‘ B4=if select EqL 1 then Hh(V1»w1+1»C»1) 

♦♦♦*509 


ELSE IF SELECT EQL 2 ThEN HH ( VI » Wl+1 »C » 2) 

♦♦♦♦510 


ELSE HH(V1»W1+1»C»3)S 

♦*♦♦511 


A= Vl+l-V* B= wl+l-W* 

♦*♦♦512 


GEOP= A+B+B1+ ( 1-A ) +B+B2+ 1 1-A ) ♦ ( 1-6 ) +B3+A* ( 1-B ) ♦B4S 

♦*♦♦513 


END* COMMENT GEOP* 

♦♦♦♦514 



♦*♦♦515 


COMMENT COMPUTATION oF GREAT CIRCLE DISTANCE BETWEEN TWO 

♦*♦*516 


POINTS BY using THE GONIOMETRIC RELATION* 

♦♦♦*517 

1 


♦*♦♦516 

1 

real Procedure geodist(lati»loni»lat2»lon2)* 

♦*♦*519 


VALUE LAT1»LAT2»L0N1»L0N2* 

♦♦♦*520 jl 


REAL LATI »LAT2»L0N1 »L0N2* 

♦*♦*521 1 

1 

BEGliJ 

♦*♦♦522 j 


local label skp* 

**♦*523 ’ 


if LATI EOL LAT2 AND LONl EOL LON2 THEN 

***♦524 

1 

BEGIN GEOUIST=0* GOTO SKP* END* 

**♦*525 


6EODIST=60*ARCCOS(SIN(LAT1+CG)+SIN(LAT2+CG) 

♦*♦*526 } 


+C0S<LAT1+Cg)*C0S(LAT2*CG)*C0S( (L0N2-L0N1 )*CG) )/CG* 

♦***527 


Skp. . 

♦***526 


END* COMMENT GEODlST* 

**♦*529 

♦*♦•530 


comment computation of coordinates of two geographical points 

****531 


ON the grid plane* 

♦***532 

♦***533 


procedure GE0MGR1D<PP»G6)* integer pp»gg* 

****534 


comment X1»X2=R0v^ indices ArjD Y1»Y2=C0L. indices * 

♦***535 


1 begin 

****536 

1 

LIS(PP)* LA1= lALAT* L01= LALONG* 

♦*♦•537 

1 

LIS(GG)* LA2= lALAT* L02= LALONG* 

♦♦**538 


Xl=(70-LAl )/2.5* Yl=(13C-L01)/5* 

1 

♦***539 




X2=(70-LA2)/2»b* Y2= ( 130-L02 ) /5S 

E.NOS> COMMENT GLUmGRID* 


COMMENT COMPUTATION qF GREAT CIRCLE DISTANCE AND TRUE COURSE 
OF A SEGMENT FLOWNS 
procedure PART2GE0M* 

BEGIN 
REAL AS 
INTEGER NS 

REAL T0TL»LA11,LA22»L011 »L022»XX1»XX2»YY1»YY2S 

Kl= ENTIERlXDS K2= Kl + l$ Ll= ENTIER(YI)* 

L2= Ll+l* K1Q= ENTIER(X2)$ K20= KlQ+l» 

LIQ= ENTIeR‘Y 2)5 L20= LlO^l* 

Nl= ABS(KIC-Kl)* 

N2= ABS(LIQ-Ll)* 

N3 = IF KIQ EOL K1 AND LlO EQL LI THEN 0 
ELSE IF N1 LEG N2 THEN N2 ELSE NIS 
BEGIN 

INTEGER I»JS 
LOCAL label LABIS 
ABtO)= OS AB<N3-»-1)= 1* 

IF N3 EQL 0 THEN GOTO LABIS If N3 EQL Nl THEN 
BEGIN 

IF Kl GTR KIQ THEN 


BEGIN 

RES1= KIS 

RES2= K2QS RES3= -1 

END 

ELSE 

BEGIN 

RES1= K2S RES2= KlQS 
RES3= 1 
ENDS 

FOR 1= ReSI step RES3 UNTIL RES2 DO 
BEGIN 

J= RLS3*< I"RFS1)+1S AB<J)= (I-Xl)/ 
(X2-X1) 

END 


END 

ELSE 

BEGIN 

IF Ll GTR LIQ THEN 
BEGIN 

KES1= LIS 

RES2= L2QS RES3= -1 

END 

ELSE 

BEGIN 

RES1= L2S RES2= LlQS 
RES3= 1 
EtJDS 

FOR 1= RE’-il step RES3 UNTIL RES2 DO 
BEGIN 


****S4l i 

♦*♦♦542 

♦♦*♦543 

♦♦♦♦544 

♦♦♦*545 

♦♦♦♦546 

*♦♦♦547 

♦♦♦♦548 

♦*♦*549 

♦♦♦♦5b0 

****5bl 

****552 

****5b3 

****5b4 

****555 

*♦**556 

****557 

****558 

****559 ; 

****560 

***♦561 

****562 

****563 1 

****564 I 

****565 

****566 

****567 

***♦568 

****569 

****570 

****571 

****572 

****573 

♦♦♦♦574 

****575 ' 

♦♦♦♦576 

****577 

*♦**578 

****579 

****530 

****581 

****582 

***♦533 

****58ti 

***♦535 

****586 

****587 

***^588 

***♦539 

****590 

****591 

****592 

*♦*♦593 


J= RES3*U-RES1) + 1* AB(J)= d-Yl)/ 

(Y2-Y1) 

END 

ENDS 

TOTL=OS 

LA11=LA1S l011=L01S 
N=N3+1S 

FOR 1=1 STEP 1 UNTIL N DO 
BEGIN 

XX2=AB(I)*X2 + (l-ABd) )*X1S 
YY2=AB<I)*Y2 ♦ (l-ABd)»*YlS 
LA22=70-(XX2*2.5)S 
L022=130-(YY2*5.0)S 

LENGTH ( I ) =6E0D I ST ( LA 1 1 » LOl 1 » LA22 » L022 ) S 
LA11=LA22S L011=L022S 
TOTL=LENgTH ( I ) +TOTLS 
ENDS 
D=T0TLS 

TC=ARCC0S( (sIN(LA 2*C6) - SIN(LA1*C6) ♦C0S(D/60*C6) ) / 
(Slrj(D/ 60 *CG)*COS(LAl*CG) ) )S 
TC=TC/CGS 

IF SINUl02-L01)*CG) GTR 0 THEN TC=360 - TCS 
IF TC LEO 90 THEN 
BEGIN A=(90-TC)*CGS 

El=COS(A)S E2=SIN(A)S END 
ELSE IF TC lEO 180 THEN 

BEGIN A=(tC-90)*C6S <*• 

El=COS(A)S E2=-SIN(A)S END 
ELSE IF TC leg 270 THEN 
BEGIN A=(270-TC)*CGS 

El=-COS(A)S E2=-SINCA)S END 
ELSE BEGIN A= < TC-270 ) ♦CGS 

El=-COS(A)S E2=SIN(A)$ ENDS 

End 

ENDS COMMENT PAKT2GE0MS 

REAL Procedure su^‘l(A» b» arr)s 
VALUE A» BS 

INTEGER A» BS ARRAY ARRS 
BEGIN integer IS REAL SS 
S= OS 

FOR 1= A STEP 1 UNTIL B DO 
S= S + ARR(I)S 
SUMl = SS 

ENDS comment SUMS 

REAL PROCEDURE SUw2(A» B» ARRl» ARR2)S 
VALUE A» BS 

INTEGER A» BS ARRAY ARR1» ARR2S 
BEGIN integer IS REAL SS 
S= OS 

FOR 1= A STEP 1 UNTIL B DO 
S = S ARRl (I )*APR2d )S 


***«594 

*•*♦596 

***♦597 

****596 

****599 

**♦♦600 

*♦♦♦601 

***♦602 

*♦♦♦603 

*♦*♦604 

*♦♦♦605 

*♦*♦606 

**♦♦607 

*♦♦♦606 

*♦♦♦609 

****610 

*♦**611 

****612 

****613 

***♦614 

*♦♦*615 

****616 

****617 

****618 

****619 

♦*♦*620 

*♦*♦621 

*♦*♦622 

*♦*♦623 

*♦♦♦624 

****625 

****626 

*♦*♦627 

*♦*♦626 

*♦*♦629 

****630 

***♦631 

*♦**632 

*♦*♦633 

♦♦♦♦634 

*♦♦♦635 

*♦*♦636 

*♦*♦637 

*♦♦♦636 

♦♦*♦639 

*♦♦♦640 

*♦♦♦641 

*♦*♦642 

*♦**643 

♦♦♦♦644 

*♦♦♦645 

♦♦♦♦646 

♦♦♦♦647 


SUM2 = SS 

tNDS comment SU(.2i 

comment computation qF a drift angle for the segments 

real procedure DRiFT(X»Y)$ 

VALUE X»Y5 
REAL X»Yf 
BEGIN 

REAL SPD»ANG»DEG*Vk* 

LOCAL LABEL BR» 

IF X EQL 0 AND Y EQL 0 THEN 

BEGIN SPD=DEG=0» GOTO BR* ENDS 
SPD=SQRTU**2 ♦ Y**2)S 
w=abs<x)s 

IF SPD EQL 0 Then goto brs 
ANG=ARCCOS lW/SpD) /CGS 

IF X GEO 0 AND Y LEG 0 THEN ANG=270‘*'AN6 
ELSE IF X GEO 0 AND Y GEO 0 THEN 
ANG=270-AnG 

ELSE IF X LEO 0 AND Y GEO 0 THEN ANG=90-^AN'G 
else if X LEO 0 and Y LEO 0 THEN ANG=90-ANG 
ELSE WRITER *WInD DIRECTION ERROR* >S 
DE6=ANG-TCS 

BR • t 

DRIFT=ANG=SPD*SIN(DEG*CG) /TASS 
C0SA=C0S(ARCSIN<ANG) )S 

ENDS COMMENT DRIFT* 



procedure METPR0C(S)S VALUE SS INTEGER SS 
COMMENT PROCESSIfJG METEOROLOGICAL PARAMETERSS 
BEGIN 

INTEGER I* NS 

real hd*gs»tass* 

ARRAY TEMP»WNDX»WNDY »ENDUR(1» ‘NS+DS 
XX1=X1S YYl=Yi* 

N=N3+IS 

FOR 1=1 STEP 1 UNTIL N DO 
BEGIN 

XX2= AQ(I)*x2+tl-AB(I) )*X1S 
YY2= AB(I)*Y2+(1-AB(I) )*Y1S 
XSTER= tXXl+XX2)/2S 
YSTER= tYYl+YY2)/2S 
TEMP(I)= -geopuster»yster»s»3)s 
wNDX(I)=GEOf USTER»YSTER»S»1)S 
WNDY ( I ) =GEOP ( XSTER » YSTER * S » 2 > S 
XXl=XXi;S YYl=YY2S 

Er;Ds 

wx=suMi ( 1 »n»v;nox) /ns 

WY=SUM1 (1»N»WNdY) /ns 
TEM=SUM1 ( 1»N»TeMP) /NS 
TAS=TASM*S0RT (273*16+TEM)S 
AnGLE=DRIFT(WX»WY)S 


•♦♦♦646 1 

♦♦♦♦649 

♦♦♦♦6b0 

•♦♦♦6b2 
•♦♦♦6b3 
•♦♦♦6b4 
♦♦♦♦655 
•♦♦♦656 
♦♦♦♦657 
♦♦♦♦65B 
♦♦♦♦659 
♦♦♦♦660 
♦♦♦♦661 * 
♦♦♦♦662 
♦♦♦♦663 
♦♦♦♦664 
♦♦♦♦665 
♦♦♦♦666 
♦♦♦♦667 
♦♦♦♦668 
♦♦♦♦669 
♦♦♦♦670 
♦♦♦♦671 
♦♦♦♦672 1 

*♦♦♦673 ' 

♦♦♦♦67K 

♦♦♦♦675 

♦♦♦♦676 

♦♦♦♦678 

♦♦♦♦679 

♦♦♦♦680 

♦♦♦♦681 

♦♦♦♦682 

♦♦♦♦683 

♦♦♦♦684 

♦•♦♦685 

♦♦♦*686 

♦♦♦♦687 

♦•♦♦686 

♦*♦♦689 

*•♦*690 

♦♦♦♦691 

♦♦♦♦692 

♦♦♦•693 

♦•**594 

♦♦♦♦695 

*♦♦♦696 

♦♦♦♦697 

*♦♦♦698 

**♦♦699 

***♦700 

*♦*♦701 




AU20. . 


FOR 1=1 STEP 1 UNTIL N DO 

begin 

TASS=TASW*SqRT( 273.16 ♦ TEMP(I))S 
GS=TASS*COSA ♦ (WNDXd)*El ♦ WNDY(I)*E2)5 
ENDURU)=36 oOO*LENGTH(I) / GS4 
ENDS 

ENDU=SUP1 ( 1 » N » ENDUR ) S 
RES1= 31+(S-1)*MS 
W I X ( S ) =SUM2 ( 1 » N » ENDUR » WNDX ) /ENDUS 
W I Y ( S ) SSUM2 ( 1 * N » ENDUR » WND Y ) /ENDUS 
TEN=SUM2 ( 1 » N » ENDUR » TEMP ) /ENDUS 
TAS=TASM*S0RT(273.16+TEM)S 

TDEV(S)= IF RESl LEO 35.332 THEN TE« ♦! .98*RES1-15 ELSE 
SS^TEMS 

TIMEtS)= ENDUS 

AIRDIST<S)= ENQU *TAS/36000S LE= DS 

HEADIN6(S)=TC+ARCSIN(ANGLE)/CGS 

HD=HEADING<S)S 

IF hD LEO 90 then 

BEGIN HD=(90-HD)*C6S 

TAX(S)=TAS*COS(HD)S TAY (S)=TAS*SIN(HD)» END 
else if HD LEO 180 THEN 
BEGIN HD=<HD-90)*C6S 

TAX(S)=TAS*COS(HD)S TAY <S> =-TAS*SIN(HD> S END 
ELSE IF HD LEO 270 THEN 
BEGIN HD=t270“HD)*CGS 

tax(s)=-tas*cos(hd)s tay(S)=-tas*sin(hD)s end 
else begin HD= (HD-270) oCGS 

TAX(S)=-TAS*COS(HD)S 

tay(s)=tas*sin(hd)s ends 

ENDS comment METpROCS 

comment ’TABLE* AND ’CLIMBCORT* ARE TABLE LOOK-UP PROCEDURESS 

integer procedure table (P»R»MA)S 

VALUE P»RS integer P»R* INTEGER ARRAY MAS 

BEGIN 

Integer g»b»ds real ro»cs 

local LABEL LAB20S 

B=IF P EOL 0 then 4536 ELSE 5S 

C=IF P EOL 0 THEtJ GRWO ELSE TDEV(P.)S 

D=IF P EOL 0 THEN 127006 ELSE -15S 

R0=lC - D)/BS 

g=entier(ro)s 

R0=R0 - GS 

IF MA(R#G^1) EOL -1 THEN 

BEGIN TABLE=-iS GOTO LAB20S 
Ends 

TABLE=R0*MA lR»G+l ) ■*■ ( 1-RO) *MA(R»G>S 
ENDS comment tables 

Integer procedure climbcort(l»ma)$ 

VALUE LS INTEGER LS INTEGER ARRAY *^AS 


*••*702 

*•**703 

****704 

***•705 

****706 

•**»707 

**•*708 

*•**709 

***•710 

*•*<711 

****712 

***♦713 

•**♦714 

*•*♦715 

*♦•♦716 

*♦♦♦717 

*♦♦♦718 

*♦♦♦719 

♦♦♦♦720 

*♦♦*721 

****722 

****723 

♦**♦724 

*♦♦*725 

***^726 

*♦*♦727 

*♦**728 

*♦♦♦729 

♦♦♦♦730 

*♦♦♦731 

**♦♦732 

*♦*♦733 

♦***734 

*♦♦♦735 

****736 

*♦*♦737 

**♦•738 

*♦♦♦739 

*•♦♦740 

*♦*♦741 

*♦♦♦742 

♦♦♦♦743 

♦•♦♦744 

♦♦♦♦745 
♦*♦•746 
♦•♦•747 
*♦♦♦746 
*•♦♦749 
♦ ♦♦*75,0 
*♦♦♦751 




♦♦♦♦752 

♦•*♦753 

*♦♦♦754 

*•♦*755 



BEGIN 

INTEGER R*C»E1»E2»E3»E4S 
REAL G»R0»R1& 
local label LAB2H. 

6=(GRW0 - 127006)/4536» 

R=ENTIER<G)» 

R0=G-R» 

R=R*3 ♦ LS 
G=(TDEV(L) ♦ 25)/5» 

C=ENTIERlG)» 

R1=G-CS 

E1"MA(R»C)$ 

E2=MA(R+3»C)S 
E3=MA(R»C+1) J 
E4=MA(R*3»C+1)S 

IF El EOL -1 OR E2 EOL -1 OR E3 EQL -1 OR E4 EOL ■ 
BEGIN CLIMBC0RT=-1» GOTO LAB21S END* 
CLIMBCORT=RO*R1*E4 ♦ (1-R0)*a-R1)*E1 

♦ (1-R0)*R1*E3 ♦ R0*(1-R1)*E2S 

END* comment CLImBCORTS 

procedure re ado ( Y 0 » Y1 » Y2 » Y3 » ARR ) * 

VALUE Y0»Y1»Y2»Y3* INTEGER Y0»Yl»Y2»Y3* 

INTEGER array ARRS 
BEGIN 

FOR I=Y0 STEP 1 UNTIL Y1 DO 

READ(PCF( •DC10F» ) , FOR J=Y2 STCP 1 UNTIL Y3 
DO ARR(I»J) * ERR* ERR)* 

end* comment read* 

procedure WRITEC(yO*Y1*Y2»Y3»ARR)* 

VALUE Y0»Y1»Y2»Y3S 
integer Y0»Y1*Y2»y3* 
integer array ARRS 
BEGIN 

FOR I=Y0 STEP 1 UNTIL Y1 DO 

WRITE(FOT»Y3*FoR J=Y2 STEP 1 UNTIL Y3 DO 
ARR(I»J)» ERR)* 

end* comment WRITEG* 

READ0d»3*l *16 »PAuGE)* 

READQ(1»3» 17 »33» RANGE)* 

READQ ( 1 » 99 » 1 » 9* TEMPTIME ) * 
REA0Q(l»99»l*9*TEl'PDIST)* 

READO (1*99*1* 9* TEmPFUED* 

READQ(l*3*l*7*Mik)S 

write(pejt)* 

»»RITEt** PERFORM, At^CE DATA FOR DCIO ♦•)$ WRITE(’ 
KRIlEl’ ’)* WRITl(’ •)* 

WRITE(’* SPECIFIC RANGE ♦•)$ WRITEC ')S 
WRITEO (1*3*1* 16 *RArJGE)* 

WHITE0(1*3*17*33*rANGE)* 


•1 THEN 


♦•••7b6 

♦♦•*7b7 

♦••♦7b8 

♦♦•*7b9 

**»*750 

*••*762 

*•••763 

••*•764 

*•••765 

*••♦767 

•••*768 

*••♦769 

*•••770 

•♦•*771 

*••♦772 

•••*773 

•••♦774 

•♦•*775 

••••776 

♦♦•♦777 

•••♦778 

*♦♦♦779 

•♦*♦700 

•••♦781 

♦♦♦♦782 

♦♦♦♦763 

♦♦♦♦784 

•♦♦♦785 

♦♦♦♦786 

♦♦•♦787 

♦•♦♦786 

♦♦♦♦769 

♦•♦♦790 

♦♦♦♦791 

♦**♦702 

♦♦♦♦793 

♦**♦794 

♦♦♦*795 

♦*♦♦796 

♦*♦♦797 

♦*♦*798 

♦♦♦♦799 

♦*♦♦800 

♦•♦♦801 

♦♦♦♦802 

♦♦♦*803 

♦♦•♦004 

♦♦♦♦0C5 

♦♦♦♦P06 

♦♦♦♦007 

*♦♦♦806 

♦♦♦*g09 


j 


>»RITE<* *)» WRITE** *>* WRITEt* *)S WRITE*' *)S 
white*'* max. weight *')» WRITE*' ')S 
WRITE0*1»3»1»7»MW)» 

WRITE*PEJT)S 

WRITE*'* CLIMB time *')S WRITE*' ')S 
WRITEG * I » 99 » 1 » 9 » tempt IME ) S 
WRITE*PEJT)S 

write*'* climb distance *')S write*' ')» 

WRITE0*1»99»1»9»TEMP0IST)S 

WRITE*PEJT)S 

WRITE*'* CLIMB FUEL *'»» WRITE*' ')» 

WR I TEO * 1 » 99 • 1 » 9 » TEMPFUEU i 


READ*PCF* »OC10F' ) » FLUR» TAXI* GRW» RESERVE* TOw* 
maxtow»maxlw»date» io*ttt*route*done» err)* 

WRITE*' ')* WRITE** *)* WRITE*' ')S 
WRITE*' *** INPUT parameters ♦**')» 

WRITF.*F02» U* FLuP*TAXI»GRW»RESERVE*TOW» 

MAXTOW f MAXLw * DATE » 1 0 * TTT * ROUTE > * 

IF ROUTE 6EQ llll THEN 

begin 

READ*PCF* 'UClOF* ) » Q» DONE* ERR)$ 

WRITE*F02* 1* G)* 

Q=Q-1* 

END 

ELSE 

BEGIN 

READ*PCF* 'DClOF* ) » ST»ST1» DONE* ERR)* 

WRITE*F02* 2* ST* STD* 

G10= 20NEI*ST)S G20= 20NEI*STD* 0= ADS*G2Q-G1Q) 

END* 

TT - TTT* 

BEGIN 

INTEGER array sTOREI »ST0REI0*ST0REG* STOREGO* E ( 0 • • 0 ) * 
F(-1..Q)* 

FORMAT F04*X63» 16* A2)* 

F05*'TRlP FUEL'* 16* AD* 

F06*'C0ST'* De.2* AD* 

REAL array DEClTIME*DECLFUEL*DECLDIST*1..W)» 

Integer procedure so*hi*h2)* integer hi*h2* 

FC= A*610 + I0*hl-1 )’^H2* 

Integer procedure jo*hi *h2)* integer hi*h2* 

JQ= F*H1-1)+H2-ST0REI*H1)*1* 

PROCEDURE DESClIMB*S)S VALUE SS INTEGER S* 

COMMENT computation OF PERFORMANCE IN CLIMB OR DESCENTS 
BEGIN 

INTEGER alts 
LOCAL LABEL »iJMS 
real TASORTS 

TAS0RT=S0RT(TAX(S)**2 ♦ TAY*S)**2)S 
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ll 

i 

I 

I 

I 

I 

LAS30* • 

1 

) ,1 
1 I' 

t LAB31.. 





ALT= 31+(S-1)*4S 

DECLTIME(S)=IF mo EQL 0 THEN 60*CLIMBCORT^S»TEMPTIME) 
/lOOO else 3«1+0*41*ALTS 

DECLDIST(S)=IF MQ EQL 0 THEN CLIMBCORT (S»TEMPDIST) 
ELSE 9+3*ALT* ^ , 

DECLFUEL(S)=IF MQ EQL 0 THEN CLIMBC0RT(S»TEMPFUEL) 
ELSE 940+a.5*ALTS IF DECLTIME(S) LSS 0 THEN GOTO JUMS 

AID4=DECLDISTCS)*60/(DECLTIME(S)*TASQRT)S 
DECLDIST(S)= SQRT( (0»75*WIX(S)-»-TAX(S)*AID4)**2^ 

(0#75*wIY(S)+TAY(S)*AID4)**2)*(DECLTIME(S>-2)/60S 

IF DECLDIST(S) GTR LE THEN OECLDIST (S) =LE* 

AIRDIST(S)= (1-DECLDIST(S)/LE)*AIRDISTCS)S 
TIMEtS)=AIRDlST(S)/TASQRT*36000$ 

LE=LE-DECLD 1ST ( S > $ 

ENDS COMMENT dECLIMBS 


PROCEDURE LINE<N1»N2)5 BOOLEAN Nl»N2S 
COMMENT PRINT-OUT PROCEDURE* 

BEGIN 

INTEGER I»J»K»L»M»N»P»Q»R»S»T S 

FORMAT Pl('TOC*» 2I5» Xl2» I4» XlO» I4» 3l6» I8» A2)» 

P2(I3» X22» I4» X10» I4» 3I6» I8» A2)» 

P3(’T0D*» 4I5» I6» I8» I6» 3I6» I8» A2)» 

P4(I3» 4I5» I6» I8» 4I6» I8» A2)S 
LOCAL label LAB30»LAB31»LAB32S- 
IF N1 THEfJ goto LA830S 
IF N2 THEN GOTO LAB31S 
I=G3* S=MOS T=B3S 
G£0MGRID(G4» G3)S 
J=HEADING(G)S 
K= 310 + (G-1)*40S 
IF NOT N1 and not N2 THEN BEGIN 
L= TDEV‘6)S 

M.= SQRTUAX(G)*»2 + TAY(G)**2)S 
ENDS 

N= (IF N1 OR N2 THEN .75 ELSE !)♦ 

(WIX(6)»E1 + WIY(G)*E2)$ 

M3= IF N1 OR N2 THEN 

DECLtIME(G)*600 ELSE TIME(G)* 

P= 100*ENTIER(M3/36000) + M0D(M3» 36000)/600* 

G= 100*ENTIE.R(FLTIME/36000) + MOD(FLTIME» 36000)/600S 
R= BURN - GRWQS 
IF N1 THEN 
BEGIN 

WRITE(P1» J»K»N»DISTANCE*P»0»R»GRWQ)S 

GOTO LAB32 
ends 

IF N2 THEN 
BEGIN 

WRITE (P2.ST1»N»DISTANCE»P»0*R»GRWQ)S 



**♦*864 

***♦865 

****866 

***♦367 

****868 

****869 

****870 

****871 

****872 

**♦♦073 

****874 

***♦875 

****876 

****077 
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GOTO LAB32 
ENDS 

IF S EQL T then 

WRITE (P3»J»K»L»M»N»LE»DIGTANCE»P»0»R»GRWQ) 

ELSE TfjP‘*»I»'J»K»L»M»M»LE»OlSTANCE»P»Q»R»GRWQ)S 

LAB32.. 

ENDS comment LINES 

procedure SPAC£0PT(FBQ»FB»BI»BG»RR)S integer fb»rrs 
boolean fbo»bi»bgs 

COMMENT FBQ DETERMINES WETHER THE FLTGHTPLAN 
COMPUTATION WILL BE PERFORMED BACKWARDS (FALSE ) OR FOR 
WARDS ( TRUE ) • 

FB DETERMINES wETHER C0STS<+1) OR FUEL(O) OR FLIGHT 
TIME(-I) will bE OPTIMIZED. RR IS TAKE OFF WEIGHT OR 

landingweight. 

BI defines WEThER the navigation regime is free in the 
HORIZONTAL^FALsE) or bounded BY one point (TRUE). 

BG defines WEThER THE CRUISING ALTITUDE IS FREE (FALSE) 

OR bounded ( TRUE 

BEGIN 

LOCAL label ITER» ENDG» CI» REPEAT. EINDS 
INTEGER AA»oUANT»II»JJ»GG»GRWP» DDG» DTG» TMGS 

boolean iteration, bclimb.bno.bns 

INTEGER array ROW.ROWQ. QUANTO (0 . .F (Q) )S 
PROCEDURE PREP(A1.B1.AA.BB)S INTEGER Al.BlS 

comment Computation of distanc-e and flight time, and 
print-out of The flight plans 
boolean aa.bbs 

comment aa refers to climb, bb to descent* 

BEGIN 

DISTANCES DISTANCE+AlS FlTIME= FLTIME+BlS 
line(aa.cb)s burns grwo 

ENDS comment PREPS 

procedure editings 
begin local label ended* 

FORMAT Fo3(X44. »CH*. 14. 16. D6.2. A2>* 

PROCEDURE OPS 
BEGIN 

IF FB EuL 1 then 
WRITE ( ’COST' ) 

ELSE IF FB EQL 0 THEN 
WRITE ( ’FUEL' ) 

ELSE IF FB EQL -1 THEN 
wRITt(*FLTlME*)S 
GOTO ENDECS 
ENDS comment QPS 


WRITEC 


WRITE(PEjT)* 

DETAILED FLIGHT PLAN ♦♦**)* 
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write** ♦>» WRITE** *>s 

**♦♦972 i 


IF ROUTE EOL 1111 THEN 

****973 


BEGIN 

**#*974 


WRITE (’NON SPECIAL MIN TRACK »)S 

****975 ^ 


QP 

****976 ; 


ENOS 

****977 


IF ROUTE EQL THEN 

****978 


BEGIN 

****979 


WRITE* *ATC I*)S OPS 

w***9ao 


ENOS 

*«W4>981 


IF ROUTE EQL 5555 THEN 

*♦**982 


BEGIN 

****963 


WRITE* *ATC II* )S OPS 

***•984 


ENOS 

****985 


IF ROUTE EOL 7777 THEN 

****986 


BEGIN 

*♦**987 


WRITE* *GRC*)S OPS 

*•*•986 


ENOS 

***•989 


IF ROUTE EQL 9999 THEN 

***•990 ; 


BEGIN 

*•**991 i 


WRITE* ’SPECIAL TRACK* )S OPS 

**••992 


ENOS 

****993 


WRITE* *MIN TRACK IN SPACE* )S OPS 

*•*•994 

ended « t 

****995 

WRITE*F03* FLUR» 0ATE» TT>S 

****996 


WRITE** *)S 

••••997 

WRITeC 

• •••998 

*N0« HEAD 

FL TMP TAS WIND OIST ACCO- TIME ACCT BURN WEIGHT 

#*•••999 

)S 

*♦♦♦1000 

ENOS comment editings 

*♦♦♦1001 

*♦♦♦1002 


PROCEDURE QQ*0)S INTEGER OS 

*♦♦♦1003 


comment computation of SEGMENT CONTRIBUTION. COST» 

♦•••1004 


TIME* FUEL* WEIGHT ETC.S 

*•**1005 


BEGIN 

****1006 


LOCAL label JM1» JM2» JM3* JM* ENOQS 

*♦**1007 

*♦•♦1008 


PROCEDURE CL*M1»M2)S REAL M1»M2S 

*♦♦♦1009 


COMMENT N*.l DENOTES SEGMENT FUEL# M2 SEGMENT TIMES 

*♦♦♦1010 


BEGIN 

***♦1011 


LOCAL LABEL JUS REAL MMS 

♦♦*♦1012 


IF bclimb then 

*♦♦♦1013 


begin 

*♦*♦1014 


MM: OS GOTO JU 

*♦♦♦1015 


ends 

*♦♦*1016 


IF iteration then 

*♦♦*1017 


BEGIN 

*♦**1018 


MM= IF G GEO GG THEN -28 

*♦*♦1019 


ELSE -700S GOTO JU 

*♦♦♦1020 


ends 

*♦♦*1021 


MM= if *0 EOL -1 AND G GEO GG) OR (0 EOL 1 AND 

*♦♦*1022 


G LSS GG) THEN 700 ELSE 28S 

♦♦*♦1023 

J0*« 

MM= Ml+ 0**GG-G)* *IF MQ NEQ B1 THEn MM ELSE 0>S 

**♦*1024 

***♦1025 


JMl** 


JM2« • 


j Jm 3 «« 
I 



GRWQ: 6RWQ 0*MMS 

GUANTs QUANT+(IF FB EOL I THEN 1083»3*M2+15*MM 
ELSE IF FB EQL 0 THEN MM ELSE 600*M2>S 
ENDS comment CLS 

quant : oUANTQiDS 

iteration^ if not FBQ and MQ EOL 0 AND 0 EQL -1 
THEN true 

ELSE FALSES if ROWd) EQL 0 THEN GOTO JMS 
GRWQ= IF ITERATION THEN ROW(O) ELSE ROw(I)S 
IF MG EOL B1 OR ITERATION THEN 
BEGIN 

IF DECLFUEL(6I LSS 0 THEN GOTO JM* 

BCLIMB= trues CL(DECLFUEL(6) » DECLTlME(G) )S 
BCLIMB= FALSES 

DDG= DECL0IST(6)S 0TG= DECLTIME<G)*600 S 
IF ED THEN PREP<DDG» OTG» TRUE» FALSEIS 
ENDS 

IF NOT FBQ then GOTO JM2S 

AI04= IF TOEV(G) LSS -10 THEN 64.08+6 
ELSE TABLE<1»6»MW)S 
IF 6RWQ GEQ AID4 THEN GOTO JMS 
IF NOT FBQ then goto JM3S 

GRWP= GRWQS GRWQ= IGRWQ + 0* AIRDIST(G)/ 2 * 

(IF G EQL 1 then 2.{3'66+0; 0000722* 6RWQ ELSE 
IF G EOL 2 THEN 0.t|5 ♦ 0;0000792* 6RWQ 
ELSE -0re3+ 0.0000&9* 6RW0)S 
AA= table (0»6»RANGE)S 6RW0= GRWPS 
IF AA LSS 0 THEN GOTO JMS 

CL(AIRDIsT(6)*10**4/AA»TIME(G)/600)S 

TMG= TIMe(G>5 

IF ED then PREP(LE» TMG » FALSE* FALSE >S 
IF NOT FBQ THEN GOTO JMIS 

IF MQ EQL B3 AND NOT ITERATION THEN 
BEGIN 

BCLIMB= 

TRUES cl(oeclfuel(g)»oecltime (g))$ 

BCLIMB= FALSES 

DDG= dECLDIST(G)S DTG= DECLTIME (G) *600 S 
IF ED then PREP(DDG» DTG» FALSE* TRU£)S 
ENDS 

GOTO ENDqS 

GRWQ= os QUANT= 64.08+6S 
ENDS 

COMMENT QQS 


ED= IF BI and BG then TRUE ELSE FALSES 
FOR AA= 0 step 1 UNTIL F(Q> DO 
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****1032 
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««t**i040 

****1041 
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*♦*♦1044 

*♦*♦1045 
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. ****io68 
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♦♦♦♦1074 
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*♦♦♦1076 
*♦♦♦1077 
*♦♦♦1078 
**♦♦1079 


[ 


1 


ROW(AA)s ROWQ<AA)= 0* QUANTQ(AA)= 64»0A+6 
ENOS 

quants os bclimbs false* 

IF FBQ THEN 
BEGIN 

R0M(0)s rRS QUANTQ(0)= 0 

END 

ELSE 

BEGIN 

ROW(F(Q))s RRS QUANTQ(F(Q))s 0 
ENDS 

IF ED THEN 

BEGIN ^ 

distances FLTIMEs 0* BURNS ROW(O)* 

EDITINGS WRITE(F04# RR)* 

ENDS 

IF FBQ THEN 
BEGIN 

Bl= OS B2= 1* B3= Q-1 

END 

ELSE 

BEGIN 

BIS Q-l$ B2= -1* B3S 0 
ENOS 

COMMENT ZONE CYCLE BEGINS* 

FOR MQS 81 STEP B2 UNTIL B3 DO* 

BEGIN 

Gl= IF FBO then MQ+l ELSE MO* 

G2S IF FBQ THEN MQ ELSE MQ+lS 
M= I0*MQ+61Q+(IF 10 EQL -I THEN -I ELSE 0)5 

COMMENT A CYCLE FOR GRAPH POINTS ALONG MERIDIAN 
begins* 

FOR JJs sTOREI(GI) STEP I UNTIL STORElQtGl) DO 
BEGIN 

(]3s SQ^G1*^J)S Us UQ(Gl^UU)* 

COMMENT A CYCLE FOR GrAPH POINTS 

ALONG NEXT MERIDIAN BEGINS* 

FOR Il= ST0REU62) STEP I UNTIL ST0R£IQ(G2) DO 

begin 

IF NOT BI THEN 

begin 

IF(M GEQ 3 AND M LEO 8) AND ABS 
(II-JJ) GEO 10 THEN GOTO CIS 
BNQS (FBQ AND 10 EQL -D OR <NOT FBQ AND 
10 EQL DS 

BN= (FBQ and 10 EQL 1) OR (NOT FBQ AND 
10 EQL -DS 
IF M LSS 3 THEN 
BEGIN 

if(bno and BK 

(M»JU»II)) OR (BN AND BK(M»II» 


♦•#•1080 
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****1064 
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*♦♦♦1132 
•**♦1133 


J 


JJ)) THEN GOTO Cl 
ENOS 

IF M GTR 8 THEN 
BEGIN 

IF (BNQ AND 

BK0(M»JJ*1D) OR (BN AND BKQ 
(M*I1»JU)) THEN GOTO Cl 

END 

END* 

G4= SQ(G2*II)S 1= J0(G2»II)S 
IF FBQ THEN G£0MGRI0(G4*G3)ELSE 
6e0M6RID(63»G4)S 
PART2GEOMS 

IF m NEQ B1 THEN GG=ST0REG(MQ - B2)S 

comment a cycle for flight levels begins* 

FOR G= STOREG(MQ) STEP 1 UNTIL 

stoRegq(mq) do 
begin 

IF NOT FBQ AND MQ EOL 0 THEN 
BEGIN 

TEL = 0* IF G NEQ 6G THEN 
GOTO EN06* 

Nl= ROW(O)* ROW(0)= GRWQ= 
MAXTOW-20000S 

TEL = TEL +1* IF 6 LSS 1 

THEN 

BEGIN 

6= G6* 6RWQ= 0$ 

GOTO ENDG 
END* 

METPROC(G)* DESCLIMB(6)S 
IF DECLFUEL(G) LSS 0 THEN 
BEGIN 

TEL= 0* G= G-IS GOTO REPEAT 
END* 

QQ(-D* IF GRwO EOL 0 THEN GOTO EN06S 
N2= ROWd) - GRWQS if N2 GTR 10 
AND ROW(O) EOL MAXTOW THEN GOTO 
ENDG* 

ROW(0)= ROW(O) ♦ (IF TEL GTR lO 

THEN N2/2 ELSE N2>S 

GRWO= ROW(0)S IF ABS(N2) LSS 10 

OR TEL GTR 10 THEN 

BEGIN 

R0W(0)= NIS GOTO ITER 

END 

ELSE GOTO REPEAT* 


END 

ELSE 

BEGIN 

METPROC(G)* 
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IF MQ EOL 0 OR MQ CQL 0-1 THEN 
BEGIN 

6RWQ= RRS OESCLIMB(G) 

ENDS 

IF FBQ THEN QQ(-1> ELSE 00^1 )S 
IF 6RW0 EOL 0 THEN GOTO ENQG 
ENDS 

IF QUANT LEO QUANTO (U) THEN 
BEGIN 

ROW<J): GRWQS quanto (J) s QUANTS 
ROWQ(J)= IS IF BI AND NOT BG 
THEN ST0RC6 <MQ)s GS 
IF ITERATION THEN G= GG 
ENDS 

END* COMMENT FLIGHT-LEVEL CYCLES 

ENDS COMMENT ALONG-NEXt-MERIDIAN CYCLES 

ENDS comment ALONG-CURRENT-MERIDIAN CYCLES 

IF BI AND not BG THEN 

STOREGQtMO) = STOREG(MQ) 

ENDS comment ZONE CYCLES 


IF ROW(O) EOL 0 THEN GOTO EINDS 

IF BI then goto EINDS J=IF FBQ THEN FtQ) ELSE F(0)S 
FOR MQ= B3 STEP “B2 UNTIL BI DO 
BEGIN 

K= IF FbQ then MQ ELSE MQ+IS 1= R0WQ(J)5 
STOREIQ(mQ ^tlF B3 EOL 0 THEN 1 ELSE Q))= I-FCk-D^ 
ST0REI(K)-1S J= IS 
ENDS 

FOR MQ= 0 step 1 UNTIL Q DO STOREI(MQ)= STOREIO(MQ)S 

TOWQ= ROW(0)S Ml= IF FBQ THEN QUANTO <F(QM ELSE 
GUANTQ(0)$ 

Lw= row(f(q))s if ED Then 
BEGIN 

WRITEi* *1$ 

WRITE(F05» TOWQ-LW)S 

WRITE(F06» IF FB EQL 1 THEN Ml/100 ELSE 
( (TOWQ-Lw>*15^1083.3*FLTIME/600)/100)$ 

ENOS 

ENOS COMMENT SPACEOPTS 

PROCEDURE FF(UU»VV)S INTEGER ARRAY UU»VV$ 

begin 

F(-l)= -IS for K= 0 STEP 1 UNTIL Q DO 
F(K)= FtK-l)+VV(K)-UU(K)+l 
ENDS COMMENT FFS 

PROCEDURE PR0CES(S1*S2»S3»PR»FFBB»BB»FACT0RISATI0N)S 
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♦♦♦♦122V 

♦♦*♦1230 

**♦♦1231 

**♦♦1232 

♦*♦*1233 

**♦*1234 

♦**♦1235 

♦♦*♦1236 

♦♦♦♦1237 

*♦♦♦1238 

**♦♦1239 

**♦♦1240 

**♦♦1241 


J 




INTEGER S1*S2»s3»FFBB$ REAL PRS 

boolean bb*factorisations 

comment si and S2 determine the limits of flight 

LEVELS. 

S3 IS A DUMMY *E)(CEPT WHEN 

factorisation IS TRUE* 

then S3 IS A datum FOR THE FLIGHT 

level 

IN WHICH THE M»F»P» IS COMPUTED* 

PR IS percentage trip fuel 
ffbb determines wether COSTSC+D. 

FUEL(O) OR TIME(-l) IS OPTIMIZED. 

BB..USE LANDING WEIGHT (TRUE) * USE 
TAKE OFF weight (FALSE) 

FACTORISATION(TRUE) INITIALISES 
PROCES IN THE HORIZONTAL ON THE 
BASIS OF FLIGHT TIME. 

FOLLOWED BY COMPILATION OF FLIGHT 
PLAN IN THE VERTICAL ON THE BASIS 
OF COSTS. FUEL OR TIMES 
begin 

INTEGER NN* 

LOCAL label ENDPR*ING$ 

PROCEDURE AaS 

COMMENT specification OF LIMITS IN HORIZONTAL OF 
GRAPH POINTS ALONG MERIDIANS 

BEGIN 

LOCAL label UTRECHTS 

STOREI(0)=STOREIQ(0)=ST - A(GlQ-l)* 

STOREI(Q)= ST0REI0(0)= STl-A(G2Q-l)S 
IF ROUTE GEQ 1111 THEN GOTO UTRECHTS 
FOR Ks 1 STEP 1 UNTIL Q-1 DO 
BEGIN 

STOREl(K)s 1$ ST0REIQ(K)= E(K) 

ENOS 

FF(ST0REI»ST0REIQ) 
ends comment AAS 

PROCEDURE AAA(F1.F2)S INTEGER F1.F2S 

comment specification of limits IN THE VERTICALS 
FOR K= 0 STEP 1 UNTIL 0“1 DO 
BEGIN 

STOREG(K) = FIS STOREGO(K) = F2S 
ENOS comment AAAS 

PROCEDURE EpS 

IF TOWQ GEO MAXTOW OR TOWQ EOL 0 THEN 
BEGIN 

WRITE (‘SUPER CRITICAL* )S 
GOTO E^)0PR* 
ends comment EPS 


INTEGER PROCEDURE SSS 

SS= GRW+RESERVE^PR* ( ( GRW+RESERVE > • 

(MINTXME*0 . oe-0 . 15) -500*M1NT1ME+7300+TAXI ) S 

INTEGER PROCEDURE SSSS 

SSSs GR)i»+RESER VE+PR* ( TOWQ-LW+TAX I ) S 

INTEGER PROCEDURE SSSSS 

SSSSs TOWO* ( GRW^RESERVE+PR*TAx I ) / < ( 1+PR ) ♦Nn-PR*TOWQ) $ 

procedure nnq<a)s value as integer as 
NN= IF NOT BB then TOW ELSE As 

PROCEDURE TW* 

COMMENT safeguarding AGAINST OVERLOADINGS 
IF TOWQ GEQ MAXTOW OR TOWQ EQL 0 THEN 
BEGIN 

WRITE(PEJT)S WRITECDECR LANDING W*)S 
GOTO ENOPR 
ENOS 


• • 


S • 


COMMENT IF FACTORIZATION IS TRUE» THEN THE OPTIMIZATION 
TAKES PLACE FIRST IN THE HORIZONTAL AND THEN 
IS FOLLOWED BY IN THE VERTICALS 
IF FACTORISATION THEN 
BEGIN 

AAS AAA(s 3»S3)$ NNQ(SS)S R«UTE= IIHS 
SPACEOPTdF BB THEN FALSE ELSE TRUE *-l* FALSE* 
TRUE»NN)S TWS 

aaa(si*S2)g goto ING 

END 

ELSE 

BEGIN 

AAS AAA(S1*S2)S NN0(SS)S 

SPACEOPTdF BB THEN FALSE ELSE TRUE »FF8 b* FALSE* 
false *NN)S TWS 

EPS NNQ(S3S)S 

SPACEOPT (IF BB THEN FALSE ELSE TRUE*FFbB* TRUE* 
FALSE»NN)S TWS 
EPS NNQ(SSSS)S 

SPACEOPT { TRUE ♦ FFBB » TRUE » TRUE * NN ) S EPS 
IF factorisation THEN ROUTE* 1000 
ENDS 

ENDS comment process 


IF ROUTE GEQ 11 11 THEN 
BEGIN 

READ(PCF(»DC10F»>* ST* ST1» DONE* ERR)S 

G1Q= Z0NEKST)S 

FOR K= 1 STEP 1 UNTIL 0-1 DO 

begin READ(pCF(*DC10F»)* STOREKK)* DONE* ERR)* 


♦***1296 

♦♦•*1297 

♦♦♦♦1296 

♦♦♦♦1299 

♦♦♦♦1300 

♦♦♦♦1301 

♦♦♦♦1302 

♦♦♦♦1303 

♦♦♦♦1304 

♦♦♦♦1305 

♦♦♦♦1306 

♦♦♦♦1307 

♦♦♦♦1308 

♦♦♦♦1309 

****1310 

♦♦♦♦1311 

♦♦♦♦1312 

♦♦♦*1313 

♦♦♦♦1314 

♦♦♦♦1315 

♦*♦♦1316 

♦♦♦♦1317 

♦♦♦♦1318 

♦♦♦♦1319 

*♦♦♦1320 

♦♦♦♦1321 

♦♦♦♦1322 

♦♦♦♦1323 

*♦♦♦1324 

♦♦♦♦1325 

♦♦♦♦1326 

♦♦♦♦1327 

♦♦♦*1328 

♦♦♦♦1329 

♦♦♦♦1330 

♦♦♦♦1331 

♦♦♦♦1332 

♦♦♦♦1333 

♦♦♦*1334 

♦♦♦*1335 

♦♦*♦1336 

♦♦*♦1337 

♦♦♦♦1338 

♦♦♦♦1339 

♦♦♦♦1340 

♦♦♦♦1341 

♦♦**1342 

♦*♦♦1343 

♦♦♦*1344 

♦♦♦♦1345 

♦♦♦♦1346 

♦♦♦♦1347 

♦♦♦♦1348 

♦♦♦♦1349 


WRITE(F02» 1* STOREKkMS 
STORE! (K)s sTOREIQ(K)s STORE! (K)-A(61Q+I0*K-1)* 

ENOS 

G2Q s ZONE I ( STDS 
WR!TE(F02» ST*STl>S 
ENOS 

GE0MGR!0(ST»STDS 

N ! NT !ME:GE00 ! ST ( L A1 • LOl * lA2 » L02 ) /465S 
CTQS 

E( 0 )s E(0)s !S FOR Ks l STEP 1 UNT!L 0-1 00 
E(K)s A(G1Q+I0*K)-A(G1Q+!0*K-1)S 

COMMENT 

1. WHEN A FL!6HT PLAN FOR THE cPT!MAL COST TRACK !S NEEOEO 
CALL *PR0CES(1» 3» 1* 0*03» 1» TRUE» FALSE)* 

2. WHEN A FL!6HT PLAN FOR THE 0PT!MAL FUEL TRACK NEEOEO 
CALL »PR0CES(1» 3» 1» 0.03* 0» TRUE* FALSE)* 

3. WHEN A flight PLAN FOR THE LEAST TIME TRACK NEEOEO 
CALL *PR0CES(1» 3» l» 0*03*-l* TRUE* FALSE)** 


PROCESU* Zt it 0.03» 
PROCESd* 3* It 0.03t 


0* TRUE* false )S 
-l» TRUEf FALSE)* 


GOTO RUN 
ENOS 

OONEtf end* 

WRITE(* *)S wR!TE(* ♦)* 

WRITE (* ENO OF RUN ♦♦♦*)* 

GO TO EOPS 


CDC8). 


WRITE {* ♦** INPUT OATA ERROR *♦♦*)* 

aSlT*0C6* 1*810528# 065734 C DC 8). 

comment KNmI/KLM 020169 SPG NAV « JNG* V 

NAVIGATION-FLIGHTPlANNING-MOOULE for the proouction of 
FLIGHTPLANS over the NORTH ATLANTIC* 

THE MAIN features 

are: 

1. inclusion of 6 meteorological parameter fields* 
temperature and GEOPOTEnTIAL fields for 300 250 AND 200 MB 
valid for two standard times 12 HOURS APART* 

2* NAVIGATIONAL GRID BASED ON REPORTING POINTS OF THE KLM* 

3* INCLUSION OF A DYNAMIC PROCESS BASED ON INTERPOLATION OF 
TWO PARAMETER FIELDS* 

4* OPTIMALISATION OF FUEL* COSTS OR FLIGHT TIME IN SPACE AS 
well as in the horizontal as in THE VERTICAL* 

5. IN - AND OUTBOUND TRAFFIC (10= 1 OUTBOUND 10= -1 INBOUND)* 
6* FLIGHTPLAN PRODUCTION FOR OPTIMUM TRACKS* SPECIFIC ROUTES 
ETC* E*G* ALTERNATIVE ATC MINIMUM ROUTES* 

7* BLOCKING E*G* SECTOR BLOCKING OVER THE NORTH ATLANTIC 
BLOCKING In AIRWAYS AND ATC RESTRICTED AREAS* 

6* STANDARD CRUISE AND USE OF PERFORMANCE TABLES* 


4***1350 
*••*1351 
**••1352 
*•••1353 
*•••1354 
**••1355 
***•1356 
••••1357 
***•1356 
*•••1359 
*•••1360 
*•••1361 
*••*1362 
SIEEOEO ****1363 
***•1364 
DEO *••♦1365 

*•••1366 
D ****1367 

*•••1368 
*♦•♦1369 
*•♦•1370 
*•♦•1371 
*•••1372 
*•*•1373 
*•••1374 
*•••1375 
*•••1376 
*•••1377 
*••♦1378 
*♦♦*1379 
*•••1380 
**••1381 
AC*>WKWAACa3a aX 
♦♦♦♦1 
♦♦♦♦2 
♦♦♦•3 
•♦♦♦4 
*♦♦♦5 
*♦♦♦6 
•♦♦♦7 
♦♦♦♦8 
♦♦♦♦9 
♦♦*♦10 
•♦♦♦11 
*♦••12 
•♦•♦13 
♦♦♦♦14 

•♦♦♦lb 

*••*16 

*♦•♦17 

•♦♦♦18 

••♦♦19 

♦••♦20 

*♦♦♦21 


constant KQ= 14 » W= 3$ 


INTEGER P1»P2»P3»01*02»Q3»S»TTT* 
61*G2»G3*64*Bi*B2*B3*MO*GlO*G2Q»Q» 
RESl»RES2*RES3»l»J»fU»N2»N3»M»G» 
M0*Ml»M2fM3*M4*MS*R6*K*10* 

KI»K2»L1»L2»K10»K2<J»L1Q»U20*KK1»KK2»LL1»L1.2» 

ENOU » TEM » T AS » FLUR » T Ax I • TO W » GRW » RESER VE » MAXTOW » MAXLW • 

DATE » ST » ST 1 » ROUTE » 0 1 stance » FLT X ME » burn » 

GR WQ • TOWO » LW * COST * LE * TELS 

REAL C0»CG»P»C1»X1»X2»Y1»Y2»X»Y»EI»E2»E3» 
GG1»GG2»0»ANGLE»LALAt»LAL 0NG»TT»LA1 »LOl »LA2»L02f 
XX1»XX2»YYI*YY2»AI01»AI02*AID3»AI04»XSTER»YSTER» 

MINTIMES 

INTEGER ARRAY HT«1..72»4. .19) » A(-l..KQ+l)» V(0..4l)» 
SR(1..3*0..19) > MW(1.»3»0*.6) » CO*CT»CF(1..3»0..17) * CTC» 
CFC(0«.7*0..8) • CDC(0**B*0..6)S 
ARRAY OD(0.*KO)Y> 

BOOlIaN array BK*0..2'1«*6»1«.6)* BKq(9* .KQ»1. *6»1. .6)S 

FORMAT F0l(E2» S10» O7.0» A5)* 

PEJT <E2» • •» Al)» 

F0T(X5* •I6» a2)» 

F02U5* *I7» A2)S 
STRING STR<X0)S 


INTEGER PROCEDURE INTERP0L(F»F1»F2»F10»F2Q)S 
INTEGER F»F1»F2*F10»F2Q* 

COMMENT LINEAR INTERPOLATIONS 

INTERPOL s <<F- FI) ♦ F2Q ♦ (F2 - FI ♦ FXO) / 


(F2 - Fl)$ 


REAL PROCEDURE SUMKA» B» ARR)S 
VALUE A. BS 

INTEGER A» BS ARRAY ARRS 
BEGIN INTEGER IS R£AL SS 
S= OS 

FOR 1= A STEP 1 until B DO 
S= S ARRtDS 
SUMl = SS 

END* COMMENT SUMIS 

REAL PROCEDURE SUM2(A» B» ARRl# ARR2)S 
VALUE A» BS 

INTEGER A» BS ARRAY ARR1» ARR2S 
BEGIN INTEGER IS REAL SS 
SS OS 

FOR 1= A STEP I until B DO 
S = S ♦ ARR1(I)*ARR2(X)S 
SUM2 = SS 

ENDS COMMENT SUM2S 


PROCEDURE WW(AA» BB) S REAL AaS STRING BBS 
MRXTEtFOl* BB* AA)S 


4**«22 
*•*•23 
*«*«24 
****2b 
***•26 
• #**27 
****28 
*••*29 

*•••30 

****31 

****32 

*•••33 

♦***34 

♦***3b 

***♦36 

*•••37 

*•••36 

*•••39 

••••40 

*•••41 

*•••42 

****43 

*•••44 

****4b 

*••*46 

***•47 

••••46 

••••49 

*••*50 

*•••51 

*•••52 

••••53 

••••54 

••••55 
*•••56 
*•••57 
••••56 
••••59 
••••60 
••••61 
••••62 
••••63 
••••64 
••••65 
*•••66 
•••♦67 
••••68 
••••69 
*•••70 
*•••71 
••••72 
••••73 
• •••74 
*•••75 


PI S 31000S P2 2 35000S P3 s 39000S 01 = 30065* •♦••76 

02 = 33999* 03 = 38662* TCL = 0* ♦•••77 

FOR Ss 0 STEP 1 UNTIU 5 00 ••♦♦78 

BEGIN ••••-79 

RESls S^12» ••♦•80 

FOR JS19 STEP •! UNTIL 4 00 ♦♦••81 

FOR IS RESl ♦ 1 step 1 UNTIL •♦••82 

RESl ♦ 12 DO READ(PCF(*NASAO»)» HT<1»J)» ERR» ERR)* ♦•••83 

END* ♦•••84 

FOR S=0 STEP I UNTIL 5 DO ••••85 

BEGIN RESl =S^12* WRITE (PEJT)» ••••86 

FOP JS19 STEP -i UNTIL 4 DO ♦•♦•87 

WRITE(F02* 12» FOR ISRESI ♦ 1 ••••88 

STEP 1 UNTIL RESl • 12 DO HT<I*J))* ♦•••89 

end* ♦•••90 

••••91 

COMMENT ARRAY *HT* 1$ REFILLED KITH INTERPOLATED PARAMETERS ••••92 

WHICH APPLY IN GRID POINTS AT FLIGHT LEVEL* 31000* ••♦•93 

35000* AND 39000 FT* ••••94 

FOR IS 1 STEP 1 UNTIL 12* 37 STEP 1 UNTIL 48 DO ••••95 

FOR J= 4 STEP 1 UNTIL 19 DO ••••96 

BEGIN ♦•♦♦97 

MOS HT(1*J)* ♦♦••98 

MIS ENTIER(MO/100)* M2= M0-100*M1* ••♦♦99 

MOS HT(I+12*J)* ♦♦••JOO 

M3= ENTIER(MO/100)* M4S M0-100^M3* ♦••♦101 

M0= HT(I>24*J)* ♦•••102 

MSS ENTIER<M0/100)* M6= M0-100^M5* -- ••••103 

HT<I*J) s 100^INTERPOL<P1*01*G2*M1»M3)^ ••••lOU 

INTERPOL (PI *01 *02 »M2»M4)» ••••105 

HT<I^12*J>= 100^InTERPOL(P2»Q2*03*M3*M5)'* ••••106 

INTERPOL ( P2 *02* 03 *M4*M6)* ••••107 

HT(I+24»J)s 100*lNTERPOL(P3*02*Q3fM3*M5)+ ••♦♦108 

INTERP0L(P3*02*03#M4*M6) ♦•♦♦109 

END* ♦•••110 

COMMENT THIS ROUTINE CONCLUDES THE PRODUCTION OF THE ♦♦♦♦111 

meteorological parameters geopotential and TEMPERATURE IN THE ♦♦••112 

FLIGHT levels* ♦♦•♦113 

♦♦♦♦114 


COMMENT relationship BETWEEN GRAPH POINT NUMBERS ANq ZONE 
A(-l)s -IS A(0)= 0* A(l)s 2* A(2)= 6* A(3>= 11* A(4)= 26* 
A(5)s 41* A(6)= 56* A(7)s 71* A(8)s 86* A<9)s 90* A(10)s 96* 
A(ll)s 102* A(l2)s io7* A<13)s 111* A(14>= 114* A(15)s 115* 
V(0)s 52290477* 

53640150* 


V(l)s 

V(3) = 

V(6)s 

V<9) = 

V(12)s 

V(15)s 

V(18)S 

V(2l) = 

V(24)s 

V(27)s 


51850115* 

51000200* 

55000200* 

52700892* 

54271005* 

53735697* 

48545856* 

53286035* 

48006000* 

51006750* 


V(2) = 
V(4) s 
V(7) s 
V(10)= 
V(13)= 
V(16)= 
V(l9)= 
V(22>= 
V(25)= 
V(2fe)s 


51990539* 

49980632* 

53500630* 

48905454* 

584T6263* 

55568000* 

54836683* 

49846439* 

52006900* 


V(5) s 53580300* 
V(8) = 51840849* 
V(ll)s 55490459* 
V(14)= 51375560* 
V(l7)s 47005800* 
V(20)s 57006000* 
V(23 )s 46156006* 
V(26)s 50506500* 
V(29)s 45006300* 


••♦♦116 
♦•♦♦117 
♦•♦•118 
•♦•*119 
••••120 
••••121 
•♦♦•122 
••••123 
•*•♦124 
•••♦125 
•••♦126 
••••127 
♦•••126 
•••♦129 


V(30)s 46176459$ V(3l)s 47506600$ V(32)s 48576626$ 

V(33)S 50007050$ V(34>s 43836606$ V(3S)s 44846667$ 

V(36)s 45327179$ V(37>- 45467385$ V(36)s 41287003$ 

V(39)s 42367099$ VUo>- 43007250$ V(4l)s 40647378$ 

COMMENT determination OF CONNECTIVITY AM0N6 GRAPH POINTS* 

falsesconnecteo* truesblockeo f 

FOR IS 1 STEP 1 UNTIL G 00 
FOR js 1 STEP 1 until 6 DO 
BEGIN 

FOR MOs 0*1*2 DO BK<NQ»1*J) s TRUE* 

FOR MQs 9 STEP 1 UNTIL KQ DO BK0(M0*I«J) s TRUE 
EN0$ 

6Ka*l*l)s BK(2*l*l)s BKU*2*2)s BK(2*2*2>s 
BKQ(9*l*l)s BKQ(10*l*l}s 6K0Cll*l*l)s BKQ( 12 *l*l>s 
SKQ(l3*l»l)s BK0U4»l*l)s BKG(9*2*2)s BKQ(10*2*2 )s 
BKQ( ll*2*2)s BKQ(12*2*2 )s BKQ<13*2*2)s BKQ<10*3*3)s 
BKQU1*3>3)s BKQ(13*3*3)s BKQ<11*4*4)s BK<2*3*4)s Bk(2*4*5)s 
BKQU2*4*4) s BKQ(10*6*6) s BK(1*1*2) s 
BK(1*2*3} s BK(1*2*4) • 8K(2*2*3) s 
BKQ(9*1*2): BK0(9*2*3’- FALSE$ 

BKQ(9*2*5)s BKQ(9*3*4 >s BKQ(9*3*5)s BKQ(9»4*6)s BKOaO*2*I): 
BKQ(10*4*3)s BKSao»5»3)s BKQ(lO*S*4)s BKO(10*6*S)s 
BKQ(11*1*2)S BK^ai*3*4)s BK0<11*5*4U FALSE$ 

BK0(ll*6*5)s BKQ(12*2a)s BK$<12*3*2): BKQ(12*4*3)S 
6K0(12*S*4)s BKC«13*l*2)s BKQ(l4*2*lis BKQ<14*3*l)s bK(0*1*1): 

BK (0*1*2) sFALSES 


BEGIN 

REAL ARRAY AB<0**10)* LEN6TH( 1 • *10) * 
TIME*AIROIST*TOEV*W1X*WIY*TAX*TAY(1..w)$ 

INTEGER procedure Z0NEK(i)$ INTEGER 0$ 

COMMENT this SUBROUTINE DETERMINES THE INDEX OF THE 
ZONE ASSOCIATED KlTH A STATI0NSNUMBER$ 

BEGIN 

LOCAL LABEL AgA1N$ INTEGER 1$ 

IS -1$ 

Again* • 

IS !♦!$ IF A(I)-0 LSS 0 THEN GOTO AGAIN$ 

ZONEI= I 

ENO$ COMMENT ZONElS 
procedure 0T(K)$ INTEGER K$ 

comment This subroutine dissects latitude and 

LONGITUDE FROM THE COMPRESSED COORDINATES* 

EAST OF GREEN4HICH THE SIGN OF LONGITUDE CHANGES$ 
BEGIN 

INTEGER K0$ 

KOS ENT1ER(V(K)/10000)$ 

LALATS K0/100$ 

LALONGs (V(K)-10000*K0)/100 
END$ COMMENT DISSECTS 


♦•♦•130 

••••131 

••••132 

••••133 

♦••♦134 

♦•♦♦135 

••••136 

♦•••137 

♦••♦138 

♦••♦139 

♦•♦•140 

♦•♦♦141 

♦•••142 

♦•••143 

♦••♦144 

♦•♦♦145 

♦•••146 

•♦♦♦147 

♦•♦♦148 

•♦•*149 

♦••♦IbO 

♦••♦151 

•••♦152 

♦♦♦♦153 

♦♦••154 

♦••♦155 

♦••♦156 

♦••♦157 

♦••♦158 

♦•♦•159 

•♦♦♦160 

♦♦♦•161 

♦♦♦•162 

♦•♦♦163 

♦•♦•164 

♦*•♦165 

♦•••166 

♦•♦•167 

•♦•♦168 

♦••♦169 

••♦•170 

••••171 

••♦•172 

♦•♦•173 

♦•••174 

••♦♦175 

♦♦♦•176 

♦•♦♦177 

•♦•♦178 

♦•••179 

••♦•180 

♦•♦♦181 

•♦••182 

♦••♦183 


procedure LXS<U)S integer us 

COMMENT DETERMINATION OF LATITUDE AND LONGITUDE FOR 
GRAPH POINTSS 

BEGIN 

IF U EOL 19 then DT(12) 
else if U leg ll THEN 

begin 

OT(U)S IF U EQL 0 OR U EGL 1 OR 
U EQL 2 THEN LALON6= -LAL0N6 

EMr 

ELSE IF U LEG 26 THEN 
BEGIN 

LALATS 3S'*‘US LALONGs 10 

END 

ELSE IF U LEG 4l THEN 
BEGIN 

LALATs 20+US LALONG= 20 

END 

ELSE IF U LEG 56 THEN 
BEGIN 

LALAT- 5+US LALONG= 30 

END 

ELSE IF U LEG ?1 THEN 

begin 

LALATS U-llS LAL0N6= 40 

END 

ELSE IF U LEG 86 THEN 
BEGIN 

LALAT= U-26S LALONGS 50 

END 

ELSE DT<U-74) 

ENDS COMMENT LISS 


procedure ctqs 

comment this SUBROUTINE PREPARES THE TIME INSTANTS AT WHICH 
IN EACH ZONE THE METEOROLOGICAL PARAMATER WILL BE 
DERIVED FROM BOTH PARAMETER FIELDS* COMPOSITE CHARTS 
IN TIME are SIMULATEDS 
BEGIN 

INTEGER J0»M0S 


real a»sum*ress 

As 0*4S 

JO = IF 10 EOL 1 THEN KG ELSE 0* SUM= OS 
FOR 1= KQ-JO step 10 UNTIL UO DO 


BEGIN 

SUMS SUM+(IF 
DD<I)= SUM 
ENDS 

MO = 20NEUST)S 

DD(M0-l)S 

FOR 1= 0 STEP 1 

Then tt-12 else 

ENDS COMMENT CTOS 


I 6TR 3 and I LSS 8 THEN 2*A ELSE A)S 

RES= IF 10 EOL I THEN DD<M0> ELSE 

UNTIL KO DO OD(I)= IF TT GTR 12 
TT+A/a+DOCP-RESS 
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REAL PROCEDURE HH(V*W»F»CH0XC6>* INTEGER F#V»WS 

commInt computation of grid point values fop geopotential 
OR temperature (Choice - true or false) s 


BEGIN 

MSS (F-l)*12$ 

M6S (F<^2)*12S 
M3= HTtV^MStV)! 

M4S HT(V^M6»V)S 
MIS ENTIER(M3/100)S 
M2= ENTlER(M4/i00)S 
IF NOT CHOICE THEN 
BEGIN 

Ml= M3-100*MI* 
M2S M4*100 *m2 


Ends 

hhs 


{(12-DD(M) )*M1’»‘00(M)*M2)/12 


ENOS COMMENT HHS 


REAL PROCEDURE 6E0P(V»W»C»CH0ICE)S INTEGER CS 
REAL V*WS BOOLEAN CHOICES 

comment computation OF geopotential or temperature IN 
ARBITRARY POINTS 


begin . 

INTEGER VI»W1S REAL A»B»Bl*B2»B3»B4S 

VIS ENTIER(V)S Wis ENTIER(W)S 

Bis IF CHOICE THEN HH( Vl»Wl»C#TRUE) 

ELSE HH(Vl»Wl»C»FALSE)S 

B2S IF Choice then hh(vi^i»wi»c»true) 

ELSE HH(Vl+l»Wl»C»FALSE)S 

B3= IF CHOICE THEN HH(V1-H»WI+1»C»TRUE) 

ELSE HH(Vl<«-l»V<l+l»C»FALSE)S 

B4S IF CHOICE THEN HH(Vl»Wl'*’l»C»TRUE) 

ELSE HH(Vl»Wl+l»C»FALSE)S 
As vl^l-VS Bs wl+l-WS 

GEOps A*B*Bl+ ( l-A) *B*B2^ ( 1-A ) ♦ ( l-B ) »B3+A* ( l-B) *B4 
ENDS COMMENT GEOpS 


procedure xytransf<lat»lon)s real lat*lons 
BEGIN 

real r»ss 

Rs (45-LAT/2)*C6S 
Ss (80-L0N) * CGS 

Xs 0»7909*CO*S1N(R) *SIN(S)^C0S(R)S 
Ys (p-CO*SIN(R)*COS(S)/COS(R))*0.7909 
ENOS COMMENT XYTrANSFS 

REAL PROCEDURE GE0OIST(OOO»LAT)S REAL ODD»LATS 
comment computation of distance on the SPHERE FROM 
EUCLIDEAN DISTANCES 

GEODlSTs 219»721*DDD*C0S( (45“LAT/2)*C6)**2S 
procedure 6EOMGRIO(PP*6(»^* integer PP»GGS 
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COMMENT COMPUTATION OF EUCLIDEAN DISTANCE AND UNIT VECTORS 
BEGIN 

LIS(PP)S LAIS lALATS LOls LAL0N6S 
LIS (66) S LA2 s lALATS L02s LAL0N6S 
XYTRANSF(LA1»L0I»* Xis XS Yis YS 
XYTRANSF(LA2.L02)* X2S XS Y2S YS 

E3= S0RT((X2-X1)**2^(Y2-Y1)*'>2)S 

El= (X2-XD/E3S 
E2S (Y2-YD/E3 
ENDS COMMENT 6E0M6RIDS 


procedure PART26E0MS 

COMMENT DETERMINATION OF LENGTH OF SEGMENT 
contributions from SUBSE6MENTSS 


BY SUMMATION OF 


Kis ENTIER(X1)S K2S Kl^lS Ll= ENTIER(YI)S 
L2S Ll+lS KIQS ENTIER(X2)S K2 Qs KIQ+IS 
LIQS ENTIER(Y2)S L2Q= LIQ^IS 
Ds 6E0DIST(E3» (LA1^LA2)/2)S 
Nl= ABSCK1Q-K1)S 
N2= ABS(LIQ-Ll)* 

N3 = IF KIQ EQL K1 AND LlQ EQL Ll THEN 0 
ELSE IF Nl LEG N2 THEN N2 ELSE NIS 
BEGIN 

INTEGER I»OS 
LOCAL LABEL LABIS 
AB(0)= OS Ab(N3+1)= is 

IF N3 EQL 0 THEN GOTO LABIS If N3 EQL Nl THEN 
BEGIN 

IF Kl 6TR KIQ THEN 
BEGIN 

RES1= KIS 

RES2= K2QS RES3S -1 


END 

ELSE 

BEGIN 

RES1= K2S RES2= KlQS 
RES3= 1 


ENOS 

FOR 1= RtSl STEP RES3 UNTIL RES2 DO 
BEGIN 

J= RLS3*(I“RES1)’^1S ab(j)= (I-Xl)/ 
(X2-X1) 


END 

END 

ELSE 

BEGIN 

IF Ll GTR LlQ THEN 
BEGIN 

RES1= LIS 

RES2= L2QS RES3S -1 


END 

ELSE 

BEGIN 
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RESls L2S RCS2S LlQS 
RES3: 1 

ENOS 

FOR 1= ReSI step RES3 UNTIL RES2 00 
BEGIN 

J= RES3*(I-RES1>+1S (I-Yl)/ 

(Y2-Y1) 

ENO 

ENOS 

FOR 1= 1 STEP 1 UNTIL N3+1 00 

BEGIN ^ , 

AI02= < ( AB ( I ) ♦AB ( I-l ) ) *LA2+ (2- ( AB< I ) 
♦AB(I-1)))*LA1)/2S 

LENGTH(I)= GE00IST(E3»CAB(I)-AB(I-1))* AI02) 

ENO 

ENO 

ENOS COMMENT PART26E0MS 

PROCeOURE METPR0C(S)S INTEGER SS 

comment PROCESSING METEOROLOGICAL PARAMETERSS 

BEGIN 

integer IS 

ARRAY TEMP* WNOX • WNDY » ENDUR ( 1 • *N3+1 ) S 
G61= 6E0P(X1»Y1»S*TRUE)S 
GG2= 6E0P(X2»Y2»S*TRUE)S 

AN6LE= C1»(GG2-6G1)/(470*D*SIN( <LA1+LA2)*CG/2) )S 

XX1= XIS YY1= YlS 

FOR 1= 1 STEP 1 UNTIL N3+1 DO 

BEGIN 

XX2= AB<I)*x2+(1-AB(I))*X1S 
YY2= AB<I)*Y2+(1-AB(I) )*Y1S 

A102= ((AB(i)+AB(I-1))*LA2*(2-(AB(I)+A6(I-1)))* 
LAD/2S 

XSTER= <XX14*XX2)/2S 
YSTER= <YYUYY2)/2S 

KK1= ENTIER(XSTER)S KK2= KKl+lS LLl= ENTIER(YSTER)S 
LL2= LLl+15 

TEMP ( I ) = -GEOP ( XSTER * YSTER * S » FALSE ) S 

AID1= C1*E3/(D*SIN(AID2*CG) )S AID2= XSTER-KKIS 

A1D3= YSTER-LLIS 

WNDX(I)= -AID1*(AID2*(HH(KK2*LL2»S»TRUE)-HH 
(KK2»LL1»S»TRUE) )+<1-AID2)*(HH(KK1»LL2»S*TRUE) 
-HH(KK1»LL1»S*TRUE) ) )5 

WNDY(I)= AlDl*(AID3*(HH(KK2»LL2»S»TRUE)-HH 
(KK1*LL2 »S*tRUE) )+ tl-AlD3)*(HH(KK2»LLl*S*TRUE) 
-HH(KK1*LL1»S*TRUE)))S 

TAS= 0 . 8034*38. 9826»SQRT(273»16^TEMP(I) )S 

ENDUR < I > = 36000*LENGTH ( D / ( TAS*SQRT < 1-ANGLE**2 ) 

<WNDX<I)*E1+WNDY(I)*E2) )S 

XX1= XX2S YYl= YY2 
ENOS 

ENOU= SUMld* n3*1» ENDUR )S 

TEM= SUM2(1» N3'*'l» ENOUR* TEMP)/ENOUS 


•••♦346 I 

••••347 1 

••••348 5 

••♦•349 

••♦•3S0 

••••3bl 

•••♦352 

••••353 

•••♦354 

•••♦355 

••••356 

••••357 

••••358 

••••359 

••••360 

••••361 

••••362 

•«** 3^3 

••••364 

••••365 

••••366 

•••♦367 

••••368 

♦•*♦369 

♦••♦370 

••♦•371 

••♦•372 

♦••♦373 

♦♦♦♦374 

••♦•375 

♦♦♦•376 

••♦•377 

•••♦378 

••••379 

♦••♦380 

♦♦••381 

•♦•♦382 

♦♦♦♦383 . 

♦••♦384 ; 

*•♦♦385 

♦♦♦♦386 

♦•••387 

*•♦♦388 

♦♦•♦389 

•♦••390 

♦•••391 

*♦♦♦392 

♦♦♦♦393 

•♦♦♦394 

•♦♦♦395 

•♦♦♦396 

••♦•397 

♦♦♦♦398 

♦♦♦♦399 


-LAB20* 




TASS 0.8034*38. 9826*SQRT < 273. 16+TEM ) S 
RESls 31*(S-1)*4S 

TDEV(S)= if RESl LEO 35.332 THEN TEM +1.98*RES1-15 ELSE 
SS+TEMS 

TIME(S)= ENOUS 

AIR0IST<S)= ENDU *TAS/36000S LE= OS 
MIX(S)s SUM2U* N3'*’l« ENOUR* WNOX)/ENDUS 
WlY(S)s SUM2(1# N3^1» ENDUR* WN0Y)/ENDUS 
TAX ( S ) = T AS* ( El *SQRT ( 1-AN6LE**2 > +E2* ANGLE ) S 
TAY I S ) = T AS* ( E2*SORT ( l-AN6LE**2 > -El* ANGLE ) $ 

ENOS COMMENT METpROCS 
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****411 

COMMENT PROCEOURES ’TABLE* AND *TABLEO* ARE FOR TABLE L00K-UPS****412 
INTEGER procedure TABLE(P»R»MA)S INTEGER P»RS INTEGER ARRAY ****413 

MAS ****414 

BEGIN ****415 

INTEGER G0»DrBS REAL RO»C$ LOCAL LABEL LAB20S ****416 

0=IF P EQL 0 THEN 68000 ELSE IF P EOL 1 THEN -10 ELSE 80000S****417 


BsIF P EQL 1 then 5 ELSE 4000 S 

C= IF P EQL 1 then TDEVCR) ELSE GRWQS 

RO= (C-OJ/BS 

60S ENTIER(R0)S 

RO= RO-GOS 

IF MA(R»G0+1) EQL -1 THEN 
BEGIN 

TABLES -IS goto LAB20 
ENDS 

TABLES RO*MA ( R » 60+1) + ( 1-RO > *MA (R».60 ) S 


ENOS COMMENT TABLE* 


INTEGER PROCEDURE TABLEQ(Q»L*MA0#MA1>* INTEGER Q»LS 

INTEGER ARRAY MAO.MAIS 

BEGIN 

integer H*S0»Z»R»S00S 
LOCAL LABEL LAb21S 

real T0» toos 
H= table<q»l»mao)s 
IF H EOL -1 then 
BEGIN 

TABLEQS -is goto LAB21 
ENDS 

25 IF 0 EQL 2 then 0 ELSE IF Q EQL 3 THEN 50 ELSE 600S 
R5 IF Q EOL 2 then 50 ELSE IF Q EQL 3 THEN lOO ELSE 1000* 
TOs (H - 2)/RS 
TOOs (TDEV<L)+20>/5S 
SO= ENTIER(T0)$ 


S005 ENTIERUOO)* 

TOS TO-SO* 

TOOS TOO-SOOS 

IF MAl(SO+l»SOO’*'l) EOL -1 THEN 
BEGIN 

TABLEQ5 -IS 
goto LAB21 
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TABLEQsTO*TOO*mA 1 (SO^l »SOO+l ) ♦( 1-TO) ♦ ( 1-TOO) * 

MA1< so » SOO ) ♦ < 1-TO ) ♦TOO*MAl (SO »SOO+l ) ♦TO* < 1-TOO) *MA1 

(SO-t^ltSOOlS 

LAB21« • 

ENDS comment TABLEQS 

procedure READQ( Yo*Y1*V2*ARR)S integer Y0»Y1*Y2* 

INTEGER ARRAY ARRS 
BEGIN 

FOR 1= YO STEP 1 UNTIL Yl DO 
FOR J=0 STEP I UNTIL Y2 DO 

read(pcf(*nasao*)» arr(i« err* err* 

END* comment read* 

procedure WRITEQ(Y0* Yl* Y2* ARR)S 
INTEGER Y0*Yl*Y2S 
INTEGER ARRAY ARRS 
BEGIN 

FOR 1= YO STEP 1 UNTIL Yl DO 

WRITE < FOT* Y2 ♦ I* FOR J= 0 STEP 1 
UNTIL Y2 DO ARRd* J)»S 
ENDS COMMENT WRITEQS 

READ0<1»3»19»SRJ$ READQ(l»3»6»MW)S READQ(1*3*17*CD)S 
READQ(l»3*l7»CT)S READ0(l»3*l7*CF)$ READQ(0*7*e*CTC)S 
REAOQ(0*7*8*CFC)S REAOQ(0*6»e*CDC)S 

WRITE(PEUT)S 
WRITE0(1*3*19*SR)$ 

WRITE(» •)* WRITe(* 

WRITEQ(1*3*6»MW)S 

WRITE(* •>$ WRITeC *>* 

WRITE0(1*3*17*CD)$ 

WRITE(» ♦)» WRITE(* *** 

WRITE0(1»3»17*CT)S 
WRITE<PEJT)5 
WRITEQ(1*3*17»CF)$ 

WRITE(» WRITE<» 

WRITEQ(0»7*8»CTC)$ 

WRITE(PEUT)> 

WRITEO(0*7*8»CFC)S 
WRITE<» WRITE<* 

WRITEO(0*6*8»CDC)S 

C6= 0.01745329251 C0= 39.6296148$ Cl= 21.47*32.808399$ 
P= CO*SIN(3l.2*CG)/COS(3l.2*CG)^ 

*^^*^** REA0<PCF(»NASA0*)» FLUR* TAXI* grw. reserve* tow* 
MAXTOW*MAXLW*DATE* lO.TTT * ROUTE* DONE. ERR)* 

WRITE(* •)» WRIT£(» *)* WRITEC* ♦)* 

WRITE(* INPUT parameters *•♦•)$ 

WRITE (F02* II * FLUR'TAX' »(iRW» RESERVE* TOW* 

MAXTOW * MAXL W * DATE » 3 0 * TTT * ROUTE > » 
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IF ROUTE 6EQ 1111 THEN 
BEGIN 

READ(PCF(»NASAO»)» Qr DONE* ERR)$ 

WRITE<F02» 1* Q)* 

Q=Q-1$ 

END 

ELSE 

BEGIN 

READ(PCF(*NASAOO» ST»STl» DONE* ERR)* 

WRITE(F02* 2* ST* STD* 

GlQs 20NEI<ST)* G2Q= ZONEKSTD* <J= ABS(G2Q-G1Q) 
END* 

TT = TTT* 


Integer array storei»storeiq*storeg*storego»e(o*.q)* 

F<-1..Q)S 

IORMAT F04(X63» 16* A2)* 

F05t*TRXP FUEL** 16* Al)» 

F06<*C0ST** D8*2* A1)» 

real array OEClTIME*OECLFUEL»OECLOXST(1..W)S 
integer procedure S0(H1*H2)* integer H1*H2S 

SO= AtGlQ+I0*Hl-l)‘»-H2* 

INTEGER procedure JQ(H1»H2)S INTEGER H1*H2* 

JQ= F<H1-1)'*H2-ST0REI<H1)+1* 


PROCEDURE DESCLIMB(S)* VALUE S* INTEGER S* 

COMMENT computation OF PERFORMANCE IN CLIMB OR DESCENT* 


begin 

INTEGER ALTS 
LOCAL LABEL JUM* 

ALT= 31-MS-1)*4* 
DECLTIME(S)= IF MQ EQL 0 
/lOOO ELSE 3.1+0«41*ALT* 
DECLOIST(S)s IF MQ EQL 0 
ELSE 9+3*ALT* 

DECLFUEUS)= IF MQ EQL 0 


THEN 60nABLEQ(3»S*CT»CTC) 
then TABLEQ(2»S*C0*CDC) 
THEN TABLEQ<4»S*CF»CFC) 


ELSE 940+8. 5*ALT* IF DECLTIME(S) LSS 0 THEN GOTO JUM* 


AID4= DECLDIST ( S) *60/ (DECLTlME < S> * 
SQRT(TAX(S)**2+TAY(S)**2) )* 

OECLDIST(S)= SORT( (0.75*WIX(S)+TAX(S)*AID4)**2+ 

( 0 • 75*WI Y ( S ) +TA Y ( S ) *A 104 ) **2 ) * ( DECLTlME ( S ) -2 ) /60* 
AIROIST(S)= (1-DECL0IST(S)/LE)*AIRDIST(S)* 
TIME(S)= AIrOIST<S)/SQRT<TAX(S)**2+TAY(S)+*2)* 
36000* LE= LE-0ECLDIST<S)* 


END* COMMENT QECLIMB* 

real procedure arctano(oi»o2>* Real oi* 02* 

ARCTANQS ARCTAN(Q2/Q1)/CG+<IF 01 LSS 0 THEN 180 
ELSE IF 02 LSS 0 THEN 360 ELSE 0)* 

procedure L0<2X* 22)* REAL 21* 22 * 
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****554 

****555 

****556 

****557 

****558 

****559 

****560 

****561 


LAB30«i 


LA831** 






COMMENT computation OF LONGITUDE FROM EUCLIDEAN 
IN rectangular grids 

LOIS 80-ARCTAN(Z1/U9-22))/CGS 

Integer procedu^^e heaoing(li*L2)s real li»L2S 

COMMENT computation OF HEADING IN FLIGHTS 
BEGIN 

REAL M1»M2>ZX*ZYS 

XSTERS (Xl+xZ)/2» YSTER= (Yi+Y2)/2S 
lo(xster»yster)s mis SIN((80-L01)*CG)S 
M2S COS< t80«L01)t‘CG)S 

headings ARCTANQ(-LI*M1+L2»M2»L1*M2 ♦ L2*Ml) 
ENOS COMMENT HEADINGS 






COORDINATES ♦♦♦♦562 




****553 

♦♦♦♦564 

♦♦♦*565 

****566 

♦♦♦♦567 

****568 

♦♦♦♦569 

♦♦♦♦570 

♦♦♦♦571 

♦♦♦♦572 

♦♦♦♦573 

♦♦♦♦574 

****575 


PROCEDURE LINE<N1»N2)S BOOLEAN N1»N2S 

COMMENT PRINT-OUT PROCEDURES 

BEGIN 

INTEGER I»J»K»L*M»N»P»QiR*S*T S 

FORMAT PIC ’TOC* » 2I5» Xl2» I4» XlO» 14# 3I6» 18# A2)# 

P2CI3# X22* 14# XIO# I4» 316# 18# A2) » 

P3(*T0D*# 415# I6# 18# 16# 316# 18# A2)# 

P4CI3# 415# 16# 18# 416# 18# A2)S 
LOCAL label LAB30#LAB31#LAB32S 
IF N1 THEN goto LAB30S 
IF N2 THEN GOTO LAB31S 
IS63S SSMQS TSB3S 
GE0MGRIDCG4# G3)S 

js heading CTAXCG)# TAY(G))S- 
KS 310 ♦ (6-1)^40S 
IF NOT N1 AND NOT N2 THEN BEGIN 
Ls TDEVCG)S 

MS S0RTCTAX(G)^*2 ♦ TAYCG)442)S 
ENOS 


N= (IF N1 OR N2 THEN #75 ELSE !)♦ 

(WIX(G)^El ♦ WIY(G)^E2)S 
M3 s if N1 Or N2 then 

OECLTIME(G)^600 else TIME(G)S 
ps 100^ENTIER*M3/36000) ♦ M00CM3# 36000)/600S 
Qs 100*ENTIER^FLTIME/36000) ♦ MOOCFLTIME# 36000)/600S 
R= burn - 6RWQS 
IF N1 THEN 
BEGIN 

write ( PI »j»k#n#distance#p#o#r»grwq)s 
GOTO LAB32 
ENDS 

IF N2 THEN 
BEGIN 

WRITE(P2#ST1#N#DISTANCE#P#0#R#GRW0)S 

GOTO LAB32 
ENDS 

IF S EOL T then 

WRITE (P3#J#K#L»M#N#LE#DISTANCE#P#0»R#GRWQ) 

ELSE WRITE*P4# I# J#K#L#M#N#LE#DISTANCE#P#Q#R»GRWQ)S 


♦♦♦♦576 

♦***577 

♦♦♦♦578 

♦♦♦♦579 i 

♦♦♦♦580 

♦♦♦♦581 

♦♦♦♦582 

♦♦♦♦583 

♦♦♦♦584 

****585 

****506 

♦♦♦*587 
♦*♦♦588 
♦♦**589 
****590 
♦***591 ; 

♦***592 i 
♦***593 ; 

♦*♦*594 ' 

**♦*595 
**♦*596 ' 

♦*♦*597 ij 
**♦*598 H 
****599 
*«*4i600 
♦***601 
♦♦♦♦602 
♦♦♦♦603 
♦*♦♦604 
♦***605 
♦♦*♦606 
♦***607 
♦♦**608 
****609 
***♦610 
♦♦**611 
****612 
****613 
♦**♦614 
♦♦♦*615 


LAB32«* 


write (» 


ENDS COMMENT LINES 


PROCEDURE SPACE0PTtFB0»FB»BI*B6»RR)S INTEGER FB»RRS 
boolean FBQ»BI»B6S 

COMMENT FBQ DETERMINES WETHER THE FLI6HTPLAN 
COMPUTATION WILL BE PERFORMED BACKWARDS (FALSE) OR FOR 
WARDS t TRUE ) s 

FB DETERMINES wETHER COSTS (♦!) OR FUEL<0) OR FLIGHT 
TIME(-l) WILL be OPTIMIZED* RR IS TAKE OFF WEIGHT OR 


LANDIN6WEI6HT* 

Bi defines wether 

HORIZONTAL (FALSE) 

6G DEFINES WEThER 
OR BOUNDED ( TRUE >S 
BEGIN 

LOCAL label ITER» ENDG* CI» REPEAT» EINDS 
INTEGER AA»QUANT»II»JJ*GG»GRWP» DDG» DTG» 
boolean iteration* BCLIMB*BNQ*BNS 
INTEGER array ROW#ROWO*QUANTO(0, *F(Q) )S 


THE NAVIGATION REGIME IS FREE IN THE 
OR BOUNDED BY ONE POINT(TRUE)* 

THE CRUISING ALTITUDE IS FREE(FALSE) 


TMGS 


PROCEDURE PREP(A1»B1»AA*BB)$ INTEGER Al*BlS 
comment COMPUTATION OF DISTANCE AND FLIGHT TIME* AND 
PRINT-OUT OF the FLIGHT PLANS 
BOOLEAN AA»BBS 

COMMENT AA REFERS TO CLIMB* BB TO DESCENTS 
BEGIN 

distances DISTANCE+AIS FLTIMESS FLTIME+BIS 
LIN£(AA*BB)S burn: GRWO 
ENDS comment PREPS 


PROCEDURE editings 
BEGIN LOCAL LABEL ENDEDS 

FORMAT F03(X44* ‘CH** 14* 16* D6*2» A2>S 


PROCEDURE OPS 
BEGIN 

IF FB EOL 1 then 
WRITE(»COST*) 

ELSE IF FB EQL 0 THEN 
WRITE(*FUEL») 

ELSE IF FB EOL -1 THEN 
WRITE(»FLTIME*)S 
GOTO ENDEDS 
ENDS comment OPS 


write(pejt)s 

DETAILED FLIGHT PLAN ♦**♦)$ 
WRITEC *)S WRITE(* »)S 
IF ROUTE EOL 1111 THEN 
BEGIN 

WRITE! ‘NON SPECIAL MIN TRACK' )S 
OP 


****616 

****617 

****618 

****619 

****620 

***♦621 

****622 

****623 

****624 

****625 

****526 

****627 

****628 

****629 

****630 

****631 

****632 

****633 

****634 

****635 

****636 

****637 

****638 

****639 

****640 

****641 

****642 

****643 

****644 

****645 

****646 

****647 

****646 

****649 

****660 

****651 

****652 

****663 

****654 

****655 

***♦656 

****657 

****666 

***♦659 

****660 

****661 

****662 

****663 

****664 

****665 

***♦666 

****661 

***♦666 

♦*** 66 ^ 




eNOEO«« 


WRITEt 
•NO* HEAD 
)S 



f 


ENOS 

IF ROUTE EQL ^444 THEN 
BEGIN 

WRITE (*ATC I»)S OPS 

ENOS 

IF ROUTE EQL 5555 THEN 
BEGIN 

WRITE(»ATC II»)S OPS 
ENOS 

IF ROUTE EQL 7777 THEN 
BEGIN 

WRITE (»6RC MS QPS 
ENOS 

IF ROUTE EQL 9999 THEN 
BEGIN 

WRITE ('SPECIAL TRACK* )S QPS 
ENOS 

write ('MlN TRACK IN SPACE MS QPS 

WRITE<F03* FLUR* OATE* TT)S 
WRITEC »)S 

FL TMP TAS WINO 01 ST ACCo TIME ACCT BURN 


♦*•*67 
♦**♦67 

4i***57 

****67.. 
****678f 
****676' 
****677! 
****676 
****679 
**** 68 (H 
****681 
****6821 
****683| 

****(84 

****685 

****686 

****687[ 

**** 688 | 

****6891 

****690 

****691 

WEI6HT*****692 

***♦693 


ENOS comment editings 
PROCEDURE OQtO)S INTEGER OS 

COMMENT COMPUTATION OF SEGMENT- CONTRIBUTION. COST. 
TIME* FUEL* WEIGHT ETC.S 

BEGIN 

LOCAL label JMl. JM2» JM3* JM* ENDOS 


PROCEDURE CL(M1»M2)S REAL Ml*M2S 

COMMENT Ml DENOTES SEGMENT FUEL* M2 SEGMENT TIMES 

BEGIN 

LOCAL LABEL JUS REAL MMS 

IF bclimb then 

BEGIN 

MM= OS GOTO JU 
ENOS 

IF iteration then 

BEGIN 

MMS IF G GEO 6G THEN -80 
ELSE -lOOS GOTO JU 
ENOS 

MMS IP (0 EQL -I AND G GEQ GG) OR (0 EQL 1 AND 
G LSS GG) THEN 100 ELSE 80S 


****694 

****695 

***♦696 

****697 

****698 ) 

****699 j 

****700 ' 

****701 

****702 

****703 

****704 

****705 

****706 

****707 

****708 

****709 

****710 

****711 

****712 

****713 

•***714 

****715 

****716 

****717 


MMS M1+ 0*(GG-G)* (IF MO NEO B1 THEN MM ELSE 0)S ****718 


GRWQs GRWQ * 0*MMS ****719 

quants 9UANT*(IF FB EQL 1 THEN 1083»3*M2*15*MM ****720 

ELSE IF *'Q EOL 0 THEN MM ELSE 600*M2)S ****721 

ENDS comment CLS ****722 

****723 








ll 



JM2«« 


I 


• • 


ii 

I * JM*. 

, I :noq«. 


quant s qUANTQ(1)S 

iterations if not FBO and MQ EQL 0 AND 0 EOL -1 
then true 

ELSE false* if ROW(X) EQL 0 THEN GOTO JM» 
GRWQS IF ITERATION THEN ROW(O) ELSE ROw(I)S 
IF HQ EOL B1 OR ITERATION THEN 
BEGIN 

IF DECLFUEL(G) LSS 0 THEN GOTO JM» 

BCLIMB= TRUE* CL(0ECLFUEL<G) » DECLTlME<G) )* 
BCLIMB= false* 

ODGs DECLOlSTtG)* DT6= DECLTXME<G) *600 * 

IF ED THEN PREP<ODG» DTG» TRUE* FALSE)* 

END* 

IF NOT FBO THEN GOTO JM2S 

AI04= IF TOEV(G) LSS -10 THEN 64.0G+6 
ELSE TABlEU*G*MW)* 

IF GRWQ gEQ AI04 THEN GOTO JM* 

IF NOT FBO THEN GOTO JM3* 

GRWPS GRWQ* GRWQS GRWQ ♦ 0* AIRDIST(G)/ 2* 
(IF G EOL 1 then 4.94 ♦050000693* GRWQ ELSE 
IF G EOL 2 THEN 0.54 ♦ 050001043* GRWQ 
ELSE -0.65* 0.000117* GRWQ)* 

A AS table (0*G* SR)* GRW0= GRwP* 

IF AA LSS 0 THEN GOTO JM$ 
CL(AIRDIST(G)*10**4/AA*TIME(G)/600)* 

TMG= TIME<G)* 

IF ED then PREP(LE* TMG .FALSE* FALSE)* 

IF NOT FBQ THEN GOTO JMIS 

IF MO EQL B3 AND NOT ITERATION THEN 
BEGIN 

BCLIMB= , , , 

TRUE* CL(DECLFUEL(G)*DECLTIME (G))* 

bclimb= false* 

ODG= OECLDIST(G)* DTGS OECLTIME(G) *600 * 

IF ED THEN PREPCODG* OTG* FALSE* TRUE)* 
END* 

GOTO ENDQ* 

GRWQS 0* QUANTS 64.06*6* 

END* 

COMMENT QQ* 


ED= IF BI AND 86 THEN TRUE ELSE FALSE* 

FOR AAS 0 step 1 UNTIL F(0) DO 
begin 

ROW(AA)= ROWQ(AA)s Q* QUANTQ(AA)= 64* 04*6 


end* 

quants os bclimbs false* 

IF FBQ then 


«««*724 

*••♦725 

••♦•726 

*♦♦♦727 

•♦♦♦728 

*♦•♦729 

*♦♦♦730 

*♦♦♦731 

•♦♦♦732 

♦♦♦♦733 

*♦♦♦734 

*♦♦♦735 

•♦♦♦736 

*♦♦♦737 

*♦♦♦738 

♦♦♦♦739 

♦♦♦♦740 

*♦♦♦741 

♦♦♦*742 

♦♦♦♦743 

♦***744 

•♦♦♦745 

*♦♦♦746 

♦♦♦♦747 

♦♦♦♦748 

♦♦♦♦749 

*♦♦♦750 

••♦♦751 

♦♦♦♦752 

*♦♦♦753 

♦♦♦♦754 

*♦♦♦755 

*♦♦♦756 

♦♦♦♦757 

'♦♦♦758 

♦♦♦♦759 

•*** 7^0 

*♦♦♦761 

*♦♦♦762 

*♦♦♦763 

♦•♦*764 

*♦♦♦765 

♦♦♦♦766 

*♦♦♦767 

♦*** 7^6 

*♦♦♦769 

♦♦♦♦770 

*♦♦♦771 

*♦♦♦772 

*♦♦♦773 

♦*♦*774 

*♦♦♦775 

*♦♦♦776 

****777 



y 


ROW(0)s rRS QUANTQ<0)s 0 

END 

ELSE 

BEGIN 

ROM(F(Q))S RRS QUANT0<F<Q) )s 0 
ENOS 

IF EO THEN 

BEGIN , 

distances FLTIMEs os burns ROWCOIS 
EDITINGS MRITE(F04* RR)S 
ENOS 

IF FBQ THEN 
BEGIN 

BIS OS B2S IS 63s Q-1 

END 

ELSE 

BEGIN 

BIS Q-lS B2S -IS B3S 0 
ENOS 

COMMENT ZONE CYCLE BEGINSS 
for MQs B1 iP B2 until B3 DO 
BEGIN 

61 s IF FbO then MQ41 ELSE MQS 
G2S IF FBQ then MO ELSE Mq^IS 
M= I0*MQ4G1Q^<IF 10 EQL -1 THEN -1 ELSE 0)S 


COMMENT A CYCLE FOR GRAPH POINTS ALONG MERIDIAN 
BEGINSS 

FOR JJs sTOREICGD STEP 1 UNTIL STOREIQ(Gl) 00 
BEGIN 

63S SQ<G1*JJ)S Js JQ(61*JJ)S 

comment a cycle for graph points 

ALONG NEXT MERIDIAN BEGINSS 
FOR 11= ST0REI<62) STEP I UNTIL ST0REI0(G2) 
BEGIN 

IF NOT BI then 


DO 


begin 

IF(M GEQ 3 AND M LEO S) AND ABS 
(II-JJ) GEO 10 THEN GOTO CIS 
BNQs <FBO AND 10 EOL -1> OR (NOT FBQ AND 
10 EQL DS 

BN= (FBO and 10 EOL D OR (NOT FBQ AND 
10 EQL -1)S 
IF M LSS 3 THEN 
BEGIN 

IF(BNQ AND BK 

(M*JJ»II)) OR (BN AND 6K(M*II» 

JU)) THEN GOTO Cl 
ENDS 

IF M 6TR 8 THEN 
BEGIN 

IF (BNO AND 


••••776 1 

••••779 I 
••••780 I 

••••782 
••••783 
••••784 
••••785 
••••766 
••••787 
••••788 
••••789 
••••790 I 
••••791 
••••792 
••••793 
••••794 
••••795 
••••796 
••••797 
••••798 
••••799 
••••800 
••••801 
••••802 
••••803 * 

••••804 I 

••••805 I 
♦•••806 I 
••♦•807 

♦••♦808 I 

••••809 ! 

••••810 ! 
••••811 
•♦♦•812 
••♦•813 
••♦•814 
•••♦815 
♦•♦♦81b 
♦♦♦♦817 
•♦♦•816 
•♦•♦819 
••••820 
♦••♦821 
••••822 
•♦••823 i 
••••824 
••••825 
♦♦♦♦826 
•••♦827 
••••828 
♦♦••829 
•••♦830 
♦♦♦♦831 


BKO(M*JU*nn OR <8N AND BKQ 

Then ooto cx 

END 

END* 

Gtfs SQC62*XX)S Xs J0<62»XX)* 

XF FBQ THEN GE0M6RX0(G4*63)ELSE 
G£0HGRX0(63*6«)* 

PART2GEOMS 

XF MO NEQ 61 THEN 6GsST0REG(MQ - B2)S 

comment a cycle for flxght levels BEGXNSS 

FOR G: STOREG(MO) STEP 1 UNTIL 
STOREGOCMQ) DO 

begin 

XF NOT FBO AND MO EOL 0 THEN 
6EGXN 

TEL s OS XF G NEQ GG THEN 
GOTO ENOGS 

Nis ROW(0)S ROW<0)s 6RWQ- 
MAXTOW-20000S 

TEL s TEL +15 XF G LSS 1 

THEN 

6EGXN 

Gs GGS GRWQS os 
GOTO ENOG 
ENDS 

METPROC(G)S 0ESCLIMB<6)S 

XF DECLFUEL(G) LSS 0 THEN 
BEGXN 

TEL= OS Gs G-IS GOTO REPEAT 
ENDS 

OQ(-l)S XF GRwQ EOL 0 THEN GOTO ENOGS 
N2= ROWCX) - GRMOS IF N2 GTR 10 
AND ROW(O) EOL MAXTOW THEN GOTO 
ENOGS 

ROW(0)= ROWCO) ♦ (XF TEL GTR 10 

THEN N2/2 ELSE N2>S 

GRWQs ROM(0)S XF ABS(N2) LSS 10 

OR TEL GTR lO THEN 

BEGIN 

ROW(0)= NlS GOTO ITER 

END 

ELSE GOTO REPEATS 


END 

ELSE 

BEGXN 

METPROC(G)S 

XF MO EOL 0 OR MO EOL 0-1 THEN 
BEGXN 

GRWQ= RRS DESCLXMB(G) 

ENOS 

XF FBO THEN QO(-l) ELSE 0Q(1)S 


*4>*«632 i 
♦•••833 ; 

••**G3i| 

••••635 

•«**636 

••♦•637 

•«**G3e 

♦•♦♦839 

•••♦64U 

♦•♦♦841 

••♦•842 

♦♦••843 

**•*844 

•**•845 

♦••♦846 

•*•♦847 

••••646 

****949 

****930 

**•*951 

♦•••852 

*•**933 

••♦♦854 

♦♦••855 

*•*•936 

•*•♦857 

**««939 

♦•••859 

*•••990 

♦•••861 

••••892 

♦•••863 

••••994 

••••995 

••••896 

••••867 

••••998 

••••869 

♦•••870 

••••871 

••••872 

*•••873 

*•••974 

♦•••875 

♦•••876 

••••877 

♦•••878 

••••879 

••••990 

••••881 

•♦•*982 

♦•••883 

♦•••884 

••*•885 


IF 6RM0 EOL 0 THEN GOTO EN06 
ENOS 

IF QUANT LEQ OUANTO<J) THEN 
BEGIN 

ROw(J)s 6RW0S QUANTO (J)* QUANTS 
ROWO(U)s xs XF BI And not B6 
THEN STOREG(mQ)s GS 
IF ITERATION THEN Gs GG 
ENOS 

ENOS COMMENT FLIGhT-LEVEL CYCLES 
ENOS COMMENT ALONG-NEXT-MERIOIAN CYCLES 

ENOS comment along-current-merioian cycles 

IF BI ANO NOT BG THEN 

storegQ^mq) s ST0REG<MQ) 

ENOS comment zone CYCLES 


IF ROW(O) EQL 0 THEN GOTO EINOS 
IF BI then goto EINOS UsiF FBQ THEN F<Q» 
FOR MQS B3 STEP UNTIL BI 00 


ELSE F(0)S 


K= IF FBO then mo else MQ^IS is ROWQ(«J)S 

storeiq(mQ ♦uf b3 eql 0 Then i else oms i-f(k-i)^ 

ST0REICK)-1S JS IS ^ 

FOR MQS 0 Step i until q oo storei<mq)s storeio(mq)s 


towqs row(o)s mis if fbo Then quanto <f( on else 

OUANTQ(0)S 

LWS R0W(F(0))S XF EO ThEN 

begin 

WRITE(* *»s 

WRITE (F05* T0W0-LW>S 

WRITE (F06» IF FB EQL I THEN Ml/100 ELSE 
( (TOWO-Lw)*15^1083.3*FLTIME/600J/100>S 
ENOS 

ENOS COMMENT sPACEOPTS 


PROCEDURE FF<UU*VV)S INTEGER ARRAY UU»VVS 
BEGIN 

F(-l)s -IS for K= 0 STEP 1 UNTIL Q 00 
F(K)s F(K-1)^VV<K)-UU(K)^1 
ENOS COMMENT pFS 


PROCEOURE PR0CES<S1»S2»S3»PR»FFBb*BB»FACT0RISaTI0N)S 

integer S1»S2»s3*F^BBS REAL PRS 

BOOLEAN BB» FACTORISATIONS ^ ^ . 

COMMENT SI and S2 DETERMINE THE LIMITS OF FLIGHT 


LEVELS* 

S3 IS A DUMMY# except WHEN 


••••887 
••••886 
••••860 
••••890 
••••891 
••••892 
••••893 
••••894 
••••895 
♦•••896 
♦•••897 
♦•••898 

••••999 

••••900 
♦•••901 
♦•♦•902 
♦•••903 
••••904 
♦•••905 
••••906 
♦•••907 
♦•••908 
••••909 
♦•••910 
♦•••911 
♦•♦•912 
♦•••913 
••••914 
♦••♦9;5 ; 

••••916 I 
••••917 
♦•••918 I, 
••••919 p 
••••920 
♦•••921 
••••922 
*•••923 1 

••••924 

••••925 

••••926 

♦•••927 

••••928 

••••929 

♦•••930 

••••931 

♦•••932 

♦•••933 

♦•♦•934 

♦•♦•935 

••••936 

••••937 I 
••••938 

♦•••939 



factorisation is true# 

THEN S3 IS A DATUM FOR THE FLIGHT 
LEVEL 

IN WHICH THE M»F»P» IS COMPUTED# 

PR IS percentage trip fuel 

FFBB DETERMINES WETHER COSTS <♦!)» 

FUEL<0) OR TIME<-1> IS OPTIMIZED# 

BB«»USE LANDING WEIGHT(TRUE) #USE 
TAKE OFF WEI6HT<FALSE) 

factorisation(true) initialises 

PROCES IN THE horizontal ON THE 
BASIS OF flight TIME# 

FOLLOWED BY COMPILATION OF FLIGHT 

plan in the vertical on the basis 

OF COSTS#FUEL OR TIME* 

BEGIN 

INTEGER NN$ 

LOCAL label ENDPR#INGS 

procedure AA* 

comment specification of limits in horizontal of 

GRAPH POINTS ALONG MERIDIANS 

BEGIN 

LOCAL LABEL UTRECHTS 

STOREKO)=STOREIQ<0)=ST - A(GlQ-DS 
ST0REI<0)= STOREIQ(Q)= ST1-A(G2Q-1 )$ 

IF ROUTE GEQ 1111 THEN GOTO uTRECHTS 
FOR K= 1 STEP 1 UNTIL Q-1 DO 
BEGIN 

ST0REI(K)= 1$ STOREIQ(K)= E(K) 

ENOS 

UTRECHT.. 

FF ( STORE I » STOREI Q ) 
ends comment AAS 

procedure AAA<F1#F2)S integer Fl»F2S 

comment specification of LIMITS IN THE VERTICALS 
FOR K= 0 STEP 1 UNTIL Q-1 DO 
BEGIN 

STOREG(K) = FIS STOREGQ(K) = F2S 
ends comment AAA$ 

procedure EpS 

IF TOWQ GEQ MAXTOW OR TOWQ EQL 0 THEN 
BEGIN 

WRITE( ’SUPER CRITICAL* )$ 

GOTO eNDPRS 
ends comment EPS 

INTEGER PROCEDURE SSS 

ss= grw+reserve+pr*((grw+reserve)* 

(MINTIME*0 . 08-0 .15) -500*MINTIME-*7300+TAXI ) S 
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*♦♦♦941 

4i**«942 

4t***943 

4>«««944 
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♦***947 

****949 
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****952 

****953 
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**♦♦957 

**«4>9S6 

*»«4i959 

*«4c4e950 

**♦*961 

**«4>9g2 

**»«963 

**♦*964 

**♦♦965 

*«»«955 

♦*♦♦967 

**««958 

***«989 

***♦970 

•»«4c97X 

**♦♦972 

**♦♦973 

***4:974 

****975 

****975 

**♦*977 

**♦♦978 

****979 

♦***960 

****931 

****932 

****933 

****934 

****935 

♦♦♦*986 

****937 

****933 

♦♦♦♦989 

****990 

♦♦♦♦991 

♦♦♦♦992 

****993 


INTEGER procedure SSSS 


SSSs GRW+RESERVE^PR* t towq-lw+t ax 1 1 s 
INTEGER PROCEDURE SSSSS 

SSSS= T0WG*(GRW^RESERVE'»‘PR*TAXX)/< (1^PR>*Nn-PR*TOWQ)S 

PROCEDURE NNQ<A)S VALUE AS INTEGER AS 
NNS IF NOT BB THEN TOW ELSE AS 

PROCEDURE TwS 

COMMENT safeguarding AGAINST OVERLOADINGS 
IF TOWQ GEO MAXTOW OR TOWQ EQL 0 THEN 
BEGIN 

WRITE(PEJT»S WRITE(»OECR LANDING W»)S 
GOTO ENOPR 
ENDS 


COMMENT IF FACTORIZATION IS TRUE* THEN THE OPTIMIZATION 
TAKES PLACE FIRST IN THE HORIZONTAL AND THEN 
IS FOLLOWED BY IN THE VERTICALS 
IF FACTORISATION THEN 
BEGIN 

AAS AAA(S3»S3)S NNQ(SS)S R0UTE= llllS 
SPACEOPTdF BB THEN FALSE ELSE TRUE* -I* FALSE* 
TRUE*NN)S TVS 
AAA(S1*S2)S GOTO ING 

END 

ELSE 

BEGIN 

AAS AAA(Sl*S2)S NNQ(SS)S 

SPACEOPTdF BB THEN FALSE ELSE TRUE *FFBb» FALSE* 

FALSE *NN)S TWS 

EPS NNO(SSS)S 

SPACEOPT tIF BB THEN FALSE ELSE TRUE*FFBB» TRUE* 

FALSE *NN)S TWS 
EPS NNQ(SSSS)S 

SPACEOPT ( TRUE ♦ FFBB ♦ TRUE * TRUE * NN ) S EPS 
IF FACTORISATION THEN ROUTE= 1000 
ends 

ENOS COMMENT PROCESS 

IF ROUTE GEO Ull THEN 
begin 

REAO(PCF(*NASAD»). ST* STl* DONE* ERR)S 

610= 20NEHsT)S 

FOR K= I STEP I UNTIL 0-1 00 

BEGIN REA0<PCF< *NASAD* ) * SToReKK)* DONE* ERR)S 
WRITE(F02* 1* ST0REI<K))S 
STOREI(K)= sTOREIQ(K)= STOREI (KI-ACGIQ-^XO^K-DS 
ENOS 

620 = ZONE I< STDS 
WRITE(F02* 2* ST*STDS 


**«*99i| 

••♦•995 

♦•••996 

•♦••99T 

*»*»99e 

« 4>«*999 
♦••♦1000 
••••lool 
•♦♦♦1002 
••♦•1003 
••••1004 
♦•♦♦1005 
♦♦♦♦1006 
♦♦♦•lOOT 
♦♦♦♦1008 
♦♦♦♦1009 
♦♦♦♦1010 
♦♦♦♦1011 
♦♦♦♦1012 
♦♦♦♦1013 I 
♦♦♦♦1014 
♦♦♦♦1015 
♦♦♦♦1016 
♦♦♦♦1017 
♦♦♦♦1018 
♦♦♦♦1019 
♦♦♦♦1020 
♦♦♦♦1021 
♦♦*♦1022 
♦♦♦♦1023 
♦♦♦♦1024 
♦♦♦♦1025 
♦♦♦♦1026 
♦♦♦♦1027 
♦♦♦♦1028 
♦♦♦♦1029 
♦♦♦♦1030 
♦♦♦♦1031 
♦♦♦♦1032 
♦♦♦♦1033 
♦♦♦♦1034 
♦♦*♦1035 
♦♦♦♦1036 
♦♦♦♦1037 
♦*♦♦1038 
♦♦♦♦1039 
♦♦♦♦1040 
♦♦♦♦1041 
♦♦♦♦1042 
♦♦♦♦1043 
♦♦♦♦1044 
♦♦♦♦1045 
♦♦♦♦1046 
♦♦♦♦1047 


ENOS 

6E0MGRID<ST»ST1>S MINTIMEs 
6E001ST(E3* (LA1'*‘LA2)/2)/<465S CTOS 
E(0)s E<Q)= IS FOR Ks I STEP 1 UNTIL Q-1 DO 
E(K)s A( 61 Q^X 0 «K)«A(G 1 Q^I 0 *K* 1 )S 


I 


COMMENT 

If WHEN A flight PLAN 
CALL *PROCES(l* 3» 1* 
2. WHEN A flight PLAN 
CALL *PROCESU» 3* I* 
3f WHEN A flight PLAN 
CALL «PROCES(lf 3* 1* 


FOR THE OPTIMAL COST TRACK IS NEEDED 
0f03» 1* true* FALSE)* 

FOR THE OPTIMAL FUEL TRACK NEEDED 
0f03* 0» TRUE* FALSE)* 

FOR THE least TIME TRACK NEEDED 
0.03*-l* true* FALSE) *S 


*♦•*1048 

****1049 

****1050 

*♦*♦1051 

****1052 

****1053 

****1054 

***•1055 

***♦1056 

****1057 

****1058 

****1059 

****1060 

****1061 


|l 



PROCES(l» 3f If 0f03»-l* TRUE* FALSE)S 

****1062 


PROCESd* 3* If 0f03» 0* TRUE* FALSE)S 

****1063 



****1064 


GOTO RUN 

♦***1065 


ENDS 

***♦1066 

OONEf f 

ends 

****1067 


WRITE(« »)$ WRITE(* *)S 

****1068 


WRITE(» *** END OF RUN ***»)S 

****1069 


GO TO EOPS 

****1070 

ERR«f 


****1071 


WRITE(* ♦** INPUT data ERROR ***«)S 

****1072 

EOPf. 


****1073 


^ ELT NWS»1»810227»039147 
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G kIMERA L. 

Thl« prugram which wa^ writian for e upar l.ncnt* uf thr^a 0 iiT)«h% lunal 
flight planning for DC8-D3 rM*craft* ov*r the North Atlantic hat 
been converted for use in UNI VAC 1108 in ALGOL 60 and will be 
enhanced to be appropriate for flight planning for OClO and B747 
aircraf c». 

Reference i» made to pages and formulas in 'Optimal Traci Selection 
and 3-Dimensional Flight Planning' written bg H. N. DE JONOi KNNl No. 93. 

INPUT /Q UTPUT DATA D ESCRIPTION 

There are six data tables for compressed meteorological data ( geopotential 

iieight(m gpdamiand temperaturelcentigrade degree) ). 

. ^ 

See p 126. Fig. 39. 

Oeopotenliai values and temperature are depicted from the grid points in a 
Cartesian grid superposed on a polar stereographir chart projection 
with standard parallel at 60 degree N. 

This grid i* part of the well-known octagonal grid prepared by 
fhe Ncttionai rie teoro 1 ug i ca 1 Center. 

The y-a«is runs parallel with the 80 degree U meridian. 

Tiie are.* covered is rectangle cumprisinq i£ t 16 grid points 
Fee 0 iit|>uL lilting. 

Those sii tables show meteorological date leach table 12 « 161 for 20C. 
2t/0* and 30'J inbar for two validity times 12 hours apart. 
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Table l Stiout specific range for flight levels 310» 3bG< and 390 
(See p.116. table 1). This is an important economic indei uhich 
is tiid air distance covered bg a turbo-jet aircraft per unit of 
fuol c«>nt.uinption 

This datum is almost independent of temperature and depends on 
the aircraft uieiqht onlg at a prefixed flight level. 

Table 2 shouis maximum weightCkg) per flight levels 310> 3t70 and 390 
(see p. 122« table 9). The maiimallg allowable weight depends 
on flight altitude and temperature deviation from standard. 

TABlE 3 shows climb distance (nautical mile) as a function of weight (kg> 
ami flight levels 310> 350> and 390. 

See H 1 IHf tjble 5 . ^ 

Table a snows climb time (1/1000 hrs) as a function of weight (kg) and 
flight level*;. bee p. 117. table 3 

'Table S shows cilmb fuel (kg) as a function of weight and flight levels 
See p. 1 IV, table 7 

Table 6 stiowa Climb time temperature correction as temperafure deviation 
and climb time (1/iOOC hrs) at standard temperature 
Sci* V. Ii7. tabli’ M. 

Table 7 siiows Climb fuel temperature correction (kg; See p llv, table f 

Table 8 s»iOu>s climb distance temperator# correction <n m). 

Sei' p llH. toble 6 


I MITI AL I 21MC INPUT PAhAflE rERS ’ 
ElUP sequi-nse number 


'^A)( I ■ taii fuel (k<|) 


leio fuel weight (kg). 
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RESERVE = rest*rve fuel (ky) 

TUW - -1 computaion starting in end-point of flight. 
ilAXTOiJ - maximum take-off uieiyht. 

H^XLW ~ maximum landing uieight. 

DATE ■= fiate/month/gear. 

10 = -1; east bound fiight> li uiest bound flight. 

TfT = time of departure. 

ROUIE = T te indicator. 

Sr = starting graph point number. 

STl = destination graph point number. 

Uote here that all input data are read from card images of PCF file 
FCF file name is NASAD. 

i .UTPUT 

In the flight plan (see output listing). 

MO. = graph point number. 

HEAD = heading in flight. 

Fl = flight level m 100 ft. 

IMF = off standard temperature in centigrade degresb 
TaS = true air speed 

WIND - kail (head) uind along the track 
DIST - distance of flight segment 
ACCD = accumulated distance fioum(n m) 

TIME = time along segment floum 

Acer = accumulated time floum (hr. min). 

B(JRN - +^uel consumed in segment (tgi. 

WEI(»Hi - me I gilt (kg) 

TOC = top of climb 
TOD = t.ip Of descent 
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i\ XlKA Cl ihMErJTB ifilH. ''He PRiiOftAri . 

4* - •»? 

ProcuOure INTERPOL computes u iinedr interpolation. 

51 - a9 

Procedure SUlil end SOIi2 compute summations. 

These pf'ocedures are used for computing average values in 
pvocedui-e METPROC. 

72 - 73 

Procedure WW is a write-procedure. 

76 - 77 

pl> )i2, and p3 are flight levels. 

qlf q<L'. and q3 are pressure altitude for standard pressure levels 
300> 250 » and 200 mbar. 

79 - B4 • 

reads in compressed values for geopotential height and temperature in 
200> 250> and 300 mbar> valid for two standard times# 12 hours apart 
and fills array HT with tiiose data, 
eb - -0 

writoi: out tfie inputed data for checking purposes 
92 - 11 2 

produces tiie meteoroiog icai pai'ainvters ^compressed geopotential and 
tempera tore data, which apply in grid points at flight levels 3lO. 350 
<Mid 3v0 Therefore the inputed data for the standard pressure altitude 
are di's troy ud 
315-118 

The graph points are labelled by numbers running fTom 0 to ll5 
iriote that the graph point number shown in Fig. 37. F llC runs from 


1 to lift). 


$ 
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'fie greiph is such that the point sets consist ot socsets of* points 
.arranged along ^meridians' (index running from O to 

Arrag Aiindexi contains the graph point number assigned to the most 
nurtherlg point on each murid ian. 

Array A is used for the computation of zone index (procedure ZONE I > 

119 - 133 

Array V has values assigned for geographic coordinates for graph points 
and check points on continents : latitude (4 digits) and longitude ( 4 

digits) in degrees and hundredths. See Fig. 37> P. llO. 

135 - 142 

fills Boolean array BK and BKQ uiith 'TRUE' : ( initialization ). 

143 - 156 

denotes blocking in continental airways (connectivity in the graph). 
Boolean array entries of BK and BKQ ind«>cate which graph points are 
connected or blocked ( TRUE=blocked and FALSE=connec ted ) . 

Array BK is used for the zone index 3 and array BkG is used for 
tlie zone index > 8. 

Bk ( a. l> I c ) ^FaLSE means point'b' on meridian 'a' is connected with 
puint 'c' on mef idian 'a+l ' 

Were point 'b' denotes the 'b ' th point from the bottom on the 
meridian. Example; see Fig. 37 Bk( 2«2>2) = FALSE means graph 

'•oint number 3 on meridian 2 is connected with graph point number 5 
Ot: meritli.in 3 
1.^3 - 1 .'2 

Procedure ZOnEI determines tiie index of the zone associated with a 
rraph point number. r'a name ter vJ = graph point number 
-or example, a yrapfi point number 49 gives tlNEI(49)=;6 
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PruCbOure uT dissetts latitudis and longitude from the compressed 
c uord in.^t<;s 

lALAT means latitude and LAlONG means longitude. 

Parameter K is an index of arraij V. 

1B5 - 216 

Procedure LlS determines actual latitude and longitude 
value for a graph point, (see procedure DT). 

The sign of longitude is changed in east of Greenuihich. 

Parameter U-graph point number. 

218 - 2rJ7 

Procedure CTCl prepares a time i<i!>cant arrag DDi zonallg specified* 
to be used 'Hor the time interpolation later on. 

A=0. 4 denotes an estimated flight time of 24 min. per continental zone. 
There are lb instant values <index 0 - #iidex 14). 

239 - 256 

Procedure HH compuces grid point values for geopotential or temperature 
( CH01Ct=TRUE computes geopotential value, ujhile CHQICE*FALSE computes 
temperature value) in flight level index F (index 1 denotes 31000 ft 
and indices 2 and 3 denote 3500u and 3*^000 ft respectively). 

Parameters V and u) denote index numbers of array HT ( meteoi o log i ca 1 > 
Interpoia ciun takes piate according tu z one-ass i gned time instants Df. ^ J. 
258 - 2. '5 

Procedure GEuP computes qeopotential (when CmQ 1C£~ fROE^ or temperature 
(when Cri01 CE=P AlSE ) value in arbitrary point, using adjacent 
grid point values from array HT (see procedure HH> ana bi-linear 
interpo latiun 

Parameters V and mill be modified to denote index numbers of array 
H.r Parameter C = flight level index 


2/7 - 2:-!4 


7 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Procedure XYTRANSF converts the geogr^ph ical coordinates (LAT.lDN) of 
a point on earth into the Cartisan coordinates (Y.X) 

See P n2—63. 

296 - 299 

Procedure GEODIST computes distance on the sphere from Euclidian 
distance DDD given in units of the rectangular grid on the map. 

LAT denotes latitude. 

291 - 3ul 

Procedure GEOMGRID computes Euclidian distance ana unit vector 
pointing along an arc in the graph Petueen two arbitrary graph points 
u/ith respect to rectangular grid system. 

PP and GG denote graph point numbers of end points. 

E3 denotes Eucledean distance • 

El and E2 denote unit vectors for i» y components respectively. 

303 - 3ih3 

Procedure PART2G£0r1 determines length of segment by summation of 
contr 1 but 1 ons from subsegments (LENGTH! I )i. 

These subsegments ai e identified by intersection of an arc with 
coordinuti^ lines in Euclidean grid. 

Normaliied coordinates (running from 0 to i > are stored in array AB 
’Jiese values serve also as weighting factors later or. 
a denotes distance computed for averaged degree of latitude on the 
sphere Arrag LENGTH! I ) contains lengths of sob segmeiiis on the 
sphere AIu 2 denotes interpolated degree of latitude. 

£ee P o4— Bf; 

3a5 ~ 4 1 0 

Procedure flETPROC processes meteorological paramaters. 


ORIGINAL PACr: 

OF POOR QUAUrv 

PeirametKi' S^fliqht level indei! (1.2 or 3) 

ChTAli.S 
270 - 3.1 

coiTiputes geopotential value at end-points of a flight segment. 

372 

computes drift angle bg using bellamg's formula (see P. 78). 

376 - 3h1 

estimates coordinates in mid point P of segment. P^lXster. Ysteri. 

334 

estimates temperature in mid point P at flight level S. 

Temperature has minus sign. 

337 - 3v2 

wind components, geostrophic wind Formulas snown in P.69-70 is not 
used (reason: not tu smooth too drastic-a41g from geopotential fields, 

and to economize in speed of computation pg computer). 

See formulas 5. 22 and 5. 23 in p 66. 

Array WNDX and (JNDY contain K-component and g-component of wind resp. 

393 

computes true air speed (TaS) for the specific type or aircraft. 

For values of cunstants see Smithsonian tables. Meteorological text OooK 
or ICaO Standard Attnosphere. 

394 - 395 

computes flight time cumponents using vector sum of wino and TaS. 
and segment length (LrfJGTH (I>) 

398 - 40? 

determines average values for meteorological ana other parameters alt 
segment under consideration using previouslg calculated contributions 
TEM=temperature. TDEV(S)"temperatura deviation 
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Procedure REaDO reads in performance data and proceoure URITEQ directly 
urites out those inputed data for checking purpose. 

499 - 500 


coefficients used in algorithms. 

502 - 507 • - 

reads in and writes out initializing input paratnecers for flight plan. 
508 - 5J0 

In case of prescribed track flight# the routine reads in Q (number of 
graph points on the prescribed route) 

In another case# it reads in ST(starting point number) and ST 1 ( ending 
point number) and decides Q value 

In the latter# (i=( nuir#ber of zones between source and destination) - i. 


5:29 - 5.-:2 

Procedure SO and J'J are conversion formulas for indices used in 
subsegment aloorithms (bg using entiles of array A(inoex) or 
array P*( i ndex > ) 

5:i4 - 5-.;5 

Procedure oESCLlnE computes performance by using the table look-uc 
f.roceduie in cl\mb(if rlu=0) or descenc(if Pi(l not* O' 


10 


ORiGtNAi pA,-- in 


AiTdy Di-.ClTiME gets cJimb time or deatent timo. 

Arrag CnCLDlST gets climb distance or descent distance. 

Arrdig Dit.ClFOEl gets climb fuel or descent fuel. 

Array AlRDIST gets climo air distance or descent air distance. 
Array TIME yets climb time or descent time. 

Note that in the descent zonei performance data table is not used. 

specific expressions by company's planning policy are chosen. 
See formula 6. 3l in P 120-121. 

5S>7 - 5ii9 


Procedure ARCTANG is for computatiun of heading degree by using 
arc tanifeiit formula 
5fel - 5.^4 


Procedure lO is for computation of longitude from Euclidean coordinates 
in rectangular grid See P. 62< formula <&. 17 and 5. l8. 

Zl and Z2 denote Euclidean coordinates. 


5a6 - 574 


Procedure HEADING computes heading in degree as required for flight 
plannintj purposes 

LI and i_2 denote x— component and y— component of true air speed. 

576 - 617 

procedure LINE uiriccs out optimized flight plan 
See the output listing of the flight plan 
619 - 9. '6 

Procedure SPaCEOkT is a keg procedure tor optimal i z<«tion. 

In this procedure crucial is that some parameters mag not surpass uppe'" 
bounds, for example aircraft uieight 

This requires the inclusion of several protective statements in order to 
ensure proper functioriing of the operational scheme in-flight 


1 1 ORIG‘NAL PAGE : - 

OF POOP 0U»‘'I.ir:' 

FjiQ determines uihettier the flight pion computat lOn u/ii] be perfurmod 
bac kUriT d ( f a 1 6c ) or forward ( true ) . 

Fb determines whether cost(+l> or fuel(u) or flight time<-l) will be 
Optimized 

Cl determines whether the navigation regime is free in the horizontal 
(false) or bounded by one point(true). 

Cti determines whether the cruising altitude is free(false) or 
bounded ( true > . 

RR is take-off weight or landing weight. 

DfeTAlLS. 

637 - 6-15 

Procedure PR£P computes distance and flight time> and prints out one 
line uf flight plan 

Alssegment distance. Qlssegment flight 'time. 

AA=true denotes 'climb' and BQ-true means 'descent' 

See the output listing of the flight plan. 

617 - 6V4 

Procedure EDITING prints out proper heading lines 
oV6 “ 7r.8 

procedure uQ computes segment contribution! fueli etc. 
bJhen o“-l computation i» for in-flight direction. 

i 

When computation is for opposite flight direction. f 

7U2 -7-2 i 

j 

pToceOure CL coon-otes the aircraft we i gh t ( GRidG > ana total tueli flight \ 
time or custs at end-point of a flown segment, (the contribution of fuel 
in sti-*P'. included). 

farameter Ml denotes segment foel(kg) and M2 segment time(min) 

In the climb phase no step contriDutiun is taken into account as 
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r«(\un-cil (uCi.lMB-'THUh*) 

Cu d«i)Ot«b tligiit ItfVfe'l index At the gr«pn point when the Cutnputetiun 
itarCe. & denotes the flight level index to be investigated, 
step c 0115 oiTt|it ion: UH»100 kg/4000ft 

DaWN=-80 kg /4000ft 
(surplus of segment fuel) 

ITERATION becomes TRUE when (1) bac k—trac k ing (2) first zone is reached 
(3) TOW is unknown and Wb is known (see P. 120). 

In arrag ROW the aircraft weight is stored. 

GRWQ denotes aircraft weight at the end point of segment. 

GOANT denotes cost* fuel or time at the end point of flown segment. 

729 - 7:»6 

computation in climb or descent bg procedure Cl. 

If ED is TRUE then it prints out the results. 

739 - 741 

safeguarding statements. 

744 - 752 

computt-s segment contribution of changes in aircraft weight 
Note that fuel consumptiofi is eipressed as loss of aircraft weight 
The table fur specific range iP 116) is used 

Formula 6.27 (P 115) requires an estimate of weight halfwag a flight 
segment Fur that put pose coefficients as shown in table 2i P 116 
required «rl. r2) 

753 

jumping back to label JNl is for safeguarding that the actuallg deiiveu 
aircraft weight does not surpass the maximallg allowed weight 
755 - 7s2 

computes t)ie contribution for tne last segment 




ORiGfNAL PAGE 13 

If £D iy. fRUE tlien it prints out the results. POOR QUALITY 

765 

CrRMQ and QUANT attain absurd values for safeguarding operational 
perform. mce 

771 

ED becomes TRUE uhen all degrees of freedom in tne horizontal and the 
vertical are lost* for example when a call of procedure SPACEOPT takes 
place in the final mode for flight plan computation along the (optimal) 
track found in previous calls. 

(when navigation regimes is bounded bg one point and cruising altitude 
is bounded). 

772 - 7^7 

Initial ization and preparation of all parameters needed for the 
algorithmic process of optimization. . 

PuQ«THUE means flight plan computation will be performed forwards* 
while FliQ-FALSE backwards 

Arrag F is filled in line 92S-932 (see procedure FF ) 

The process can be activated for arbitrag begin- and end-points. 

In order to reduce storage space a renumbering is made for all subset 
points So arrag F is defined and functions completeig analogous to 
arrag A for whole graph point set 

Array SlOHEl and STOREIG contain the assigned nofuber* of labelled 
•?nd-pnii>ts on meridians in the subgraph 
"line cycle runs from line 7v9 to vOa. 

Tne zone cycle repeats as many times as the nun^uer of zones. 

A cycle for graph points along meridian rons from line 606 to 9oC. 

This cycle repeats aS many times as the number of graph points in 
rtie current meridian 
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OF POOR QCALITY 

A cgci« fur gr^ph points dlong udxt it<vriili4n runt from lin« 611 to 6vB. 
This cgclv repeats es m«ng times es the number of gr.iph points in 
ttie ni'Kt meriiiiw^n. 

A cgcie fur flight levels runs from line tiA2 to 69m 

This cgcie repeats as mang times as the number of flight levels used. 
810 

conversion for Indices used in Subsequent algorithms. 

815 - 835 

If the navigation is free <B1>FALS£) in the horiiontal* it checks 
u'hcther a segment (II« JJ) is blocked or not. 

BK or BkQ * TRUE means airuiag is blocked. 

If blocked* it skips computation for the segment. 

836 - 8-»0 

conversion for indices (point (>4>. . ^ 

determines geographic elements for a segment between points 03 and U4. 
See procedures (2E0M0R1D and PAR(2QEUIi. 

846 - 875 

In case of backtracking in a flight plan computation an iteration 
procet-s i*. put into action in the climb zone. 

This IS dune in order to determine the (unknouini taM-*off weight bg 
iteratiun In such a wag that the climb parameters do match with the 
parameters found during backtracking when arriving ir. tne climb zone 
The iteratiun starts with a take-off weight 2000C kg below the 
(known) maKimum take-uff weight. 

Arreg ROW contains the (optimal) value of aircraft weight found during 
optimal 1 zatiorw using a zune-cgcle. 

Tlie iteration finishes when either the weight difference at the point 
Of matching is less than lO kg or the number of iteration steps 
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oicttdfi 10 (poor convrrqoocr) 

In toch stop tak»-off wvight is adjusted* s»o lino fc(6‘^-87S. 

8/9 - 8m7 

similar computations but not backtracking cast or not take-off 
weight .idjustment case. 

698 - 8v4 

These statements are crucial optimaliiation criteria. 

909 - 9;>5 

Preparation of elements which are required for a subsequent call for 
procedure SkACEUPT with lessened number of degrees of freedom 
including a final call for computation and presentation of flight 
plan data. In this final stage no degrees of freedom are left. 

Procedure SfACEUPT then merelg operates along a predefined track 
(the solution of previous calls). . ^ 

9:?8 - 9:^2 

Procedure FF determines a zone index arrag P for a subset of graph 
points analogous to arrag A. (arrag F is a sub-zone arrag). 

934 - 103S 

A tail for prcceuure PHOCES results in the production of a flight plan 
along whatever route is desired 
details 
9.‘s2 - 9 3 

Spec 1 f 1 cat 1 on of limits in hurizoncai of graph points along mtTidians 
Arrag SiOkEI and arrag STuREIU store lower limits atto upper limits 
respectivelg for the number of graph points on each meridian. 

If route >■' 1111. then production of flight plan alur.g a presribed 
r uute 

In this Case the procedure ttPACtOKT works in a degenerate mode 
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■975 - 9hO 

Spec j P 1 c»jt ) on of limits in tne vi»rtii.ai, 

Array SlOKE^i and array STURECU store lower limits and opptrr limits 
respectivelg tor the number of flight levels on each meridian 
982 - 9k7 

Procedure tP indicates that the flight planning computation blows up 
989 - lOOO 

See P. 123-124 for procedures SS» SBS and SSSS 
These procedures compute landing weights. 

NN contains landing weight or take-off weight. 

1002 - lOoB 

PruCedore I'W is for safeguarding against overloading. 

1011 - 1035 

activates process ns follow. 

If FACTORIZaIuN is TRUE, then the op t ima 1 i rat ion tatet place first m the 
horizontal and is followed by an op t imai i zaion in the vertical. 

The following steps occur 

(a) optimal 1 ration in ttio horizontal ba«-ed on time 

(b) optimal 1 zation m the vertical, using the trad solution found m 
(a) and based on fuvl< tiifie or cos tv 

(c) computation of tlignt plan along voiotion found in <D). 

If the optimaliiation takes place in fr^e space, thv sequence of 
results in 

fi) Op 1 1 mal 1 za 1 1 on in Spaie bast'd on time, tuel or costs 
iii) op limal i ration in the vertical throoqh honzonai track found in 
:i>. b^sed on time, fuel or costs 
dll) compilation flight plan along trick solution found m *ii) 

Note that 111) could U€ bypassed as the solution is ^.jun.j aJ'eody 
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;n (i). BuC (ii) cnn be generaccd uiitti a slightlg different i<andinq 
weight. The compilation of a flight plan along a prescribed route 
passes through all three procedure calls. 

This means that in fact the computation is repeated threefoid< however 
with properly tuned landing or taKe-off weights. 

In order to protect against subsequent calls of procedure PROCES> RUUrE 
is assiuned by 1000 
1037 - 1043 

Read staienient in case that a flight plan compilation is desired along 
a prescribed track. 

This tract is specified by graph points indicated by their numbers 
1049 - 1052 

Computation of an estimate of flight time to be used for tlie estimation 
of other parameters • 

The same with distance. 

Procedure C I'Q makes time instant array DD for each lOne (mer id lan ) . 

Array E is for stor?nq the number of graph points used on a meridian 
1062 - 1 vO,4 5 

perfuTiiis a process for the prouuction of a flight plan. 

Actual pa r ame te r s when calling procedure < 
use landing weight iBm=TRUE). 
cost op l imai i la t i on <hF>iB=l). 

FrBB=0 :iieans fuel opt. FFBB=-1 means time opt. 
regularity percentayo (PR-0 03) 

vligh levels used between Sl=l for 31000 and S2-3 for 29000. 

S3=li tiie value is immaterial if PACTuRIZaTION is FaLSE. 

For requited flight plsns ; 

i. when a fiii,ht plan for the optimal cost traci< is needed. 




call 'pruCfti ( 1 / 3. 1 » 0 03. i I ThUh!. FALSt ) ' . 

2. when a flight iilan for the optimal fuel track is needed, 
cail 'proces ( 1 . 3. 1 . 0. 03. 0. TrtUfc. FALSt ) ' . 

3 . when a flight plan for the least time track is needed, 
call 'proces ( 1. 3. 1. 0. 03. -J . TRUE. FALSE) 

For the flight plan 1. 2 or 3 above. FACTOR I ZAT10 N=FAlSE. 
For other flight planning simulations. FACTORlZATlOr\|=TRL'E. 
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